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Abstract 


A description is given of experiments carried out during 1949 and early 1950 that were 
designed to produce some phenomenon exhibiting effects of the resultant horizontal Coriolis 
force on a symmetrical vortex that results from the latitude variation of the Coriolis force. 
A circular cap of dense liquid at the base of a small hemispherical shell filled with liquid was 
forced to rotate about a vertical axis independently of the rate of rotation of the remaining 
fluid and containing vessels. Within the ranges of 20 to 150 r.p.m. for the vessel and o to 120 
relative r.p.m. for the cap vortex, the center of mass of the polar cap never showed appreci- 
able displacements away from the polar axis exept for a range of sufficiently high relative 
rotations when the vortex was anticyclonic in sense. Numerical checks are given with a simpli- 
fied calculation in spherical coordinates of the torque equilibrium for the vortex as a whole. 
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1. Introduction 


In two recent papers Prof. C.-G. RossBy 
(1948, 1949) has advanced the idea that a 
certain type of instability of anticyclonic vor- 
tices may have important dynamic and syn- 
optic applications. He obtains a resultant force 
(or better a resultant torque) acting on various 
vortex models in a direction away from the 
pole as a consequence of only somewhat 
restricted hypotheses concerning the pressure- 
wind field relation. The same computation, 
though with some differences in results, has 
been carried out by Davies (1948). Such an 
effect, if not obscured by some of the processes 
not fully taken care of, would have many 
consequences for the pole-to-equator transfer 
of air within the anticyclones and for related 
mechanisms in the general circulation. 

In view both of the importance of the implica- 
tions of this result and of the delicacy of con- 
clusions which depend on finer aspects of the 


1 This research has been supported in part by the 
Geophysical Research Directorate. Air Force Cambridge 
Research Laboratories, Air Materiel Command, U. S. 


Air Force. 


pressure-wind relation, Prof. Rossby suggested 
in the spring of 1949 that an attempt be made 
to get experimental observations of some 
phenomenon as nearly analogous as possible. 
This has been done by making modifications 
of the rotating hemispherical shell equipment 
used for an earlier investigation (FULTZ 1949) 
and the following constitutes a preliminary 
report on results. 

A little reflection soon shows that it would 
be impossible to set up an experimental baro- 
tropic vortex with complete dynamic and 
geometric similarity to the case contemplated 
by Rossby. It was therefore decided to attempt 
to use a two-fluid system in the available appa- 
ratus and modify the theoretical analysis to 
correspond. A lower circular cap of denser 
fluid (shown in cross-section in fig. 1) is caused 
to rotate about the vertical axis of the system 
by a fan which is driven independently of the 
rotating containing shells. The lower cap is 
overlain by water that fills the hemispherical 
space to the equator. The fan is thus able to 
produce a relative rotational motion of the 
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Fig. 1. Cross-section through the hemispherical shell 
showing the approximate position of the cap of denser 
liquid. 


lower cap of fluid with respect to the upper 
fluid and the hemispherical bowls. The two- 
fluid arrangement was adopted in order to help 
confine the relative motion to the lower cap 
and to assist in identifying the cap. 

The simplest case calculated by Rosspy 
(1949) was that of a circular vortex of incom- 
pressible fluid of unit depth and radius R 
rotating rigidly with angular velocity w in 
a resting environment. He found for a vortex 
of infinitesimal size centered near the pole, 
that the displacement b of the center from the 
pole would obey approximately an equation 
of the form! 


ah Re 
= (Lrok)t (x) 


where fy is the Coriolis parameter at the pole 
and a is the radius of the earth or sphere. The 
term on the right is an average force per unit 
mass on the vortex that results from the varia- 
tion of the Coriolis parameter. It tends to 
displace the vortex as a whole either toward 
or away from the pole along a meridian. 
Equation (1) is seen to imply that a cyclonic 
vortex ( > 0) will undergo stable oscillations 
across the pole while an anticyclonic vortex 
(wo < 0) will experience an exponential flight 
away. 
In the experiment described there are several 


1 Among other approximations to be discussed later, 
it may immediately be noticed that, since the fluid is 
frictionless, the assumption of solid rotation within the 
vortex is incompatible with the vorticity equation: f + & 
=individual constant, except when the center is located 
at the pole itself. Thus the real problem being considered 
in this case is that of a rotating solid disc in a frictionless 


fluid. 
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factors that do not correspond with Rossby’s 
calculation where the Coriolis force makes the 
only contribution to be calculated: 

1) The two-fluid system of lighter above 
denser liquid introduces a hydrostatic stability 
not present in the atmospheric vortex. 

2) In terms of the equations of motion 
relative to the rotating bowls the centrifugal 
force terms do not drop out as they do in 
the ordinary meteorological equations. 

Both of these differences could theoretically 
be eliminated by setting up a vortex relative to 
a rotating paraboloidal shell of fluid that is 
rotating at such a rate that the sum of the 
horizontal components of the gravity and 
centrifugal forces vanish. We contemplate 
carrying this out with apparatus which is be- 
ing designed at present, even though the 
problems of generating a more or less isolat- 
ed vortex will be more difficult to overcome. 

Qualitatively in the present instance, we 
attempted to find whether there was any 
difference between cyclonic vortices and suffi- 
ciently strong anticyclonic ones of the type 
to be expected from the theory: namely, any 
tendency for the center of mass of the experi- 
mental vortex to be thrown away from the 
pole beyond some critical anticyclonic rotation 
value. Such a phenomenon was actually ob- 
served over certain ranges of conditions, with 
some accompanying difficulties of interpreta- 
tion which are considered at some length 
below. 


2. Theoretical considerations 


In the first attempts at making an analysis 
corresponding to Rossby’s for this experimen- 
tal case the centrifugal term was neglected. 
The gravity term was estimated as approxi- 
mately that for a simple pendulum of length a 


with gravity reduced by the factor y = Ore 


where o is the density of the lower fluid and 
0’ that of the upper fluid. (r) then becomes 
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where g is the acceleration of gravity and the 

other terms have been defined previously. 
However, in order to get a reasonable sized 

experimental vortex it was necessary to carry 


ee‘ 


EXPERIMENTAL STUDIES OF A POLAR VORTEX I 139 


the lower cap up to as far as 40° latitude. This 
made it obvious that some attempt would have 
to be made to perform the calculations without 
restriction to infinitesimal size and also if 
possible without restriction to infinitesimal 
displacements of the center from the pole since 
only rather pronounced ones are visible. At 
first we tried analytical integration without 
approximation of RossBy’s equations (1949, 
p. 170). These are written for a transforma- 
tion plane tangent at the pole with projection 
from the center of the sphere. These approaches 
were only partially successful and Mr R. R. 
Long of our laboratory, going further, dis- 
covered that an exact integration of these 
terms could be carried out for the sphere 
under any of several assumptions as to the 
velocity field in the vortex. The remaining 
uncertainties in this calculation are mainly of 
a kind associated with what Rossby calls the 
“escape motions” of the environment fluid. 

The above mentioned analysis by Mr Long 
may be put in the following form: consider a 
circular vortex of incompressible fluid and 
of unit depth centered at colatitude « of a 
rotating spherical shell. This vortex is to be 
imbedded in a layer of incompressible fluid 
of density o’ which temporarily we consider 
at rest relative to the shell. Choose coordinates 
as in fig. 2 where 


are rectangular coordinates with z in 
the direction of gravity and also that 
of the angular velocity © of the 
spherical shell. 

x’,y',z’ are auxiliary rectangular coordinates 
as shown with z’ passing through 
the center of the vortex. 

are spherical polar coordinates having 
the usual relation to the x, y, z frame 
(© colatitude). 

r’,@’,A' are a similar set corresponding to 
the F7 bare: 


2% 


r, O, À 


a is the colatitude of the center of the 
vortex in the r, ©, A system. 
0,’ is the angular radius of the vortex. 


The essential point noticed by Mr Long is 
that for displacements of the vortex as a 
whole the constraints are such that the appro- 
priate equations of motion are the torque- 
angular acceleration equations taken about the 
proper axis through the center of the sphere. 
These equations for moments about the origin 


Fig. 2. Schematic diagram of the coordinate systems 
used. The primed coordinate system is fixed relative to 
the vortex. 


written relative to the system rotating about 
the z axis with angular velocity Q are 


ESS rxevdr=ffnxprda— 


SSI 1X (22 Xov)dr— 
— SS S tX[OX(QXN) ede + 
if Jef 120g (3) 


where r is the radius vector from the origin, 
v is the relative vector velocity inside 
~ the vortex, 
n is the outward normal at the bound- 
ary of the vortex, 
is the pressure, 
is the gravity acceleration vector, 
is the density of liquid in the vortex, 


sn Pots 


and the volume integrals are taken over the 
entire volume of the vortex while the area 
integral is over its boundary. 

The left-hand term is the total rate of change 
of moment of momentum of the vortex about 
the origin and the right-hand terms are the 
moments of the applied forces neglecting 
friction. In the simple case we are considering 
all these terms on the right will be found by 
symmetry about the xz-plane to give net 
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contributions only to moments about the 
y-axis. Consequently we may take scalar 
products with the j unit vector along the y-axis 


and treat only the scalar moments about this 
axis. 

Also it may be noted that all components 
of Coriolis forces, etc., normal to the shell 
will make no contribution since, regardless of 
their magnitude, they have zero moments 
about the origin and consequently about the 
y-axis. 

At least two of the right-hand terms in (3), 
the moments of the centrifugal and gravity 
forces, are functions of the geometrical con- 
figuration only and may be simply calculated. 
For the centrifugal force moment 


Che PO RÉSISTER ES 
Se 9 ee Sera 

= ff fo@xzdr (4) 

which transforming first to r,@, À and then 


to r', O', A’ and taking the volume element 
for unit depth, becomes 


9) 22 
022 
iy, CF —=0 (sin ©’ cos A’ cos « — 
7 fe) o 
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cos @’ sin a) (cos ©’ cos « — sin ©’ cos A’ 


sin a) a? sin ©’ di’ dO’ (s) 


where a is the radius of the spherical shell. 
This gives 


0? a4 x sin 2x 


oe a (cos 9, sin? Oy’) (6) 


My, Ch=6 


The gravity force moment may be simply 
calculated by finding the center of mass of 
the vortex and calculating the moment of 
the force considered as acting at that point. 
The total mass of the spherical cap of unit 


depth is 
M = 0 2a? (1 — cos O,') (7) 


The center of mass is located on the radius 
vector from the origin to the center of the 
vortex at a distance (a/2) (1 + cos @o’) from the 
origin. Thence x 


my, G=—gona sin (1 — cos? @,') (8) 
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The resultant ee ‘force moment is 
m, C=— [ff 1X (e2Xov)dı (9) 


For the case analogous to Rossby’s we take 
v=o a sin ©’ À’ where w is the relative angular 
velocity of the vortex (a rigid rotation) and 
A’ is the unit vector in the direction of the A’ 
coordinate. Proceeding as before 
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(11) 


[2 — 2 cos 9,’ — cos O,' sin? 9, ] 


The most difficult term to calculate is the 
pressure term over the boundary 


m,,p = [[j-(mxpr)dA (12) 


Only the lateral boundaries enter since 
both top and bottom give zero contributions 
to the moment. The principal problem is 
that of assigning proper conditions outside the 
vortex. Once these are assigned a straight- 
forward computation from the appropriate 
Bernoulli equation for the relative frame will 
give the pressures and consequently the total 
moment. For the Rossby assumption of an 
environment at rest relative to the rotating 
shell a very simple consideration gives this 
term. In this case if the vortex cap is imagined 
to be replaced by a cap of the environment 
(water) density o’ at relative rest, the pressures 
outside will be unchanged and the entire 
system will be in equilibrium. 

Consequently the pressure force moment 
over the boundary is just sufficient to equilibrate 
the centrifugal and gravity moments for a cap 
of density 9’, these being the only other terms 
present then. The pressure term therefore 
contributes terms equal to the negatives of 
my, CF and my,G with o replaced by 9’. 
When the entire right-hand side of (3) is divid- 
ed by the total moment of inertia of the cap, 
involving the factor o, the gravity and centri- 
fugal terms combined with the pressure terms 
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give the same terms as above except with o 
replaced by the factor 
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(13) 


If I, is the total moment of inertia of the 
vortex about the y-axis 


at : 
Lee 23 [4 — 4 cos Oy’ — cos 9, sin? ©, ] 
and (3) becomes 
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ment of inertia, 


The angular acceleration of the vortex as a 
whole about the y-axis 
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(14) 


This is entirely analogous to Rossby’s result 
in equation (1), particularly if multiplied 
through by a to give a linear measure pro- 
portional to x. This equation is exact for 
evaluating equilibrium conditions (« fixed) but 
da dx 
dt?’ dt 
as the inertial effect of the environment liquid, 
the Coriolis force for the motion of the cap 
as a whole, or the Coriolian pressure field 
associated with the escape motions of the 
environment. Other factors will be discussed 
in section C in connection with evaluation 
of the experimental results. 

In particular, at the instant of starting from 
rest under the influence of torques, two 
independent analyses to be published else- 
where by Mr Long and Mr Phillips of the 
University of Chicago show that the effective 
inertia of the vortex in the spherical and in a 
plane case is twice that implied in equations 


72:0) 


for == 0 does not allow for such things 
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(1), (2), or (14). This result arises from an 
additional contribution to the pressure force 
on the vortex due to the inertial reaction of the 
environment to the attempted acceleration of 
the vortex. This pressure torque in the above 
case equals one-half the resultant of all the 
others acting on the vortex! so that effectively 
the vortex appears to possess twice the mass 
relative to the original torques. The result 
might have been anticipated from the com- 
parison with that for the plane motion around 
a cylinder given by Lamb (“Hydrodynamics” 
6th Edition, p. 77). Part of a numerical differ- 
ence between RossBy’s (1948) and Davies’ 
(1948) result may be due to this factor since 
Davies carried out directly a series expansion 
for the pressure on the boundary. Since 
we wish to discuss mainly the equilibrium 
criterion for no motion of the vortex, this 
particular factor will not affect the conclusions. 

It will be noted that in this form the Coriolis 
term reduces to Rossby’s expression (equa- 
tion [1]) for Oo’ of the first order of small 
quantities. Thus, substituting the series expres- 
sions for sine and cosine and retaining all terms 
leading to 4th and 2nd order quantities, 
respectively in the numerator and denominator, 
one obtains 


EIN Pe i 9 (0,')2 
3 (HH... fee ae 
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Rewrite the equation (14) as 


da 
dt? 


& 


> sin | 20 Ds —37y< 
a 


> Sy =F 


(15) 


+ 32? s, cos x 


where the size factors s,, Sg, s; depend only 
on Oh. 

The transition from a stable to an unstable 
motion as a whole occurs at equlibrium (at 


1 The one-half factor applies where © = 0’ and is 
altered to 0’/(@+@’) in the case we are actually consider- 
ing. In the spherical case there 1s a*size restriction on 
the vortices to which this result is applicable. 
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least with d«/dt =o also) where the factor 
in curly brackets equals zero and independently 
of the inertial pressures of the environment. 
This condition for a critical relative angular 
velocity ®. of the vortex may be written as 


Wc chew wale SN Oop 
DO NET gto? (®) er (2) RR) 


Thus the boundary of instability for this case 
is a straight line in the w./y 2, g/a 2? plane 
(see fig. 14). In this form the rough criterion (2) 
becomes 
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Table 1. See equation (16) 


, E Linear shear from 
Solid rotation 
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SE 230 |, —~ 225 
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1522929 28.30 
20° 16.40 15.41 
25° 70.506 9.522 
30° 7.291 | 6.314 8.284 PRLS 
as 5.352 4.384 8.065 6.606 
40° 4.093 3.136 7.751 5.938 
45° 3.230 2.284 7.400 5.233 
50° 2.610 1.678 7.044 4.528 
55° ANON 1.234 6.694 3.840 
60° 1.800 .900 6.362 3.181 
65° 1.526 644 6.055 2.559 
70° 1.306 .446 5.778 1.976 
Ty Uc .292 5.537 1.433 
80° 1.004 270 5.341 1.1187 
85° .856 .074 5.208 4536 
90° .750 [6) SNS T o 

3 So 3 53 à : 
The factors = — and = . are given in table 
51 2 Sy 


1. The slopes derived from (17) are almost 
exactly equal to those in the table but the 
intercepts of course are zero in (17). This 
table also gives the corresponding factors 
2 27 ß 4 

2 2 and = for a case which will be important 
256 HS 

for the later discussion. In the actual operation 
of the experiment the fan that generates the 
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vortex motion is usually smaller in radius 
than the polar cap of liquid. While it is quite 
efficient in generating velocities approaching 
solid rotation out to its outer radius there is a 
more or less rapid decrease in relative velocity 
out to the interface between the two liquids. 
In the water itself the relative zonal velocity 
is generally observed to be quite small so that 
the boundary condition of relative rest for the 
top liquid at the interface is probably not far 
from correct. : 

An observationally plausible assumption 
about the vortex velocity field is to take a 
solid rotation to colatitude ©,’ (where ©,’ 
is approximately the angular radius of the 
fan) followed by a linear decrease to zero 
at the boundary colatitude @o’. sı in equation 
(15) and equation (16) is then replaced by 
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yi = yee Th ae ’ 
4—4cos 9,’ — cos, sin? Oy 


{2 —2 cos 9,’ — cos 9, sin? ©, '] 


ce 
=sin Où 
= 


sin ©,’ — sin 9,’ 
[sin ©,’ cos ©,’ (2 sin ©,’ — 3 sin O,’) 
+ 3)sin Oy’ (O,'— @,') + sin? @,' cos 0, 


+ 4 cos O,'—4 cos Zu 
in the above factors. ©,’ is taken as 28° for 
the experimental arrangement and in table 1. 

Further discussion of the significance of 
these relations and the approximations used 
is given in section 3 in connection with the 
experimental results. 


3. Experimental results 


A general view of the apparatus used, 
slightly modified from its earlier form (FULTZ 
1949), is given in fig. 3. A hemispherical part 
of a s-liter round-bottom boiling flask is 
fastened to a vertical hollow shaft. Inside it a 
half-section of a 3-liter flask is floated giving 
a space 1.6 cm deep in which liquid is placed. 
The liquid has the form of a hemispherical 
shell of mean radius to cm. The vertical shaft 


is connected by V-belt pulleys to a split-phase 
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Fig. 3. General view of the apparatus showing the 

hemispherical bowls in the center, mirror above the 

bowls, camera to the right rear above, motor driving the 

bowls to the right rear, Rotoscope to the right for- 

ward, and the motor driving the fan shaft in the center 
front. 


electric motor which rotates the shaft at a 
variety of constant rates depending on the 
pulley combination. (Occasionally, it was con- 
venient to use a series-wound motor whose 
speed could be altered continously by varying 
the impressed voltages, though in many circum- 
stances this arrangement does not give steady 
enough angular velocities.) 

We used various liquids that were immiscible 
with water and of slightly greater density to 
form an isolated cap at the pole. The rest of 
the hemispherical shell was filled with water as 
shown in fig. 4. In order to produce relative 
motion in the polar cap, fans of various forms 
fixed on a separate shaft concentric with the 
hollow shaft were used (fig. 5). This inner 
shaft passes through a packing gland at the 
top of the outer shaft that prevents liquid from 
leaking out of the shell. It is driven by pulleys 
and oval leather belting from a series-wound 
motor operating off a controlled voltage 
source. By these means we were able to work 
in ranges of up to 125 rpm for the outer shaft 
and glass shells and up to 200 rpm in either 
direction for the inner shaft and fan.! 


1 Since both the absolute and relative motions will 
have to be discussed below the following terminology 
will be used: 

“cyclonic absolute” for an absolute rotation in the 
same direction as the hemispheres, 

“anticyclonic absolute” for an absolute rotation in 
the opposite direction, 

“cyclonic and anticyclonic relative” for the same rela- 


Fig. 4. Close-up view of the hemispherical bowls showing 

the lower cap of darker fluid reaching about to 42° 

latitude. Notice the inner bowl which is held down by 

the cross-rods. The space above the cap to the equator 

is filled with water. A contact switch for the revolution 
counter is on the right. 


Some of the data used has been obtained by 
taking sequence photographs from above along 
a line of sight parallel to the axis of rotation 
and at the same point in successive rotations 
of the hemispheres. However, much the largest 
mass of numerical observations has been taken 
visually with the aid of a ‘“Rotoscope” 
patterned after that designed by D. Thoma 
(FISCHER 1931) and used in the investigation 


Fig. 5. Close-up of fan inside the outer bowl with inner 

bowl and liquid removed. The fan shown is that used 

in most of the earlier investigations and, e. g., for fig. 
12. The fan used for fig. 13 was much larger. 


tions of a relative rotation to the absolute motion of 
the hemispheres. | 

Cyclonic absolute rotation at a rate less than that of 
the hemispheres will thus be a relative anticyclonic 
rotation. 
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Fig. 6. Detail view of the Rotoscope without magni- 

fying eyepiece. The image of the hemisphere can be 

seen both in the mirror and in the reversing prism. 
The variable speed motor is on the right. 


of the relative flow in centrifugal pumps. 
This instrument utilizes the fact that if a 
mirror, observed at almost glancing incidence, 
is rotated about the line of sight the virtual 
image seen in the mirror will rotate uniformly 
at twice the rate. In the instrument, a Dove 
reversing prism to give a wider field of view 
replaces the mirror. Our version (fig. 6) con- 
sists of nothing but the prism, a rotatable 
mounting cylinder, and a small variable speed 
motor which drives the cylinder by string and 
pulley or by gears. 

For best results the line of sight of the 
Rotoscope must be along the axis of rotation. 
This requires a rather careful adjustment but, 
once it is attained, the apparatus when rotating 
at any speed can be reduced to apparent rest by 
rotating the Rotoscope at half the rate in 
the proper direction. Furthermore, any wave 
pattern, for example, which is moving around 
relatively unchanged in shape at a steady rate 
can similarly be reduced to rest. Thus the 
relative rotational motion of such a feature 
can be determined either by reducing the glass 
shells to a conveniently low apparent rota- 
tion and clocking its travel between fixed 
meridians or by reducing the feature to rest, 
measuring the required rate of rotation of the 
Rotoscope, and computing the relative motion, 
knowing the rate of rotation of the hemi- 
spheres. Both of these methods have been 
extensively used wherever one or the other 
was convenient. This device has also been 
invaluable in correcting mistaken impressions 
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gained from direct observations of the shells 
at moderately high rotations and from the || 
rather widely separated photographs taken | 
only once per revolution. 

A number of liquids were tried in the 
attempt to find some with a favorable value 
of y, 1. e., with densities close to that of water. 
It is apparent from either (2) or (16) that, within 
practical limits, small values of y will give the 
better chance of obtaining critical angular 
velocities and unstable motions of the polar 
cap as a whole. Several mineral oils and, in 
particular, castor oil with a y of about 1/30 were 
tried, but were found generally to be too viscous 
to give any well defined motions. As a result 
of a suggestion by Dr Gardner of Sinclair 
Research Laboratories we attempted using a 
mixture of monochlorobenzene and toluene 
which turned out to be almost ideal. Both 
compounds have low solubilities in water and 
are completely miscible with one another. By 
mixing in the proper proportions any densi 
between that of toluene (0.86) and that of 
chlorobenzene (1.1) can be obtained. The 
mixtures have viscosities which are even lower 
than that of water and by using a water- 
insoluble dye (Sudan III) the cap can be made 
to stand out clearly from the rest of the 
liquid. If the above liquids are used alone, a 
relatively rapid contamination of the interface 
occurs which causes sticking of the lower cap 
to the glass. The liquids used had to be changed 
for each experiment and the glass very care- 
fully cleaned until we discovered that a pinch 
of laboratory detergent in the water would 
completely prevent the sticking. The water 
must be arranged to wet the whole glass 
surface at least once and from then on a surface 
film prevents any sticking of the organic 
liquids to the glass. 

The experimental procedure varied some- 
what depending upon the type of data being 
taken but in all cases was quite straightforward. 
For instability point determinations, the bowls 
were adjusted to the speed desired and the fan 
then slowly changed from zero relative motion 
to cyclonic or anticyclonic relative velocities 
until the phenomenon wanted was reached or 
until the greatest practicable value was attained. 
Simultaneous angular velocity determinations 
for the bowls and the fan were made at any 
significant points. If the two liquids were then 
unduly beaten together the fan was returned 


Fig. 7. Three-crest pattern observed when fan is rotating 

faster than bowl at a moderate rate. The crests progress 

relatively cyclonically at a rate equal to about .16 of 

the bowl rate of absolute rotation. (Apparent rotation 

of bowl counterclockwise, bowl absolute rotation 

56.4 rpm, fan rate +30.1 rpm relative cyclonic, y = 
0.096 with chlorobenzene as lower liquid.) 


Fig. 8. Two-crest pattern observed when the fan is 

rotating faster relatively than in fig. 7. The crests progress 

relatively cyclonically at a rate equal to about .20 of 

the bowl rate of absolute rotation. (Apparent rotation 

of bowl counterclockwise, bow] absolute rotation 56.4 

rpm, fan rate +41.4 rpm relative cyclonic, y= 0.096 
with chlorobenzene as lower liquid.) 


Fig. 11. Photograph taken one bowl revolution later 

than fig. 10. The off-center crest has moved clockwise 

to the top of the figure corresponding to a relative 

anticyclonic rotation equal to about 0.4 of the absolute 
bowl rate (conditions as for fig. 10). 
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to zero relative velocity until settling out had 
occurred and the experiment was then con- 
tinued. If necessary, the liquids were changed. 

Qualitatively, the results concerning in- 
stability of the polar vortex as a whole were 
clearcut. For cyclonic relative motion of the 
fan, no indications of asymmetry about the 
pole or throwing of the center of gravity of 
the cyclonic vortex away from the pole were 
observed at any velocity conditions which 
could be reached with the apparatus. On the 
other hand, for anticyclonic vortices, at every 
velocity of the bowl except the lowest ones 
there was some relatively sharp value of 
anticyclonic fan speed at which the center of 
gravity of the polar cap was thrown visibly 
away from the pole (see figs. 10 and 11). 
Quantitatively there are difficulties, discussed 
below, which are still not satisfactorily 
accounted for. 

When the fan type pictured in fig. 5 is 
used to drive the vortex the phenomena which 
occur are roughly as follows (the general 
features are insensitive to the number of fan 
blades) : 

A. For relative cyclonic fan motion (e—o’ 
~ 0.1— 0.01 gm cm?) as the speed of the 
fan is increased: 

First, no distinct motion of the interface is 
seen in the Rotoscope though flashes of up and 
down motion can be seen from the side. Then 
a distinct and very stable set of 3 crests is seen 
(fig. 7) which moves around cyclonically rela- 
tive to the bowls but more slowly than the 
fan. The rate of this motion increases with 
the fan speed as appears in fig. 15. At a quite 
sharp fan speed the number of crests suddenly 
becomes 2 (fig. 8). It will be noticed that the 
interface is very smooth and well defined in 
these figures. As the fan speeds are increased the 
number of lobes remains 2, their relative rate 
of travel increases, and ultimately the centri- 
fugal action of the fan tosses the pole free of 
the lower liquid and draws water down in its 
place (fig. 9). Throughout, the center of 
gravity of the denser liquid remains sensibly 
on the pole. Bubbling of the two liquids is 
apparent in fig. 9 and ultimately the polar cap 
is beaten into a froth which, however, remains 
symmetrical in form about the pole. 
 B. For relative anticyclonic fan motion, as 
the speed of the fan is increased from zero 
relative: 
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Fig. 9. Two-crest pattern with center partly thrown clear 
of chlorobenzene and bubbling occurring at high cy- 
clonic fan rates. Note that the pattern is still quite symmet- 
ric about the pole. The crests progress relatively cyclo- 
nically at a rate equal to about .34 of the bowl rate of 
absolute rotation. (Apparent rotation of bowl counter- 
clockwise, bowl absolute rotation 56.4 rpm, fan rate 
+99.6 rpm relative cyclonic, y= 0.096 with chloro- 
benzene as lower liquid.) 


First, the same sequence as in the cyclonic 
case of passage through no discernible motion, 
four or three crests, and finally two-crest 
patterns, is seen. These crests move relatively 
in the direction of the fan but at a lesser rate 
(fig. 15). Very soon the interface is less smooth 
and shows much more irregular, short wave- 
length motion than in the cyclonic case. 

Then, in a range of about $ rpm or less of 
the fan speed, the two-crest pattern suddenly 
goes over into a single-crest pattern of more 
or less regularity. This corresponds to a very 
distinct displacement of the vortex center of 
gravity away from the pole (figs. 10 and 11). 
As the fan speed increases further, the violence 
of the motions increases and also the dis- 
placement of the center of gravity from the 
pole. Even after the liquids have begun to be 
beaten together and up to the speed limits 
of the apparatus the same general feature is 
present. 

The single crest, or center of gravity, moves 
anticyclonically relative to the bowls at about 
20 rpm over the range of speeds seen in fig. 
12. This immediately involves dynamical 
factors which were omitted in the analysis of 
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Fig. 10. Cap thrown off-center when fan is run at a 
high anticyclonic rate. The lower left-hand side is much 
farther out from the pole than the other side, the amount 
being even decreased from its actual magnitude by the 
perspective appearance from above. The fan rate at 
which this appears for the present bowl rate is about 


—95 rpm relative anticyclonic. (Apparent rotation of 

bowl counterclockwise, bowl absolute rotation 56.4 

rpm, fan rate —146.2 rpm relative anticyclonic, no- 

minal y = 0.096 with chlorobenzene as lower liquid. 

The effective y is smaller because of the bubbling 
which occurs under these conditions.) 


the previous section such as effects of the 
escape motions of the environment. The com- 
parison with the theoretical curve in fig. 12 
plus other data brings out several points: 

1° The observed fan speed for transition 
away from the pole is almost constant with 
bowl speed instead of strongly falling off. 
Except at the lower end where, as the bowl 
speed is reduced below about 45 rpm, the 
critical fan speed takes a very sharp rise and 
(with y = 0.10) for smaller bowl speeds the 
cap is badly torn up before a well-marked 
critical throw-off point is reached. 

2° The throw-off point of this type seems 
to be very much less sensitive to the density 
difference than it was expected to be. In fact, 
there is almost no real change in the character 
and location of the curve as the y-factor is 
varied from o.1 to less than 0.01. 

3° The critical fan speed curve shows even 
the wrong slope beyond about 80 rpm bowl 
speed. This may be connected with a circula- 
tion of the escape motion about the cap and 
is also probably related to the fact that the 
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interface begins to climb rapidly up the surface 
of the outer bowl at about this bowl speed. 
This flattens and lengthens the cap edge and 
makes it depart radically from the two- 
dimensional assupmtion in the calculations. 
In addition, the whole outer rim is less efficient- 
ly driven by the fan with a consequent 
decrease of the total unstabilizing Coriolis force. 
The very strong effect this might have is 
shown in fig. 14 by how much the stable 
area is increased under the theoretical critical 
curve when a linear shear correction is applied. 
Qualitatively then, this phenomenon is of 
the expected order but quantitatively the dis- 
crepancies are disturbing, though of course 
understandable, particularly in view of the 
motions of the center of mass. We looked 
around then for possibly some other mani- 
festation of the same effect which might give 
better quantitative comparisons. A very slight 
tendency had been noticed for the cap to be 
very slightly higher on one side than the other 
at moderate and high rotations. This had been 


thought to be due to lack of sphericity of the 
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Fig. 12. Curve showing values of relative anticyclonic 
fan speed where the throwing-off phenomenon of figs. 
10 and 11 first occurs. Off-center conditions occur 
everywhere above the curve except for very turbulent 
conditions to the left of the end of the indicated curve. 
It will be noted that figs. 10 and 11 and the curve to the 
left of 60 rpm bowl speed in this case correspond to 
absolute anticyclonic fan rates while to the right the fan 
rate is relative anticyclonic but absolute cyclonic. Note 
also that the cap radii are different for this figure and for 
figs. 10 and 11. Theoretical curves from equation (16) 
for solid rotation and linear shear assumptions for a cap 
reaching 50° latitude are shown. Since the linear shear 
correction looks too strong, but corresponds best to the 
velocity pattern, it appears definitely that one of the 
escape motion torque factors is involved in this pheno- 
menon. 
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outer bowl which was very carefully centered 
in the attempt to combat it. 

Meanwhile, it was decided to try fan blades 
which reached almost to the interface and had 
leather attachments reaching across almost the 
. full depth of the shell in order to drive the 
cap as efficiently as possible. With this fan the 
violent motions associated with the earlier 
anticyclonic throw-off point had much less 
chance to develop so that this phenomenon 
and also the various waves were only poorly 
developed. However, the tilting seemed to be 
quite regular. Scales were therefore placed 
around the outer bowl near the interface 
edge which enabled its position to be measured 
to within less than a !/; mm. The differences 
of level between the highest and lowest points 
of the interface were then observed as functions 
of bowl and fan speed, the cap being found to 
remain very closely circular. At any given 
condition, the high and low points were found 
to remain at the same longitude indicating that 
the center of gravity was stationary with 
respect to the bowl. This of course better fıts 
the elementary analysis of the previous section. 

The results of one such determination with 
y = 0.10 is shown in fig. 14 where the figures 
entered are the approximate maximum diffe- 
rence of level at the edge of the cap in milli- 
meters. The curves drawn are for 1, 5, and 
10 mm difference, respectively. The general 
form of these curves is very encouraging, 
particularly as to slope, except for two difhi- 
culties: 
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Fig. 13. Curves of equal tilt (in mm difference in level) 

for throwing-off occurring with large fan. The theo- 

retical curve is calculated from equation (16) for assumed 

solid rotation in the vortex. The general appearance of 

these curves is encouraging but a definite reason for 

the throwing-off to begin before it should is still to be 
satisfactorily established. 
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Fig. 14. Theoretical throwing-off curves from equation 
We 2 


yO a @? 
centered near the pole (cos a ~ 1). Throwing-off should 
occur in the upper left-hand region relative to each line. 
The solid lines show the critical curves for caps in solid 
rotation reaching from the pole to the indicated latitudes. 
The dashed lines show the critical curves for caps with 
solid rotation to 62° latitude followed by a linear shear 
to the indicated latitudes. The dotted curve is the cri- 
tical one from equation (17) (ie omitting centrifugal 
terms) for a cap extending to 40° latitude. 


(16) in the 


plane for vortices of various sizes 


a) there is little reason to expect the diffe- 
rence of level to be a regular function of fan 
speed since equation (16) indicates very little 
dependence of the equilibrium near the pole 
on «, the colatitude of the vortex center; 

b) there is no reason in the analysis so far 
to expect any tilting at all until the theoretical 
curve is passed, whereas there are definite tilts 
very much below it.! 

Our present conjecture is that this be- 
haviour is caused by a net surface tension 
torque arising at the interface due to the fact 
that the two bowls are not completely accu- 
rately aligned. It is definitely known that this 
factor exists and is appreciable because at low 
density differences (y = 0.001) the cap can be 
made to climb all the way up one side of the 
shell by shifting the inner bowl slightly toward 
the other side so that the mean curvature of 
the interface is higher on the side where there 
is least distance between the bowls. Work is 


1]t will be noticed that the slope turns up beyond 
about 90 bowl rpm. This appears definitely to be asso- 
ciated with the flattening and rising of the interface 
edge which begins strongly at about 80 rpm and was 
mentioned above. 
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a /°2 attempt to go as far as possible with the 
pr ene cae general case of a moving vortex (disc) and 
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Fig. 15. Relative crest (wave) speeds versus relative fan 
speeds in non-dimensional measure for a bowl velocity 
of 56.4 rpm and nominal y of 0.096. The cap extends 
to about 42° latitude and is being driven by the small 
fan of fig. 5. The figures by the observed points indicated 
the number of crests. Note apparent jump in speed as 
wave number goes from 3 to 2 on the anticyclonic side. 
The transition from 2 to 1 crest is identical with the 
throwing-off phenomenon of figs 10 and 11. The last 
few points toward the left are uncertain. 


underway to track down this effect by varying 
the interfacial tension between the two liquid 
phases. This will be done by varying the con- 
centration of a pure detergent in the water.. 
The interfacial tension will be directly meas- 
ured and it is expected that a close relation 
between its value and the run of the curves of 
equal tilt will be found.! 

We also intend to rotate solid obstacles 
which approximate to the cap form in the 
shell of liquid. By making streak photographs 
of the environment motion it should be possible 
to form a definite estimate of the type of en- 
vironment motions going on (especially for 
the phenomenon of figs. 10 and 11) to check 
and supplement further theoretical analysis 


1 Some of this work has been completed (July 
1950) and has indicated the following: It was dis- 
covered that the inner bowl was actually riding too 
high and was leaving a space of about 2 cm between 
the bowls at the pole instead of the 1.6 cm observed 
at the equator. As a result, the surface tension effect 
seems to be predominantly stabilizing. This may pos- 
sibly be connected with result 2° above, since the sur- 
face tension is relatively independent of the solution 
used and of the density difference. Almost the only 
apparent destabilizing effect then left to explain the 
small early tilts is lack of sphericity of the bowls. 
These results will be collected and discussed in Mr 
Long’ succeeding contribution. Their importance ob- 
viously lies primarily in the problem of estimating how 
far they adversely affect the geophysical implications 
to which the experiments may give rise. 


4. Conclusion 


Qualitatively speaking, the results of the 
previous section can, I think, hardly fail to be 
related to the effect discussed by Rossby and 
to be in some measure a confirmation of his 
prediction that anticyclones can exhibit such 
an instability. Quantitatively, there will remain 
little more to ask for in the simple case if 
the discrepancies of fig. 13 can be definitely 
related to a surface tension or other effect. 

No attempt will be made here to carry the 
discussion of definite meteorological applica- 
tions further, as this should wait on the more 
extensive analyses of environment effects to 
be undertaken in later work. It should be 
mentioned, however, that the occurrences in 
these experiments have relations to so-called 
“dynamic or inertia instability”, shearing in- 
stability, and wave motions in zonal currents 
which should probably be investigated further. 
The crest systems are very regular and show very 
well-marked steady features which appear at 
certain definite speeds. For example, the curve 
of wave speeds in fig. 15, and other results not 
included here, show very intriguing breaks at 
several points where the wave number changes. 
It would be very feasible to investigate photo- 
graphically the internal motions associated with 
these waves and their possible association 
with these changes. Concerning. instability 
phenomena a relatively unexpected result 
is that the interface is roughest with the anti- 
cyclonic vortex for which the shear across the 
interface is cyclonic, while it is extremely 
smooth for the cyclonic vortex where the 
shear is anticyclonic and might be expected 
possibly to reach inertia instability values. 

One additional rather wild speculation may 
be of some interest. Call the critical angular 
velocity, at which the throwing-off of figs. ro 
and 11 occurs, ws. Then if we accept at face 
value in a dimensional analysis the experimen- 
tal independence of y and of the bowl speed 
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Q that is exhibited by ws, © would be pro- 
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The principal parameter re- 


maining would be the gravity parameter A 
so that: 
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For the same sized cap on the earth, taking 
the experimental anticyclonic w,=65 rpm 
and a = 10 cm, the critical angular velocity for 
an earth analogue to this throwing-off pheno- 
menon would be 8x 10-4 sec-1. Unfortu- 
nately, or perhaps fortunately, this is far on 
the absolute anticyclonic side and would not 
be realizable on a large scale by ordinary 
atmospheric processes. 

I do not seriously suggest that this computa- 
tion has any real meaning in the present case 
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but offer it as an-example of one of the sorts 
of things that may eventually be possible to 
do in the way of real meteorological dis- 
covery arising from experimental work. 


Acknowledgments 


Many persons have been involved in the 
work reported in this paper of whom the 
first, of course, is Prof. Rossby to whom I 
am very grateful for support and encourage- 
ment. I am principally indebted also to Mr 
R. R. Long for theoretical discussions and the 
carrying on of much of the experimentation 
and to Messrs W. Bohan and G. V. Owens 
who have done part of the experimental 
measurements and: construction. Much of the 
computing was done by Mr W. D. Kleis and 
the drafting was carried out under the direction 
of Mrs Della B. Friedlander. 


REE ER EINIGES 


Davies, T. V., 1948: Rotatory flow on the surface of 
the earth — Part I. Cyclostrophic motion. Phil. 
Mag. (7) 39, pp. 482—491. 

FISCHER, K., 1931: Untersuchung der Strömung in einer 
Zentrifugal-Pumpe. Mitteilungen des Hydraulischen 
Inst. d. Tech. Hochsch. München 4, pp. 7—27. 

Furtz, D., 1949: A preliminary report on experiments 
with thermally produced lateral mixing in a rotating 


hemispherical shell 
1733. ER 

RossBy. C.-G., 1948: On displacements and intensity 
changes of atmospheric vortices. J. Marine Res. 7, 
pps 175 1.87. 

Rosssy, C.-G., 1949: On a mechanism for the release of 
potential energy in the atmosphere. J. Meteer. 6, 


pp. 163—180. 


of liquid. J. Meteor. 6, pp. 


The Lowest of the Ionized Layers of the Upper Atmosphere 
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Abstract 


The ““Stockwerk-Gliederung”’ of the Earth’s atmosphere introduced by PENNDORF and FLOHN 
is a promising method to describe the different atmospheric regions from the nearest to the 
ground up to the highest levels. The region above the upper inversion (at ca. 75 km) up to 100 
km height is generally called the D-layer which means the lowest of the ionized regions in the 


upper atmosphere of the Earth. 


The behaviour of the D-layer, its most reasonable theory of formation, and its function with 
regard to the propagation of electromagnetic waves are discussed. Features related to solar- 
terrestrial. phenomena are shown. Theories and ideas by various authors are compared with 


each other. 


The complete investigation from which this paper gives an extract is being published else- 


where. 


1. The three principal known ionized-layers 
of the Earth’s atmosphere (D-, E-, and F-layer) 
show different behaviour owing to their 
various altitudes and origins. It is reasonable, 
though not generally accepted, to ascribe the 
uppermost layer (F) to the shortest wavelengths 
of an ultraviolet spectral band between soo and 
1300 A, while with increasing penetration 
into the lower regions of the atmosphere longer 
and longer wavelengths of that band of solar 
spectrum are responsible for the formation of 
the lower layers. Photochemistry has revealed 
to us a number of dissociation and ionization 
processes acting on atmospheric oxygen and 
nitrogen (table 1). Thus the layers formation 
is clear in principle, although many difficulties 
as to quantities still have to be overcome. 
Often there must be several processes operating 
simultaneously. 

All reactions with nitrogen and atomic 
oxygen occur at wavelengths below 1 000 À 
as can be seen from table 1. The dissociation 
and first ionization potentials, however, of 


molecular oxygen lie above 1 000 A. Accord- 
ing to PENNDORF (1949 b) there is practically 
no O2 above the height of 115 km, ice. all the 
latter reactions take place in the E-layer and 
below it, requiring a radiation which will 
not be appreciably absorbed down to that 
level. Between 115 and 90 km, rather disruptive 
at 108 km, oxygen is dissociated. 

2. In this paper the photochemical processes 
of the D-layer are of interest. Reactions 7) and 
8) are at our disposal. The absorption of 
oxygen is very weak between 910 and 1 300 
A, likewise the weakening by nitrogen is 
negligible. This allows a considerable portion 
of this spectral band to penetrate down to 
100 km and still lower in order to ionize the 
oxygen molecule and further to dissociate it 
into O* via excited states of the oxygen atom. 
The excitation-energy is used for the ionization 
of an oxygen atom reacting by three-body 
collision. Assuming a probable model of the 
atmosphere (temperature-height function) and 
a reasonable effective recombination coefficient 


"2 


nn. 


THE LOWEST OF THE IONIZED LAYERS OF THE UPPER ATMOSPHERE ISI 
Table 1 
1) N + hy>N,t+e A< 661 À a=? F,-layer 
2) © + kw>OÖt +e < 667 Ä =? F,-layer 
3) O + hy>Ot +e 722 IN = F;-layer 
4) O, + m>0,+ + Gag IN = 500 E-layer 
5) N +moN,t+e EN — 10—50 E-layer 
6) O + hr OT +e = oro À = 260 F,-layer 
7) O, + hyvn>O,+ + e < 1010 A = 0.3 D-layer 
8) O, + hy>O* O < 1330 À ON O atom layer—D-layer 
9) O3 + hy>O* + O < 171$ À = 200 O atom layer 
10) ©, + hy+O O 2124 À = 103 O,-layer 
Table 2. Electron density in the D-layer obtained by processes 7) and 8). — Model atmospheres 
I) (©) (D). ne [cm?]. i 
(A) (B) (©) (D) 
h [km] 
Pa = Ch 40° 85° 80° 40° 80° 
| 
130 4.5, 102 5.0 - To? 1.4 * 102 Sele ro” 6.6 + 10? Ie aoe L ; 
125 6.3 * 10? 7.2 * 102 1.9 + 103 8.9 : 102 I.I 103 1.8 + 103 | Contri- 
120 1.0 108 Te Serko: 2.8 ° TO? 1.0 : 10° 210 103 3082108 au 
IIS 1.9 + 108 2.3 * 10 4.6 + 10% 2 20H00 3.8 - 103 7.9, 10° to the 
IIo AIO ASS 103 SAC 6.3 “T0 ED TO 1.6 - 10% Saye: 
105 7.8 - 108 8.7 - 103 1.4. TO 1-2 010% I.3 : 10% 1.8 : 104 
100 TA rot 1.6 : 104 1.3 104 1.5 : 104 1.7: 104 I.I * 104 
95 1.8 - 10% 1.9 10% ESKETOS 1.5. 10% 1.4 : 104 3.4: 10% 
90 I.9 > 104 2.0 : Io? 2.3 10°. 1.1.1048 Tater Os §.2 * 102 
85 FOTOS T:4 "10% To 4.0 * 108 5-087 102 1.8. ro! 
80 5.6 > 108 3.6 - 108 — PORTO 6.0 > 102 — 
75 4310 1.8 : 102 — — AS O0 _ 
70 Duy OÙ (x 10°) — — — = 
Max. Dens. 90 90 103 97 IOI 107 km 


Model atmosphere (A) means an average distribution of temperature with height taken from a number of 
various papers, while (B) represents a mean atmosphere in temperate latitudes in winter, and (C) and (D) a mean 


atmosphere in temperate latitudes in summer. — It is referred to original paper. 


with these two processes of ionization and 
dissociation the vertical distribution of electrons 
can be calculated if the solar energy and the 
absorption coefficients in the spectral band un- 
der consideration are given. The result is 
shown in table 2. 

Diurnal and seasonal variations are affected 
by changes of the Sun’s zenith distance; and, 
moreover, seasonal variations are influenced 
principally by the temperature changes in the 
upper atmosphere. Therefore calculations of 
D-layer at the same solar height, but with the 
summer temperatures in the one case (D) and 
with winter temperatures in the other (B) 
produce quite different results, while the 
difference between maximum solar height 
(noon) in winter (B) and summer (C) becomes 
rather small. This difference is evident only 


x = solar zenith distance. 


at the lower border of the layer. The upward 
movement of the isobaric surfaces in summer 
causes increased absorption, already at higher 
altitudes, of the energy flux of solar radiation 
although this flux is greater at smaller zenith 
distances. It is quite striking that large annual 
variations of the maximum electron density 
in the D-layer are suppressed; the lower 
border, however, makes a marked yearly 
cycle. The “penetration frequency” of the 
D-layer at winter noon is not much lower at 
summer noon, when it may occasionally 
approach 2 Mc/s (in middle latitude). This 
was determined — though approximately 
of course — from measurements of the attenua- 
tion limit of pulsed high frequency signals 
with continually changing wavelength (multi- 
frequency record). The diurnal variation, 
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however, of this attenuation limit is quite 
impressive. Normally there is no damping 
at all at D-layer heights during the night. 
The electron density decreases to very low 
values. 

How far can the presumptions involved in 
the calculations be supported by other observa- 
tions or theories, and what are the uncertain- 
ties? The UV-spectral band of solar radiation 
cannot be measured direcly by a radiation 
‘receiver on the ground. The energy distribu- 
tion in the band below 2 900 A was obtained 
indirectly by investigating various photo- 
chemical processes in the stratosphere and 
ionosphere. The ozone formed in the upper 

art of the ozone layer can be measured be- 

cause of the absorption of solar radiation, and 
the energy necessary for its formation can be 
estimated. Similar reasoning can be applied 
to the E-layer and F-layer, the maximum 
electron densities of which are measured by 
the echo-method. The study of formation 
processes in the ionospheric layers clearly 
shows that the UV-continuum of a 6 000° K 
radiator does not suffice for all the processes 
involved. In other words the radiation of the 
Sun’s photosphere is not capable of building 
up and maintaining the Earth’s ionosphere, 
or even the ozonosphere, as they are known 
from observation. It is indispensable when 
considering the different processes or their 
spectral bands, to assume different energy 
sources (ALLEN, 1948) and to place them 
at different heigts in the solar atmosphere 
(chromosphere and corona) where especially 
hydrogen, helium and calcium and other 
elements, emit in lines and bands. The 
energy flux densities applied in this investiga- 
tion are in good agreement with those of 
other papers (DÜTSCH, 1948; KIEPENHEUER, 
1945, 1948; KUIPER, 1949; NICOLET, 1945; 
SCHRÖER, 1947). For processes 7) and 8) the 
results require quantum flux densities of a 
continuum radiator (black body) with 5 740° 
resp. 5 200° K. | 

The absorption path in the atmosphere 
depends not oly on the Sun’s zenith distance, 
but also on the gas density along the path. 
Besides the composition of the atmosphere 
it is essentially the temperature distribution 
with height which determines the atmospheric 
model to be applied. Up to so km it seems 
to be quite well known, above that height 
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uncertainty increases. The temperatures used 
in calculations are in good agreement with 
values obtained from pressure measurements 
(PENNDORF, 1949a) made during German 
V2-rocket ascents performed in New Mexico 
(USA). A result important for the D-layer 
was the low temperature of about 180° K 
found in the upper inversion at 75 km height. 
Such a low value has not yet been published 
by any other author. 

The absorption coefficients wo for oxygen 
and nitrogen in the region between 910 and 
1 300 À demand more experimental investiga- 
tion. At the present time no more can be 
stated but that they must be very small in 
order to pass sufficiently intense radiation of 
that spectral range down to 100 km and less. 
Other theories, as for instance the ionization 
of sodium or the sudden reinforcement of the 
D-layer by X-rays, will not be discussed here 
(see end). They require quite different spectral 
suppositions and, moreover, can not be 
accepted without objections. Measurements of 
the absorption coefficient of oxygen have 
advanced to the spectral region in question 
(LADENBURG u. VOORHIS, 1933) so that appreci- 
able changes of the numerical magnitudes 
and, in turn, serious corrections of the structure 
of the D-layer are scarcely expected. 

Finally, one must discuss the influence of 
the effective recombination coefficients at 70 
to 1io km which were used for the calculations. 
Unfortunately, one must depend upon inter- 
polation. Measurements during various solar 
eclipses have yielded a rather safe value for 
the E-layer, ie. an altitude of 115 km. Also 
three-body collision factors for the ozone 
layer are known from which an idea of the 
reaction speed may be determined. Thus the 
effective recombination coefficient in the D- 
layer is limited to a certain range which is 
rather definitive. Of all the assumptions made, 
that dealing with the recombination coefficient 
is of the greatest uncertainty. 

Fortunately, however, there exist several 
observational results confirming the calcula- 
tion. Records of very long waves (kilometer 
waves) serve as a very good indicator of the 
behaviour of the D-layer, since these waves 
are reflected from its lower border. Thus the 
height of this lower border is known, and it is 
possible to observe the increase in height of 
this lower border at sunset on undisturbed 
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days and after shortlived reinforcements of 
the D-layer (sudden ionospheric disturbances). 
From both we get recombination coefficients 
in good agreement with each other. The 
following values have been used for calculating 


the D-layer. 


Table 3. 

h km kp [em sec-1] 
130—IIO 10 $ (measured) 
110—100 TON! 

100— 80 7:10 7 (measured) 

80— 70 SCO 


The infrequent observation of fade-outs on 
medium wavelengths (1 149 kc/s) during a 
chromospheric eruption gives a ‘“‘penetration 
frequency” of the D-layer of about 570 kc/s, 
corresponding to a maximum ionization of 
Ne =4 10% [cm]. This value has the same 
magnitude as the average noon value (5 - 103 
to 104 [cm |). The observation was made in 
the afternoon at a zenith distance of the Sun 
Gere = 65°. 

Investigations of short distance fading of the 
wireless station at Beromünster by GERBER and 
WERTHMULLER (1947) also confirm the pre- 
ceding calculations. At vertical incidence its 
frequency lies not far from the noon “penetra- 
tion frequency’ of the D-layer. After its 
disappearance in the evening the echoes re- 
turn from the lower E-layer. 

In this connection it may be emphasized that 
in the morning the D-layer distinctly emerges 
from the lower region of the E-layer while 
later on it represents a separate layer with 
maximum ionization at ca. 90 km height, 
depending on season and latitude. Dissolution 
during the evening, of course, begins from the 
lowest regions leaving electrons after a time 
in the lowest E-level (roo km), until a night 
there is only a single layer left at that height. 
Recombination in the E-layer is not suffi- 
ciently rapid to make it disappear during one 
night. This is possible only during the polar 
night, and only occasionally has it been 
observed that the electron density decreased 
below 10? [cm?]. 

Good records of a pulsed fixed frequency 
sometimes show a laminated structure of the 
E-layer at 100 to 110 km in the morning hours, 
whereafter the echoes gradually move down- 
ward becoming invisible at levels below 
95 km in height. Here the influence of attenua- 
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tion has become so strong as to make the 
return of D-layer echoes impossible. Indeed, 
one can characterize these very low echoes as 
coming from the D-layer during formation. 
The gradient of charged particles at the lower 
border of the layer is, at the onset, very high 
and therefore favours reflection of the electro- 
magnetic wave. 

3. The undisturbed D-layer follows, as 
can be seen, the daily solar course practically 
without delay. Short periods of slightly 
fluctuating ionization are superposed on it 
and probably are due to corresponding changes 
in the lower chromosphere (granulation). 
These fluctuations of the D-layer reveal 
themselves very distinctly in the records of 
the field intensity of kilometer waves. Erup- 
tions on the Sun (WALDMEIER, 1937) are 
observed to be of all intensities from the very 
frequent, but weak, fluctuations to the rare, 
but very strong, ones whose effect can be 
detected on the Earth in various ways. Most 
striking is the sudden, often complete, fade-out 
of short-wave transmission on certain com- 
munication channels during such a chro- 
mospheric eruption of high energy. During 
such occurences the resonance emission of 
solar hydrogen can rise for a short while up 
to 106-times the intensity of the black-body 
radiator at 6 000° K in the spectral band under 
consideration. Simultaneously the propagation 
of very long waves can often be observed to 
improve. The magnetograms of the geomag- 
netic observatories on the daylight side of the 
Earth frequently show quite typical variations 
(impetus of some gammas with gradual return 
to the normal state after some minutes de- 
pending upon the duration of the disturbance). 
They can be detected very easily when they 
rise above the general activity of the magnetic 
field. Therefore, they can be most easily 
detected on magnetically quiet days. 


Chromospheric eruptions on the Sun can 
be observed in the light of the BALMER series 
(He, B, y) of hydrogen. Simultaneous emission 
in the far ultraviolet must be assumed though 
it can not be observed directly; only its 
photochemical effect upon the upper at- 
mosphere of the Earth (oxygen molecule) may 
be detected. Ly-« lies just in the gap of low 
oxygen absorption so that energy of the 
hydrogen emission line at 1 216 A with a 
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Table 4. Ionization of the D-layer during a strong sudden ionospheric disturbance. 
ne [em 2]. — Model atmosphere (C). — x = 40°. 
ee SER So a I ere un a 


Processes 7) + 8) | 


h (km) Process’7) ie 5) Pro) 104-fold contin. + 
10%-fold continuum Ly-system ee, 
130 8.1 - 104 2.0.2 10% 1-20 | 
125 158 Lo 6.6 : 104 2.0 * 105 Contribution 
120 2.2 10° 7222792 : 3.3 * 105 to the 
115 3.8 - 10° 22 MO 5.6 + 10° | E-layer 
IIo 6.9 * 10° 4.1 TOP TECÉTOP 
105 120100 72 OP 2.0 108 
100 1.0 = 10° T.12 10° 2.8 OS 
95 HS 2708 1.301708, 2.8 - 10® 
90 17.2708 1.3. : 10° DESE TON 
85 GS 100 TT TO 2.6. 10° 
80 4.0 - 104 3.2" 10° 3.6 16° 
WS 6.1: 102 TS ATOS 1.525107 
70 10 ao 8.0 Iol 
Max. Dens. 99 95 94 km 


width of some tenths of Ängströms can Sun probably comes from the higher levels of 
produce a considerable amount of ionized ‘the chromosphere, as can be deduced from the 
oxygen atoms down to 70 km, by way of fact that limb darkening is only very faint. 
dissociation of excited oxygen molecules. In fact, short wave fade-out was observed 
It has, however, been stated that not every when eruptions occurred at the limb of the 
strong increase of the D-layer ionization Sun. In most cases the emission centres are 
is accompanied by emission in the BALMER found near sunspots or active centres which 
series. Since emission in He, 6, y and in the soon became visible as sunspots. It must be 
LYMAN system are supposed to occur simul- supposed, however, that the continuum radia- 
taneously, it can be concluded that not every tion originates from the granulation, that is 
increase of D-layer ionization is due to Ly-« from the lower levels of the solar atmosphere 
(process 8), but that an intensified emission of where it can sometimes increase considerably 
the continuum is also able to produce an in limited regions covering say 1/1000 of 
abnormal D-layer with similar features to those the solar disc. It may be that these areas of 
described above. Table 4 shows the electron higher energy production represent the pre- 
density of the increased D-layer produced by decessors of later chromospheric eruptions, 
process 7) + 8) [with a continuum radiator where according to KIEPENHEUER (1948) the 
supposed to have been intensified by 10%], turbulent movement of ionized gas produces, 
by process 8) with Ly-x alone, and finally by its retardation in the strong local-magnetic 
y the cooperation of both processes. The field, sufficient energy to make possible the 
continuum extends from 910 to 1300 À heavy hydrogen emission occurring in the 
containing ca. 3 * 101% quanta per cm?—sec of upper regions of the chromosphere. Pe? 
extraterrestrial radiation. NICOLET (1945) gave later stage in this enormous energy transfer 
an estimation of maximum possible LYMAN a rising temperature of the corona, accom- 
energy which is involved in this calculation. panied by increasing penetration frequency 
A very strong emission in Ly-« yields an of the F2-layer, can be observed. ALLEN (1948) 
extraterrestrial quantum flux density of 3: 101° reports on the chronological sequence of solar 
[hy cm“? sec]. Thus a weak UV eruption and terrestrial events with regard to the relative 
in the chromosphere produces in Ly-x twice sunspot number R. The author was able to 
the energy of that produced by the granula- show that an intensified D-layer caused by a 
tion continuum between 1 000 and 1 300 A! chromospheric eruption would occur some 


The excessive hydrogen emission of the days before the zero day for R. The maximum 
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of R is found on the average § days after a 
fade-out effect. Unfortunately, the material 
was not sufficient for statistical investigation 
when divided into cases with and without 
emission of Ha, that would, so to say, separate 
LyMaN and continuum processes, but it would 
require the observation of the Sun simulta- 
neously with the fade-out in each single case. 
If a separation of the two processes were 
possible in this way any delay element existing 
between the phenomena could be seen in 
terms of R. The continuum should be expected 
somewhat before the zero day of R. 

The calculated effect of the two processes 7) 
and 8), when excessive UV radiation is de- 
parting from the Sun, differs less for the ioniza- 
tion maximum of the layer than for the height 
of the lower border which moves down- 
ward quite remarkably. Ly-« penetrates deeper 
into the atmosphere causing a greater displace- 
ment of the lower border of the D-layer. The 
experiments already mentioned in connection 
with ascents of German V2-rockets in New 


Mexico (NEWELL, 1948) were a remarkable: 


success, in that field intensity measurements 
from two transmitters on two different wave- 
lengths were obtained just on a day with 
worldwide fade-out (MÖGEL-DELLINGER effect). 
The lower border of the D-layer, where on 
both frequencies very strong, but different 
absorption suddenly began, was found rather 
distinctly at 64 km height. Qualitatively this 
observation agrees with the calculation though 
the latter, representing a general case for 
average conditions, differs from the observed 
special case quantitatively. The place of 
ascent is situated at 37° N, the start was in 
March. The model atmosphere should be 
fitted to this special experiment. 

An interesting question may be raised: 
During a chromospheric eruption are there 
frequencies, below the penetration frequency 
of the intensified D-layer, but suffering so 
little absorption as to be reflected from the 
layer (though perhaps only sporadically) so 
that echoes can be received at the ground? In 
other words the frequency must be sufficiently 
higher than the collision frequency of the 
ionized gas particles with the molecules. In 
order to make this effect possible a rather high 
ionization of the D-layer with regard to its 
lower border is necessary. The calculation 
yields maximum figures ohne 3,10% (cin. | 
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at 95 km altitude. This maximum, however, 
is rather extended covering a height interval 
of 10 to 15 km. At greater Sun heights the 
absorption of frequencies of 10 to 15 Mc/s 
is complete at 85 km (collision frequency 
at summer noon about 2 - 10% [sec~4]), 
and there is no measureable reflection of 
energy. This, however, does not hold in the 
early morning hours (or shortly before sunset) 
when the lower border of the D-layer is 
high, i.e. at a level where the collision fre- 
quency is still sufficiently low (summer 
morning about 5 : 105 [sec-1| at 100 km). A 
very strong Ly-« eruption would for a short 
instant produce a penetration frequency of 
10 Mc/s in the D-layer at that height. Echoes 
should thus be possible if the ionization does 
not reach farther down. The author can refer 
to such a case, though not of such extreme con- 
ditions, which agrees well with the calculated 
data of the intensified D-layer. Occasionally 
during the partial solar eclipse in the morning 
hours of April 28, 1949, there were echoes 
from a height of 87 km, approximately, on 
the 3 and 5.8 Mc/s records, lasting about 5 
minutes (ionization ca. 3 : 10° [cm *], colli- 
sion frequency ca - 2 - 105 [sec-1]). The lower 
border of the layer owing to the great solar 
zenith distance was still very distinct, and 
the echoes must be supposed to have been 
favoured by unusual irradiation conditions in 
the ionosphere that morning. 

The transfer of ultraviolet radiation energy 
in the spectral band (910—1 300 A) under 
consideration, which takes place in the at- 
mosphere between 100 and 75 km height, into 
kinetic energy of the molecules by means of 
the ionization and dissociation products causes 
the temperature of this lowest part of the 
ionosphere to rise but negligibly under normal 
conditions (temperature in the upper inver- 
sion 180° K). The increments of energy input 
of short duration during chromospheric erup- 
tions, can, however, produce a temperature 
rise of 30°C in half an hour (maximum of 
layer at 95—100 km) in winter. In summer the 
influence is considerably weaker. Recent 
papers (Baur, 1949, DULL, 1948) point out 
the importance of these sudden disturbances 
of the baric state on the daylight side of the 
Earth on atmospheric movement. Statistical 
investigations of these phenomena will, no 
doubt, reveal interesting results in future. 
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4. Twenty years ago CHAPMAN (1929) 
mentioned the effect on solar radiation pro- 
duced by absorption by matter emitted from 
disturbed solar regions. RICHARDSSON (1944) 
and later Brück and RUTLLANT (1946) could 
observe a very weak self-absorption of the 
hydrogen radiation (Ha) in the corpuscular 
stream envelopping the Earth when sweeping 
over it. From the Dopprer effect (violet 
shift) they could calculate the velocity of the 
particles. 

In quite another way the author (STRANZ, 
1950) succeeded in demonstrating the absorbing 
effect of corpuscles envelopping the Earth on 
their way from the Sun by utilizing half a 
year’s fade-out observations' together with 
solar and geomagnetic data (RYDBECK and 
STRANZ, 1949). Using the method of sample 
days of BARTELS an equal number of days 
with and without MÖGEL-DELLINGER effect 
(SID) were compared, and it is possible to 
show that the effect preferably occurs on days 
with a mean magnetic character somewhat 
lower than average. In other words the 
probability of fade-out is reduced slightly 
due to the weak absorption of active UV ra- 
diation when the magnetic character is higher 
than normal. The countertest used in checking 
the reality of the phenomenon (the selection 
of an equal number of days free from sudden 
ionospheric disturbances chosen from amongst 
the different months when the SID-days 
occurred) proved to be positive without any 
doubt. Higher magnetic activity, i.e. increased 
corpuscular emission from the Sun towards the 
Earth, caused a bit more intense absorption of 
short wave radiation (hydrogen resonance 
radiation) and less probability of disturbances 
in the D-layer of the kind described in this 
paper (SID). 

In fact, extraordinarily impressive fade-out 
phenomena were observed on days when 
magnetic activity was absent, i.e. when 26—30 
hours after the MÖGEL-DELLINGER effect K — o 
or 1. Not until after so hours did the percentage 
of K = 5 rise again in the group of 130 cases 


ı The observations were carried out by the author 
when working at the Chalmers Institute of Technology, 
Gothenburg. The author wishes to express his gratitude 
to the Stockholm Observatory at Saltsjöbaden and to 
the Geomagnetic Division, Kgl. Sjökarteverket, Stock- 
holm, for their kind permission to use solar and geo- 
magnetic data which they have collected. 
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treated statistically. K = 5 was chosen as a 
reasonable limit between the less disturbed 
days and those with a corpuscular stream 
directed towards the Earth worth while 
considering. 

From an analogous investigation, using 
sunspots distributed over the visible hemisphere 
of the Sun and divided into 12-day sectors, a 
similar result is obtained. The correlation 
seems to be lower, because one must not 
expect each spot in sector + 1 (west of central 
solar meridian) to throw matter in the direc- 
tion toward the Earth. Magnetic effects give 
particular security that the Earth will be en- 
velopped by corpuscles from the Sun which the 
day before had departed from it. The pa- 
rallelism between the results obtained from the 
experiments with magnetic K-indices from 
the Geomagnetic Division of the Sjökarte- 
verket, Stockholm, and from the solar data 
from the Saltsjöbaden Observatory gives a 
particularly good confirmation of the hypo- 
thesis developed concerning the absorption 
effect. 

5. A. and E. Vassy (1948) have ascribed the 
formation of the D-layer to the ionization of 
Na detected in the upper stratosphere or lower 
ionosphere. According to their opinion Ly-« 
will be consumed already at altitudes above 
100 km so that the formation of the D-layer 
in this way will not be possible. On the 
contrary they believe radiation of wavelengths 
greater than 1 800 A is responsible for the 
ionization of Na; this radiation is said to be 
the emission of Cat. In fact, the D-layer 
usually shows a cloudy structure which could 
be identical with Na-clouds. PENNDORF (1947), 
however, found from calculations of the border 
of the ozone shadow that sodium luminescence 
observable at twilight comes from the E-layer 
(r1c—120 km). What occurs during the day 
cannot be stated from these considerations. 

SIEDENTOPF (1948) recently published a new 
theory of formation of the D-layer and postu- 
lates an X-ray radiation of 1—2 À originating 
from the Sun’s corona as the causitive agent. 
The very hot corona can deliver enough 
electron energy to produce X-rays. According 
to this author this radiation only is capable 
of ionizing molecules in the atmosphere at 
the level of the D-layer. All other types 
of radiation are not able to penetrate as far 
down as this altitude. 
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The statistical result, that D-layer disturb- 
ances preceed corona disturbances (ALLEN, 
1948) when referred to sunspot number R 
zero-day, makes this explanation questionable. 

In addition to these theories which deviate 
essentially from the one described here, there 
were earlier attempts of other authors to ex- 
plain the D-layer formation which agree in 
general with this one. Only certain differences 
in the photochemical processes and in the 
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amounts of energy assumed in the calculations 
exist, but these produce no significant differ- 
ence in the results. 

This investigation, described in full else- 
where, should contribute some essential ideas 
and observations toward a better under- 
standing of the D-layer which has been only 
indirectly accessible to experiments. For a 
more complete description one is referred to 
the original paper. 
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The Wind-driven Circulation in Ocean Basins of Various Shapes 


By W. H. MUNK! and G. F. CARRIER? 


(Manuscript received 6 March 1950) 


Abstract 


The ocean circulation induced by zonal winds has been derived for a triangular ocean basin. 
The method, which involves the ‘‘boundary layer’’ technique, can be extended to a more 
general wind system, and to ocean basins of arbitrary shape provided the boundaries are not too 
irregular. In low and mid-latitudes the principal effect of an inclination, ©, of the western bound- 
ary relative to a north-south direction is to widen the western current (Gulf Stream) by the 
factor (sec @)1/3 and to reduce the current intensity by a similar factor. A variation in the 
circulation pattern associated with the variation in the value of the lateral eddy viscosity is noted. 


Introduction 


The wind-driven circulation in a rectangular 
ocean basin has been discussed by HıDakA 
(1949) and by one of us.3 The actual shapes of 
the North Atlantic and North Pacific Oceans 
deviate greatly from that of a rectangle, but 
resemble somewhat more closely triangular 
or semicircular basins, the two cases consid- 
ered in this paper. The solutions obtained indi- 
cate the effect on western currents, such as 
the Gulf Stream and Kuroshio, of the orienta- 
tion and curvature of the coast line. 


The equation of mass transport and boundary 
conditions 


The stream lines of mass transport, y, obey 


the equation (WOC) 


(aa Ep g 


Ox 


) y = — curl,t, (1) 


! Institute of Geophysics and Scripps Institution of 
Oceanography, University of California. Contribution 
from the Scripps Institution of Oceanography, New 
Series, No. 483. This work has been supported by the 
Office of Naval Research. 

2 Department of Applied Mathematics, Brown Uni- 
versity, Providence, Rhode Island. 

3 W. H. Munk (1950), On the Wind-driven Ocean 
Circulation, published in the Journal of Meteorology. 
This article will hereafter be referred to by the symbol 
(WOC). 


where A is the kinematic eddy viscosity per- 
taining to lateral stresses, V4 = d4/dx4 + 
+ 204/dx?dy? + 9*/oy* is the biharmonic ope- 
rator, 


nth ine & 
p= =: cos p ©) 


is the rate of change northward (positive y) of 
twice the vertical component of the earth’s 
angular velocity 2, f the Coriolis parameter, R 
the earth’s radius, g latitude, and curl;r the 
vertical component of the wind stress curl. 
The positive x-axis extends eastward. 

Equation (1) is essentially an integrated 
version of the vorticity equation, and expresses 
a balance between three torques: the lateral 
stress torque, the planetary vorticity, and the 
wind curl. The equation has been derived 
(WOC) by integrating the equations of mo- 
tion (including horizontal and vertical tur- 
bulent stresses) from the surface to some 
depth where the motion and horizontal pres- 
sure gradients essentially vanish, and then per- 
forming the curl operation. Observations have 
shown that the velocity of the wind-driven 
currents at depths of one or two thousand 
meters equals but a small fraction of the sur- 
face velocity. 
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Fig. 1. Dimensions of triangular ocean basin. The wind 

stress 7 varies from — J’ at y = o (p= 15° N), to + I" 

at y =a (p= $0° N). A unit distance equals 35° 

latitude/z, or 1032 km. In these units r’ is the latitudinal 

width of the ocean basin, and y§ the distance from the 
western boundary. 


Two considerable simplifications result from 
working with the vertically integrated equa- 
tions. In the first place, we can examine the 
case of a baroclinic ocean without having to 
specify the vertical distribution of density 
and current; secondly, only the vertical stresses 
at the upper boundaries need to be specified, 
and those can be computed from the known 
wind field over the ocean. 


For a wind system we set 
El COS Ny, t —0 (3) 


giving maximum easterlies at y = o (latitude 
15°) and maximum westerlies at y = x/n (see 
fig. 1). The method can easily be extended to 
an arbitrary distribution of zonal winds, and 
under certain conditions to a general wind 
stress field (WOC). 

Assume a triangular occan basin with the 
dimensions shown in fig. 1. For the boundary 
conditions along AC and BC we choose 


(42,8) 
where » is normal to the boundary. Along AB 


(22 y/0v?)bary = 0. (Sa, b) 


(0 v/ov)bary == 6 


YVbdry — O, 


Wodry == (0) 


Equations (4a) and (sa) state that the bound- 
ary ABCA is a stream line. According to 
(4b) no slippage takes place along the land 
boundary ACB, whereas according to (sb) no 
lateral shear exists along the sea boundary 
AB. The boundary conditions are equivalent 
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to those imposed on a rigid plate clamped 
along the coast line ACB and supported along 
the latitude line AB (WOC). 


Solution by boundary layer method 


Introducing the non-dimensional parame- 
RSR ern == ny Van, 
n4V'4 = V4, and writing k? = ß/A, equation 
(1) becomes 


(ve) psy (6) 
OX 

where y = n/k is the ratio of the wind wave- 
number # to the “Coriolis friction” wave- 
number k. The distance x/n between the easter- 
lies and maximum westerlies is roughly 3500 


km; it will be shown that 


27 (sec O)4/3 
V3 k 
is roughly the latitudinal width of the western 
currents, say 350 km. Setting O= 60° for the 
Pacific gives y = .035. 

Inspection of equation (6) now shows that 
the problem belongs in the “‘boundary layer” 
category. The leading term (which contains 
higher order derivatives than all other terms) 
has a coefficient y3< 1. The asymptotic 
solution of equation (6) (i.e., the solution which 
becomes more and more exact as y — 0) 
can therefore be anticipated to give a good 
approximation to the exact flow except in 
the vicinity of the boundary intersections. 

The conventional procedure for obtaining 
this solution is the following. We find a parti- 
cular solution of equation (6) by neglecting 
the term containing y?. Such a solution is 


y'© == — x" sin y‘ + d sin y’. (7) 


In order to obtain the homogeneous boundary 
layer solutions which lead to the determina- 
tion of d and the satisfaction of the boundary 
conditions, it is convenient to introduce the 
following coordinate system (fig. 1): 


nn tan. ©), 
Y 


ue 


C = Ve (8) 
so that 


py’ = (—y &€—C tan 9 +d)sin£ (9) 


ae 
(10) 


= = oe 
y a) 


d 
= Di (p 


tan 
D : Po) p 


(11) 


is a correction involving the second derivative 
of the Coriolis parameter. Here @ is the lati- 
tude at which y = € = 0, about 15° N. The 
correction is small at low latitudes and amounts 
to 25 per cent at so? latitude. 

It is seen that variations in the y-direction 
involve € explicitly as well as implicitly 
through £. Since y < 1, it follows that only 
the latter need be considered, so that 

oa 


Mn nA p—8 __ 


(12) 
where 


pri DENE) (13) 
With these substitutions the asymptotic 
form of the homogeneous equation becomes 


ot 0 ; 
(ri) 0 


with the solution y’ = y’ + w’E, where 


EEE us 
Cae Sine 


Pl = — (a + ib) 


(15) 


A) 


(16) 


It is seen that p’™ decreases exponentially with 
distance from the western boundary and y" 
+ y'W is therefore the western boundary 
solution; similiarly p’° + w’F is the eastern 
boundary solution; whereas y‘ is the solu- 
tion in the central portion of the ocean. 

The boundary conditions (5) along AB are 
already satisfied. The constants a, b, c and d 
can be evaluated from the boundary condi- 
tions along ABCA. Along the western bound- 
ary & = o, and equations (4 a, b) give 


| = 
SIC 


—CtanO+d+ta=o, 2y/p+a—V3b= 0. 
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Along the eastern boundary & = r’/y, where 
r = ru’ (1—Ë/x) is the (non-dimensional) 
latitudinal width of the ocean, and r’, = 
(tan @ + tana) the width at £ = o. The east- 


ern boundary conditions give 
—r'—ltanO+d+c=o, —y/pt+c=o. 


The boundary conditions can be combined to 
yield 


a=—r' (I—e, c=r'é, 


Vb=—r' (1— 3 &), d=r (1—e) + £tan ©, 
(17) 


where e=y/pr’ is a small quantity. The 
approximation involved in the non-homo- 
geneous term (7) can now be estimated by 
performing the operation indicated in equa- 
tion (6) on p’°: 


2 No x) y = sin y’ [ty d—x')]. 


The exact solution requires that the right side 
equal sin y’. Since d—x’=r’ (1—e)— y, or r’ 
at most, the neglected term is less than y?r’, 
that is less than 1 per cent. 

The solution becomes 


5 io) 
’ ar ve 3 ee 
y =F smc me 1 =¢elees pcos 
> 


(18) 


For north-south boundaries, 9 = 0, p = 1, 
and the equation reduces to the rectangular 
case discussed in (WOC). The first term in the 
bracket gives the variation of y with x near 
the western boundary, as shown in fig.2. 
It indicates a series of exponentially “damped” 
boundary vortices, with a strong north-Howing 
current and countercurrent in the western 
vortex, both of latitudinal width 27/3. pk. 
The principal effect of the inclined western 
coast line is to widen the western current 
(measured normal to the boundary) by the 
factor 


p-! cosO = (sec O)13 
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ce) 5 10 15 


WESTERN : COUNTER 
CURRENT CURRENT 


WESTERN VORTEX 


Fig. 2. East-west variation of y’ near the western bound- 

ary for an infinitely wide ocean [equation (18) with 

€= 0]. The arrows indicate the direction and intensity of 

the long-shore component of the current. Compare 
with Figures 4 to 6. 


and to reduce the current intensity by a similar 
factor. The approximation refers to low lati- 
tudes; in high latitudes the modification associ- 
ated with an inclination of the coast line is 
reduced, or even reversed. The total transport 
remains approximately thesame. The inclina- 


WIND STRESS 100° 140° 180° 
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tion of the eastern boundary is important only 
insofar as it involves the width r of the oceans. 


Circulation of the Pacific Ocean 


Fig. 3 shows a chart of the Pacific Ocean on 
which north-south and west-east distances 
have been preserved. The basin north of 
15° N is approximated by an isosceles triangle 
with a base of 13 units and altitude of 4 units 
(T unit = 1032 km, see fig. 1). A rectangle 
is also drawn, which has the proper width at 
a latitude of 32.5° N midway between the 
easterly and westerly winds. 

Figs. 4—6 show the computed stream lines 
in the triangular and rectangular ocean basins 
for the following numerical values of the 
parameters: 

Dives Fig. 6 
rectangle triangle 
DE Calor 108 

nis? lach gi senlat: 


Fig. 4 
Shape of basin triangle 
vaViunGane GX! 2752102 
zn 332ulat. 


At 30° N latitude: 


k in cm1 20710 1120 107 11.25 - 10" 
y .040 0.040 .064 
p 445 1.00 445 

DIE. 

— in km 410 182 650 
V3pk 

270 : 

—— cos Oinkm 215 182 341 


V3pk 


140°W 100° 


6000 


4000 


2000 


8000 


km.W of 180° 


4000 6000 
km. E of 180° 


Fig. 3. Chart of Pacific Ocean. The outer scales give distances in km west or east of 180°, and north of equator. 
The circulation induced by a zonal wind stress distribution (shown to the left of the figure) has been computed 
for the cases of the triangular and rectangular oceans drawn in the figure. 
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Fig. 4. Computed transport stream lines in the triangular ocean basin for A= 2.5 ° 


The lines 


Io? cm? sec, 


represent equal values of (ß,/]')y. The curve to the left represents the variation of y with latitude. Unit distances 
along the x- and y-axis equal 1032 km. The section between x = 1.5 and x = 4.0, y= 0.5 and y = 2.0 has 
been redrawn on an enlarged scale in Figure 8. 


The last two lines give the latitudinal and 
actual widths of the western current and coun- 
tercurrent at the latitude 30°N where these 
currents attain their maximum intensity. 

The stream lines are drawn! for equal values 


of 


Poy == are SI) 


[F(pé, ©] (19) 


where B, is the value of B at y = o. The func- 
tion in the bracket is the one in equation (18). 


1 The approximations involved break down near the 
northern vertex of the triangle where the 2 boundary 
strips overlap. In this region the stream lines have been 


The expression preceding the bracket is a 
function of y only, and shown to the left of 
figs 4 and 5. In the case of the triangular 
ocean this function vanishes at latitudes 15° N, 
50° N and 60°N, and accordingly divides 
the circulation into an anticyclonic and a 
cyclonic gyre. These gyres have been named 
the subtropical and subpolar gyres (WOC), 
and correspond roughly to the regions of the 
subtropical anticyclones and of the cyclonic 
storms, respectively. In the case of the rec- 


estimated. Because the real and abstracted boundaries 
differ greatly in this small region, it seems hardly worth- 
while to devise special methods to compute the exact 
position of the theoretical stream lines. 
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Fig. 5. Computed transport streamlines in the rectangular ocean basin for A = 2.5: 107 cm? sec-1, 
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Fig. 6. Computed transport stream lines in the triangular ocean basin for À = 108 cm? sec—}. 


tangular ocean only part of the subpolar gyre 
is included. 

The location and intensity of the gyres on 
figs. 4 and 6 compares favorably with oceano- 
graphic measurements, as presented in fig. 7. 
Note that the maximum transport of the 
western (Kuroshio) current takes place at 
30° N latitude; that the ratio of mass transport 
of the Kuroshio Current to that of the Oyashio 
Current is of the order of 10:1. The eastern 
Pacific gyre and the California Current are 
the result of local meridional winds (WOC) 
and are therefore not present in our model 
for zonal winds. 


Discrepancy between observed and computed 
transport 


Computations based on the available wind 
data and the best available information con- 


crie 
SE 


cerning the relationship between wind speed 
at anemometer level and the stress exerted 
on the sea surface gives [’= 0.65 dynes cm”?. 
(WiOG), Setting. $5: = 222110 ema! seca! 
one finds that each unit of ß, y/l'as drawn on 
figs. 4—6 represents a transport of 3 -10!? g 
sec-!, or 3 million metric tons sec-1. This 
is one half the contour interval in fig. 7. We 
conclude that the computed pattern of water 
transport resembles the observed pattern quite 
closely, whereas the computed quantity of 
transport is about half the observed quantity. 

This discrepancy cannot be ascribed to the 
value of the eddy viscosity, for, according to 
equation (18), the transport is almost inde- 
pendent of A. Nor can the discrepancy be 
ascribed to our assumption of no slippage 
against the western boundary. Taking the 
extreme case of free slippage against the western 


Vn 


Fig. 7. Transport chart of the North Pacific as derived from oceanographic observations. Transport between 

adjacent stream lines are six million metric tons sec”!. The stream lines are based on Fig. 205 of Sverdrup et 

al (1942), the arrows on the “Weltkarte zur Übersicht der Meerestromungen”, Tafel 22, Deutsche Seewarte, 
Ann. d. Hydrogr. u. Mar. Meteor., 1943. The projection is the same as the one in Figure 3. 
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Fig. 8. Computed Pacific “Sargasso Sea” circulation in 

the triangular ocean basin for A = 2.5 +107 cm? sec7}?. 

This figure is an enlargement of the south-west portion 
of figure 4. 


boundary by replacing boundary condition 
(4b) with #y/dv?=o at £ =o, yields 
b= + r/V3 (instead of b = —r'/Y3), and 


Y'max = 1301 sind at p& = 4n/3 V3 
instead of 


Y'mx = 1177 sin ar p& =a2n/ V3. 


Thus with maximum velocity occurring di- 
rectly against the boundary the western current 
has two-thirds its former width, and the trans- 
port is 20 per cent higher. 

A more important modification may be 
related to the effect of the earth’s rotation on 
lateral stress. Following the customary proce- 
dure we have, from analogy with the Navier- 
Stokes equations, introduced lateral stress into 
the integrated equations of motion through 
the term AV?M. (WOC) Since the Navier- 
Stokes equations were derived for a resting 
coordinate system it would seem appropriate 
to replace the relative mass transport M by 
the absolute value M + iQRh cos y, where 
iis the unit vector pointing eastward, and h 
the effective depth of the current. Perform- 
ing the curl operation! one obtains an addi- 
tional term 


a AQh 
om (QRh COS p) = Re 


to be added to the wind curl on the right side 
of equation (1). Setting A = 2.5 - 107 cm? sec-t, 
= 30°, h = 105 cm, gives 2-r07° 9 sec? 


1 Here again we employ a plane coordinate system. 
In spherical coordinates the modification is several 
times larger. We hope to return to this problem in greater 
detail when considering a basin on a spherical earth. 
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compared to a wind curl of about 9-10 g 
sec-?. The effect would be to strengthen the 
transport in the subtropical gyre by about 
25 per cent, to weaken the transport in the 
subpolar gyre, and to displace the boundary 
between the gyres northward. 

The physical basis of this modification 
forms the substance of RossBy’s (1936) dis- 
cussion concerning the merits of Taylor’s 
vorticity-transport theory and Prandtl’s mo- 
mentum-transport theory. If angular momen- 
tum is to be conserved, then any meridional 
interchange of eddies must induce a relative 
current setting eastward at high latitudes, and 
westward at low latitudes. 

The discrepancy may also be ascribed to an 
underestimate of the surface stress for a given 
wind speed, especially at wind speeds of Beau- 
fort four or less (WOC). 

These remarks apply particularly to the 
western currents. In the case of the north 
equatorial current (of which only half is 
shown on figs. 4—6), the equatorial counter- 
current, and the south equatorial current, 
SVERDRUP (1947) and REID (1948) have ob- 
tained satisfactory agreement between the 
transport computed from the wind stress 
and that computed from oceanographic ob- 
servations. 


The western boundary vortices (Sargasso 
Sea) and the value of 4 


An interesting feature is the double vortex 
system near the western boundary on figs. 4 
and $, which is shown on an enlarged scale on 
fig. 8. In the case of fig. 6 the small eastern 
vortex is absent, but its effect can still be 
noticed by the bulging of the stream lines. 


70° 60° 50° 


Fig. 9. Idealized topography of the sea level of the North 
Atlantic Ocean according to DEFANT (1941). The units 
are dynamic centimeters. 
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Whether the eastern vortex will or will not 
be present depends only on the width of the 
ocean and the value of the eddy viscosity, as 
will be seen from an inspection of equation (18). 
The curve in Fig. 2 represents a graph of the 
function f(p£) in the bracket of equation (18) 
for an infinitely wide ocean (e = o). The 
function oscillates about the value of 1, and, 
strictly speaking, the vortex system extends to 
infinity, though the amplitudes of the vortices 
die off exponentially with distance from the 
western boundary. In the case of an ocean 
of (non-dimensional) width r the function 
oscillates not about the horizontal line f(p €) =1, 
but about the line f(pé) = 1—epé& The 
eastern vortex exists only if the slope of the 
first term in the bracket of equation (18) at the 
distance (p&)m = 14x/V3 = 8.45, where it 
reaches its third maximum, exceeds the slope 
— e of the asymptotic line; that is, if 


> 


RZ Le 65m sin = (09. | = OOS: 
V3 2 


The corresponding conditions on y and A are: 


7 (0.015 pr’)? B 
3% SOOT, BUS lic ar 


For the Pacific at latitude 30° N, p = .445, 1° = 
un do oc sec. +ar/n = 3 000 
km, which gives A < 8.25 107 cm? sec}. 
For the Atlantic at the same latitude, 9 = 50°, 
D ROC, and A < 7.8: 107 em "sec" 
in order for the eastern vortex to exist. 

In terms of the dynamic topography of 
the sea surface the double vortex system re- 
presents two humps in ‘the sea surface, the 
western hump being much larger than the 
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Fig. 10. Stream lines in the north-west Atlantic from the 
surface drift of ships, according to Felber (1934). 
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Fig. 11. Distribution of the sperm whale based on log- 

book records dating from 1761—1920 (according to 

TOWNSEND 1935). Each point represents the position 

of a whale ship on a day when one or more whales were 
taken. 


eastern hump (fig. 8). The two humps are 
separated by a trough running parallel to the 
western current with a saddle point directly 
between the humps. Observations in the 
Pacific Ocean are inconclusive with respect 
to the existence of the eastern hump. The 
corresponding area in the Atlantic, the Sar- 
gasso Sea, is one of the most intensely studied 
areas in the world, but even so the conclu- 
sions are uncertain. DEFANT’S (1941) analysis 
of the dynamic topography based on observa- 
tions from the Meteor expedition, and previous 
expeditions, is reproduced in fig. 9. It shows a 
small hump in the Bermuda area, but accord- 
ing to Defant (personal communication) the . 
data was not quite convincing in this respect. 
At Defant’s suggestion FELBER (1934) had pre- 
viously analyzed the extensive observational 
material based on surface drifts of vessels, 
and his stream lines show a saddle point 
southwest of Bermuda throughout the year. 
His figure for April is reproduced in fig. ro. 

Defant’s and Felber’s analyses are consistent 
with respect to the existence of a saddle point 
at about 70° west, 30° north, but they differ 
regarding the Gulf Stream countercurrent. 
This discrepancy has to do with the different 
methods of averaging employed by the two 
authors. The large number of observations? 
collected since Defant summarized the data 
in 1941 have made increasingly clear the varia- 
bility of the circulation in the Bermuda area, 


1 We are indebted to C. O’D. Iselin, F. Fuglister, and 
A. Worthington of the Woods Hole Oceanographic 
Institute for having made these observations available 
to us. 
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Fig. 12. Dimensions of semi-circular ocean basin. 


and the shifting in position from season to 
season, even from week to week, of the Gulf 
Stream and related features (ISELIN and Fuc- 
LISTER, 1948). Defant’s analysis is based on 
synoptic profiles, each of which reveals a 
countercurrent seaward of the Gulf Stream, 
occurting however at different locations for 
different cruises. Felber’s averaging at fixed 
geographic location would naturally tend to 
obscure secondary features adjacent to’ strong 
but wandering primary features. The differ- 
ences are somewhat similar to the differences 
between synoptic weather maps and climatic 
charts. 

No suitable methods have yet been de- 
veloped for combining results from different 
cruises to bring out the existence of secondary 
features such as the countercurrent and the 
eastern vortex. For this reason observations 
collected since 1941 have not added materially 
to our understanding of the Sargasso Sea 
circulation, but most of these observations 
indicate a north-flowing current west of Ber- 
muda. A curious clustering in the distribution 
of sperm whales may perhaps be interpreted as 
evidence for the existence of the two vortices 
and the trough. Each point in fig. 11 repre- 
sents the position of a whale ship on a day 
when one or more whales were taken, accord- 
ing to logbook records dating from 1761 to 
1920 (TOWNSEND, 1935). Relatively few points 
are found in the region of the trough, yet 
there appear to have been no reasons for 
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whaling ships to avoid- this region. These 
troughs, according to Felber’s presentation, are 
regions of convergence. In such regions the 
concentration of nutrients is small, and one 
may expect a relatively small concentration of 
plankton, and of sperm whales which feed 
on plankton. It is not impossible that the two 
vortices represent optimum conditions for the 
existence of plankton, these being the areas 
where freely floating organisms would re- 
main for the longest time under uniform en- 
vironmental conditions. 

In a comparison of figs. $—8 one must 
keep in mind that fig. 8 applies to the Pacific 
Ocean, whereas figs. $—7 apply to the Atlantic. 
As a whole, the evidence, though far from 
conclusive, seems to support the existence of 
the eastern vortex, and would therefore indi- 
cate that the lateral eddy viscosity in this 
region is less than 7.8 : 107 cm? sec-1. 


The effect of curvature in the boundaries 


It has been stated previously that the present 
method can be extended to an ocean basin of 
arbitrary shape provided the boundaries are 
not too irregular. To express this in quanti- 
tative form, consider a semi-circular basin of 
diameter (not radius) r’, (fig. 12). For a 
coordinate system we choose 


5 
S 


where, as previously, y& is the non-dimensio- 
nal distance from the western boundary, @ its 
inclination relative to a north-south direction, 
and r’, is the width of the basin at y =o. 
Differentiation leads to the same transforma- 
tion as given by (10) except that © now de- 
pends on ¢. Neglecting variations in y, we 
obtain 


Fo JEA 


72 > 72 
er sed — ) 


The second term is a.correction for the curva- 
ture of the coastline. This correction is small 
provided the radius of curvature 


I 
= 19 > y sec ©. 


As an extreme example of irregular bound- 
aries consider the sharp change in the orienta- 
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tion of the eastern scaboard between Jackson- 
ville, Florida, and Savannah, Georgia. Setting 
© = 0.04 would require the radius of curva- 
ture to be large compared to .04 units of 
distance, or 40 km. The actual radius of cur- 
vature is 160 km. 


Conclusions 


The present method appears to be suitable 
for a numerical computation of the gross 
features in the ocean circulation, taking into 
account the outlines of the actual boundaries. 
The simple models discussed here give an 
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indication of the results of such a computation. 
Changes in the orientation of the coast line 
will lead to a widening or narrowing of the 
western current, but the general features 
should be more or less the same, irrespective 
of the boundaries. This is in agreement with 
the fact that the general circulation is similar 
in the various oceans, whereas the ocean 
boundaries are not. 
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Abstract 


Observations indicate that the density of ions increases during meteorologically undisturbed 
summer nights. This typical increase is discussed and it is suggested that the probable cause 
is to be found in the accumulation of emanation exhaled from the ground. The accumula- 
tion occurs in weather situations which weaken or suppress convection. 


Introduction 


It is well known that the density of ions in 
atmospheric air is highly variable. This fact is 
shown also by measurements of the variation 
characteristics of the density of small-ions 
carried out at the Institute of High Tension 
Research of the University of Uppsala (No- 
RINDER a. SIKSNA 5, 6). This variability is in 
part due to the diversity of basic processes 
which lead to the creation and annihilation of 
such ions in atmospheric air. It is not always 
possible to isolate these effects in a clear manner 
and they are often superimposed. To these 
processes must be added the fact that the 
character of the variation of ion-density is 
dependent also on local environmental factors. 
It appears in fact possible to state that the ion- 
density depends quite strongly on certain 
meteorological factors and on the location 
selected for such measurements. It is not always 
possible to state with certainty the character 
of this dependence, but it is clear that the causes 
of the variation of the ion-density in atmos- 
pheric air cannot be identified until we have 
succeeded in establishing the undisturbed form 
of this variation. 

In the publications referred to above we have 


discussed some influences that have been 
generally accepted by other authors. Certain 
variations, especially the increase of ion- 
density, can be explained through the advec- 
tion of air masses with higher content of ion 
producers; such air masses may be enriched 
with producers of ions by travelling great 
distances over land. The increase of the den- 
sity of negative ions during rainy weather is 
possibly caused by the negative charges pro- 
duced by spraying water. The increase of posi- 
tive ion-density during a snowstorm may per- 
haps be explained through the electrification of 
the snow-particles by collission and shattering. 

A series of observations showing decreasing 
ion-density during winter could perhaps be 
due to the transformation of small-ions into 
heavier ones, and this process may be connected 
with humidity conditions in the air. In the 
winter measurements no relation between 
ion-density and air temperature was found. 
Summarily it can be stated that variations of 
the ion-density during the winter measure- 
ments were mainly determined by certain pro- 
perties of the air-masses advected, such as 
content of ion producers, humidity etc. 


VARIATIONS IN THE DENSITY OF SMALL-IONS 


In this article we will discuss a type of ion- 
density variation observed by us during the 
summer of 1949 which is especially typical of 
certain meteorological conditions. 


Variation on the ion-density during 
undisturbed summer nights 


During the first seven days of the uninter- 
rupted series of summer measurements a kind 
of variation was observed which was regularly 
repeated, and which could not possibly be 
explained in terms of air-masses advected 
from other regions. The form of these varia- 
tions is as follows: 

During the late hours of the afternoon or 
in the evening the ion-density gradually in- 
creases, reaching a maximum value of more 


than 3 000 n/cm? at about 4 AM the following . 


morning. During the following 3—7 hours 
the ion-density decreases to the average diurnal 
value of c. 1 200—1 500 n/cm?. Such varia- 
tions of the ion-density are illustrated in fig. 
1—8. The positive ion-density is denoted with 
heavy full lines (—) and the negative ion- 
density with thick dashes (— — —); the time 
used is Central-European. In order to indicate 
the meteorological conditions favourable for 
this effect some of the meteorological factors 
are also shown in these figures. The meteoro- 
logical data were obtained from the Meteoro- 
logical Institute of Uppsala University. The 
direction of the wind is indicated at the top of 
each figure for every hour, but it is marked 
only at the points of change so that for a 
particular hour the previous designation is 
valid if no designation is given.. The velocity 
of the wind (m/sec) is also represented graph- 
ically in the figures with thin full lines ( ). 
The cloudiness is represented between the two 

arallel lines at the top of each figure; for a 
cloudless sky the space between the lines is left 
blank, for a clouded sky it is filled. Rain is 
designated by / / / /. In addition to the ion- 
density the temperature of the air in °C is 
given by thin solid lines ( ), the water 
vapor pressure in mm Hg is denoted with thin 
dashes (— — —), and the relative humidity 
in % with thin dashes and points (— - — - —). 

From the figures it appears that two of 
these meteorological factors are related to 
the observed increase of ion-density during 
night-time. These factors are the temperature 
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of the air and the velocity of the wind. The 
effects of water vapor will be discussed later. 
The maximum ion-density very nearly coin- 
cides with the minimum temperature. When 
the ion-density increases the velocity of the 
wind decreases gradually, and when the ion- 
density decreases again, the velocity of the 
wind increases. 

It seems most probable that the temperature 
and the velocity of the air are the factors that 
are responsible for the variations of ion-density 
here discussed. From this a probable mechanism 
for the origin of the increased ion-density 
observed during the time period under con- 
sideration may be deduced. The additional 
ions necessary to produce the increase of ion- 
density could not have been conveyed to the 
place of observation from another region 
because the horizontal motion of the air is 
feeble (in some cases calm prevailed). At the 
same time the vertical air motion must have 
been negligibly small, since with minimum 
temperature at the ground turbulence must 
have been almost completely supressed. After 
the lowest temperature is reached, as the air 
again begins to warm up, vertical convection 
appears mixing results and the ion-density 
decreases. The cause of the increased ion- 
density must be found at the place of observa- 
tion itself. 

It is generally accepted that ions in the lower 
layers of the atmosphere are produced mainly 
by cosmic rays, by the ionization of air 
particles through radiation from radioactive 
substances in the earth and from such ones ex- 
haled from the earth’s surface and spread in 
the atmosphere. In our case the most probable 
prime cause of the increased ion-density might 
be expected to be the accumulation in those 
air layers nearest the surface of emanation 
exhaled from the soil. 


Exhalation of emanation from the soil and 
its accumulation in the lower layers 
of the atmosphere 


Direct measurements of the exhalation of 
emanation from the soil have been carried 
out repeatedly. The first attempt to obtain 
such measurements was made by EBERT (1), 
but it was not possible to evaluate his results 
quantitatively. More exact measurements were 
carried out by JoLLY and SMYTH (2) and SMYTH 
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(tr) at Dublin and by Wricur and SMITH (13) 
at Manila, but the techniques used were not 
handy, and much time was required for the 
measurements. A somewhat different method 
was used by Kosmatu (3) at Liebenau. An 
essentially simpler method, similar to the 
one used by Scum (8) to determine the 
contents of emanation in cellar-rooms, was 
proposed by Hess, and measurements of soil 
exhalation using this method have been made 
by Zupancic (14) at Innsbruck. All the me- 
thods listed above have the disadvantage that 
the emanation is collected over a fairly long 
time period; thus the results obtained repre- 
sent averages over a relatively long time inter- 
val. With such techniques it is impossible to 
obtain the diurnal course of the exhaled 
amount of emanation, even though variations 
over longer time intervals may be determined. 
Zupancıc has found variations in the exhala- 
tion with ratios of 1: 100. In summer, and at a 
high ground temperatures, the exhalation 
reaches a maximum. 

The exhaled emanation does not remain in 
the layers of the atmosphere near the ground, 
but will generally be spread through the whole 
atmosphere. An accumulation of the emanation 
in the lower layers of the atmosphere can 
appear, however, under certain meteorological 
conditions. Already WRIGHT and SMITH (12) 
have observed an increase in the amount of the 
emanation in the lower layers of the atmos- 
phere at night. The same thing is found by 
SCHMID (7). MESSERSCHMIDT (4) has shown 
a strongly marked diurnal variation in the 
emanation content of the atmosphere. This 
variation is in opposition to the diurnal varia- 
tion of temperature. 

MESSERSCHMIDT explains the observed max- 
imum in the emanation content of the atmos- 
phere during the morning hours, and the 
minimum observed in the afternoon in terms 
of variations in thermal convection in the 
atmosphere. A decrease of 1° C per 100 m 
(adiabatic lapse rate) is sufficient to originate 
vertical convective motion in the air. As the 
air near the earth’s surface is warmed the 
emanation will be carried up from the ground, 
and in this way the emanation content de- 
creases until the time of maximum temperature 
during the afternoon hours. During the gradual 
cooling after sunset, in the absence of wind, 
the same air mass will remain near the ground, 
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and will be enriched there by emanation. This 
state of attains continues until the addition of 
more surface heat during the following 
morning produces further convection. On 
days with strong radiation heating these 
thermal convection currents may bring warm 
surface-air to a height of several kilometres. 
Even on those days with only a small tempera- 
ture variation, vertical motion of the air can 
occur, and the ground-air enriched with ema- 
nation is thereby removed from the vicinity 
of the surface. 


Application of the proposed mechanism for 
the accumulation of emanation in explanation 
of the increase of ion-density during calm 
summer nights 


We shall now attempt to apply the me- 
chanism just discussed for the accumulation 
of emanation in the lower layers of the atmos- 
phere in explanation of the observed variations 
of ion-density. As before mentioned, in our 
case we observe an inverse relationship between 
the ion-density and the temperature. Also, the 
velocity of the wind was small or calm when 
the higher ion-densities were recorded. This 
was particularly noticeable. during the last 
hours of 10.6 (fig. 1) and the first hours of 
11.6 (fig. 2); during this time the velocity of 
the wind was extremely low, and for a while 
it was even calm; similar conditions prevailed 
during the last hours of 12.6 (fig. 3) and the 
first hours of 13.6 (fig. 4), atlhough we note 
at midnight an increase in the velocity of the 
wind which causes a temporary diminution of 
the ion-density. A similar course of events 
can be noted on 15.6 (fig. 6) and 16.6 (fig. 7). 
On 17.6 (fig. 8), on the other hand, conditions 
are disturbed by rain. Also, the velocity of 
the wind is greater and, the ion-density is low 
and without great variations; the small in- 
creases in the density of negative ions at 4}, 
110, 18% and 23% may be explained through the 
production of negative charges in spraying 
water. Analogous conditions can be found on 
other days during our uninterrupted measure- 
ment series. 

Thus it is shown that during periods of 
increasing ion-density on calm summer nights 
the meteorological factors observed favour an 
accumulation in the surface layers of the 
emanation exhaled from the soil and therefore 
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Fig. 1—8: Variation of the small-ion density besides some meteorological factors. 
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that the mechanism of this accumulation may 
be the one given by MEssERSCHMIDT (4). For 
these reasons we may assume that the observed 
increase of ion-density is caused by an accu- 
mulation of the emanation exhaled from the 
soil. 

Concerning the humidity we may state in 
general that in summer during undisturbed 
weather, the absolute humidity follows a 
trend similar to the temperature, and that the 
relative humidity follows an inverse trend. 
MESSERSCHMIDT has shown that for different 
seasons of the year there is no relation between 
the relative humidity and the emanation 
content, but that, on the other hand, the abso- 
lute humidity and the emanation content run 
parallel courses. MESSERSCHMIDT’S opinion is 
that the relation between the emanation content 
and the absolute humidity is fortuitious, and 
results simply from the fact that the absolute 
humidity varies in the same sense as the 
temperature. We shall accept this opinion, 


since when the increase of ion-density is 
caused by the accumulation of exhaled emana- 
tion, the transformation of newly formed ions 
into heavier ones cannot take place very 
rapidly, and the absolute humidity cannot 
therefore play an important röle. 

It cannot be claimed that the phenomenon 
here considered must appear everywhere, 
whenever the’ meteorological conditions are 
appropriate. In our case special circumstances 
can be pointed out which may favour the 
occurrence of the phenomenon. The investiga- 
tion of radioactive substances in the air carried 
out at the Institute (SIKSNA 9, 10) has shown 
that especially in the cellar-rooms of the In- 
stitute building are high contents of emana- 
tion and its decay products found. The source 
of these radioactive substances must be searched 
for in the soil surrounding the Institute; thus 
it may be assumed that the action of these 
sources may be of considerable importance 
also in the open air. 
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The Circulation of the High Troposphere over China in 


the Winter of 1945—46 
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A Abstract 


The paper studies the structure of the high troposphere over China in the autumn of 1945 and 
winter Of 1945—46. The principal feature is the existence of two belts of maximum westerlies, 
one flowing around the southern and the other around the northern edge of the Tibetan Plateau. 
The southern jet stream (the main one) is extremely stationary in position. Its speed increases 
downstream beyond the edge of the Asiatic continent. The onset of this southern jet stream is 
abrupt in the middle of October over central and southwestern China and advances down- 
stream at a speed of about 3° longitude per day. 

South of this main jet stream is a belt of extremely uniform absolute vorticity which is zero 
in December and January. In spite of the existence of zero absolute vorticity the circulation 
above 15,000 feet is strikingly stable. Below 10,000 feet is the regular procession of warm 
troughs and cold ridges. 

Comparison with the conditions along 76° E and the east coast of United States has also 
been made. The difference in the structure of the basic current along the east coasts of the two 
continents suggests an explanation of the observed formation of winter typhoons off the Chinese 


coasts versus the non-existence of such storms in the Caribbean Sea in the cold season. 


In the past, our knowledge of the upper-air 
circulation in the Far East has been limited 
to the low troposphere. CHu (1934) and Tu 
(1939) first attempted to describe the wind 
structure in this region. More recently, Kao 
(1948) studied the upper winds up to 15,000 
feet with the aid of pilot-balloon observations 
of 28 stations in China and adjoining coun- 
tries. He published monthly streamline charts 
at 10,000 feet. 

In the latter part of World War II, espe- 
cially during 1945, the Armed Forces of the 
United States maintained numerous radio- 
sonde and rawin stations in China (fig. 1). 
The observations from these stations make 
possible an extension of the previous investiga- 
tions to include the upper troposphere. This 
report will describe the structure of the upper 


1 Present address: Geophysical institute, Nanking. 


troposhere during the autumn of 1945 and 
during the winter of 1945—46. In addition 
to the rawin data, use will be made of 300-mb 
charts prepared at the University of Chicago. 
Although the observations are limited to one 
season, it will be shown that the principal 
features encountered probably are valid in 
all years. 


Onset of winter circulation 


Figs. 2 and 3 show the changes in the wind 
field at Kunming and Chihkiang during Oc- 
tober, 1945. The striking feature at both 
stations is the sudden appearance of strong 
winds at high levels about the middle of the 
month. Prior to October 13—14, winds in 
excess of 50 mph are not observed at all and 
most of the reported speeds are much less 
than this (Kunming about 20 mph, Chinkiang | 
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110 120 
Fig. 1. Distribution of rawin stations. The dotted line 
at the upper left corner outlines the contour of 3 km 
height of Tibetan Plateau. 
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30 mph). After this date wind speeds above 
25,000—30,000 feet consistently exceed so 
mph and frequently reach 75 mph and more. 
In contrast, there is little or no change at low 
levels. Such increase, as can be observed, lags 
in time compared to the strengthening of 
the flow in the high troposphere. 
Coincident with this abrupt change, the 
daily 300-mb charts showed a distinct differ- 
ence in pattern prior to and following the 
middle of October. In the earlier part of this 
month, the contour spacing was wide and the 
direction of the contours varied irregularly 
with the passage of troughs and ridges aloft. 
Later there was a pronounced crowding of the 
contours south of latitude 30° N in the region 
where the rawin observations had shown the 
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large increase in wind speed. The shape of 
the contours was such as to indicate that 
suddenly a strong upper flow circling the southern 
rim of the Himalayan Plateau had set in. 

Fig. 4 is a typical example of the flow 
pattern at 300-mb after the beginning of the 
winter circulation. Although the contour 
configuration over the Himalayas is uncer- 
tain as there are no stations, the adjoining 
observations indicate with reasonable certainty 
that the westerly current is split into two 
branches, one situated north and the other 
south of the mountains. This map pattern is 
representative for the entire winter season. 
It is of interest that the width of the southern 
current decreases downstream (northward 
displacement of subtropical ridge) and that 
the maximum wind speed is reached far 
beyond the boundary of the high plateau. 
Two lee eddies are shown, one cyclonic and 
one anticyclonic. Both the data of November 
4 and continuity serve to support the analysis 
as drawn. Eddies that develop in the lee of 
the mountains and then travel eastward, were 
observed on numerous occasions. In partic- 
ular, formation of cyclonic centers near 
longitudes 100°—105° E appears to be a suffi- 
ciently regular event that it is possible to 
speak of a mean cyclonic circulation in this 
region. 

The rapid change of circulation pattern in 
the middle of October provides for an in- 
teresting comparison with the analysis of the 
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Fig. 3. Vertical time-section at Chihkiang, October 1945. The representation of the wind speed is the same as in fig. 2.. 


onset of the summer monsoon over India 
and Burma made by Yın (1949). He indicated 
how a small northward displacement of the 
upper westerlies over Eurasia could lead to a 
rapid transition from winter to summer 
flow pattern south of the Himalayas as the 
westerlies suddenly began to circle around 
the northern rather than the southern periph- 
ery of the plateau. The present data suggest 
an equally abrupt change in the opposite 
direction at the end of the summer period. 
In this connection it is of interest to note 
that Chaudhury! observed that the first pro- 
nounced invasion of the middle-latitude west- 
erlies into northwest India during October, 
1945 took place roughly one week earlier 
than shown in figs. 2—s. 

This eastward progression of the onset of 
the winter circulation can be traced farther 
downstreain. At Okinawa a rapid increase of 
the westerly flow above 25,000 feet occured 
on October 21 and at Iwo Jima as late as 
October 27. Thus we observe a rather regular 
downstream progression of the beginning of 
the winter season aloft over eastern Asia 
during 1945. The strong westerlies advance 
from longitude 70°E to 140°E in 20 days, 
or at a rate of 31° long. day". This is much 
less than the upper wind speed itself. The 
propagation of change from summer to 
winter conditions differs from what may 
1 A. M. CHAUDHURY, manuscript, University of Chi- 
cago, 1949. 


have been expected, namely a simultaneous 
southward intrusion of the westerlies at all 
longitudes. It is plausible to suggest that the 
downstream propagation is peculiar to eastern 
Asia and that it is produced by the effect 
of the Himalayas on the upper circulation. 


The zonal wind distribution, December 1945 
and January 1946 


We shall now turn to the structure of the 
upper wind field in the middle of the winter 
season. The discussion will treat the zonal 
wind component, since the meridional flow 
is very weak above 15,000 feet. Fig. 5 gives the 
regional distribution of zonal flow. A distinct 
maximum is in evidence in southeastern China. 
The existence of this maximum is consistent with 
the presence of a narrow but intense branch of 


Fig. 4. 300-mb map, 4 November 1945; 0400Z. 
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Fig. 5. Distribution of the mean zonal wind speed in 

mph (solid lines) of December 1945 and January 1946 at 

40,000 feet and mean precipitation pattern (dashed lines) 

in winter. The figures are the mean speeds and the heavy 
line shows the position of the jet stream. 


the westerlies that flows around the southern border 
of the Himalayas. As seen from the down- 
stream increase of wind speed, this current 
intensifies even after leaving the vicinity of 
the mountains. Its center does not closely 
adjoin the mountains but lies at some distance. 
Wind speeds in the immediate lee appear to be 
quite low. This region corresponds to the 
stagnation point in the potential flow around 
a smooth cylinder. Since the atmosphere is 
not an ideal fluid and since its motion is not 
irrotational, it is very likely that lee eddies 
of small diameter exist in the region of stagna- 
tion as shown in fig. 4. There may even be 
light easterlies close to the mountains on the 
average. 

The observation at Peking requires special 
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mention. This station lies sufficiently far north 
to be affected by the westerly current that 
circles the Himalayas on the north side. Evi- 
dence for the existence of this current will 
be given later. On the other hand, the resultant 
wind at Peking is obtained from a small 
number of observations some of which appear 
very doubtful. For this reason, the analysis 
on fig. 5 has not been extended to northern 
China. | 

Wind steadiness: The preceding discussion 
has indicated that the winter-time westerlies 
aloft over central and southern China are 
derived from a branch of the circumpolar 
westerlies that flows around the southern edge 
of the Himalayas. If this is true, we should 
observe a high degree of steadiness of the 
upper flow. Such steadiness previously has 
been noted over Argentina in the upper 
circulation that crosses the Andes Mountains 
(Born, 1948'; RIEHL, 1949). Wind steadiness 
(S) in percent, is defined by the formula 


R 


S =. 


100, 


where R is the speed of the vector resultant 
wind and V the average speed regardless of 
direction. Table 1 gives the wind steadiness 
as a function of height for the rawin stations 
for all available periods and fig. 6 shows the 


1J. A. Born, “Effect of the Andes Mountains on the 
General Circulation over the Southern Part of South 
America’, unpublished part of the M. S. thesis, Uni- 
versity of Chicago, 1948. 
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vertical time section of upper winds at Han- 
kow, January 1946, as a synoptic example. 
Clearly, the wind steadiness at high levels is 
very great, in support of the suggestion on the 
origin of the upper current. At lower heights, 
the steadiness is considerably less, as can be 
seen from the table as well as figs. 2, 3, and 6. 
Near the ground a procession of troughs 
and ridges moves eastward. These disturb- 
ances die out with height and can hardly be 
observed at 15,000 feet. Therefore the low- 
pressure areas must be warm, and the high- 
pressure areas cold, relative to their surround- 
ings. This same procession of thermal disturb- 
ances can also be noted farther east at Oki- 
nawa and Iwo Jima. 


Table 1. 
The wind steadiness as a function of height 
for the stations in China proper 
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It is surprising that the wind steadiness at 
2,000 feet, well within the layer of monsoonal 
flow, is as low as given in the table. Chihkiang 
is the only station that conforms to the pattern 
that was expected: high wind steadiness near 
the ground and in the upper troposphere, with 
an intermediate region of low steadiness. 

Jet stream and rainfall distribution: The fore- 
going data on wind steadiness aloft, together 
with the indicated role of the orography in 
producing the upper current, suggest that the 
wind distribution noted in the winter of 
1945—46 should hold also for other years. 
It is of interest, therefore, to compare the 
upper flow pattern with the mean precipita- 
tion pattern of winter (Fig. 5, U. S. Weather 
Bureau’). A maximum of precipitation practi- 
cally coincides with the center of the upper jet 


1See E. R. Bret, Climate of China, report to U. S. 
Army Air Service. 
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Fig. 7. Mean number of days with rain in October. 


stream. Rainfall decreases northward and south- 
ward from there, but much more rapidly on the 
south side. This rainfall distribution is exactly 
like that found by STARRETT (1949) relative 
to the upper jet stream in the United States. 
It agrees with the distribution to be expected 
from the vertical circulation cell postulated 
for the jet stream region (University of Chi- 
cago, 1947). Such remarkable agreement 
between the precipitation relative to the jet 
in the United States and the average regional 
rainfall pattern of China relative to the average 
position of the jet over China is possible only 
if the latitudinal fluctuations of the jet stream 
over eastern Asia are restricted to very narrow 
limits, not only from one day to the next, 
but also from one year to the next. 

It is of interest to add that charts showing 
the number of days with rain also exhibit a 
very pronounced maximum that runs east- 
west along the upper jet. In autumn, at least, 
there is also a secondary minimum of number 
of days with rain to the south of the upper jet 
stream (fig. 7) and, as computed by Liu’, a 
region of minimum rainfall variability ex- 
tends along the jet stream axis. 

Meridional cross section: Seven of the rawin 
stations operative during December 1945 and 
January 1946 were sufficiently close to the 
meridian 120° E that they could be projected 
on this meridian following the mean upper- 
air flow above 15,000 feet, which was due 
west. It was possible then to construct a ver- 
tical space cross section of the distribution of 
west wind speed (fig. 8), although the data 
are not as homogeneous as could be desired. 
The number of observations at the different 


2 see footnote 1. 


178 


(1000 Ft.) 
45 


i in 
PEKING SHANGIHA! T'YINAWA CANTON LEOAG MANILA TACLOBAN 


Fig. 8. Meridional cross section of mean zonal wind in 
mph of December 1945 and January 1946. 


stations ranges from 40 to 55 for the two 
months’ period. Moreover, the frequency of 
observations decreases with height, so that 
several stations have only 10—1$ reports at 
40,000 feet. Finally there is a large gap in 
the important region between Shanghai and 
Peking where the analysis as inferred has been 
entered with dashed lines (cf. also figs. 5, 13). 
Nevertheless, the principal features of the 
cross section should be representative in view 
of the high degree of steadiness of the upper 
wind field stressed before. 

Fig. 8 shows two centers of west wind 
concentrations. The position of the northern 
maximum is not given by the data presented. 
It can be inferred from the analyzed daily 
300-mb charts mentioned in the introduction 
which show that the center of the westerly 
current passing over Siberia, reaches the Pacif- 
ic coast between Peking and Vladivostok 
(cf. also fig. 9). The strength of this current 
at Peking may not be as great as indicated. 
As mentioned in the discussion of fig. 5, 
some of the Peking data are quite doubtful. 
All that can be said with confidence at this 
time is that there are two westerly maxima. 

The position of the southern jet stream 
center is closely given by the data of fig. 8, 
together with figs. 5, 9 and 10. Its intensity 
can be estimated as about 150 mph at the 
level of strongest wind (highest measured 
average wind 120 mph at 35,000 feet at Shang- 
hai). The vertical slope of the axis of the south- 
ern jet apparently is northward, since the lower 
and middle-troposphere wind speeds are high- 
est at Okinawa. However, this feature may 
not be entirely realistic, since Okinawa lies 
a few degrees east of the main group of sta- 
tions and the intensity of the jet stream in- 
creases downstream (fig. 5). 

On the south side of the southern jet the 


IP Oy ONSUSINNG, Ve le Is! 


~~ 


isolines of equal wind speed are almost vertical 
above the lower levels; a high degree of 
barotropy must prevail in this region. The 
subtropical ridge line intercepts the ground 
near 30° N. At first, it slopes strongly toward 
the equator with height until it reaches lati- 
tude 16° N at about 15,000 feet. From there 
on, the ridge line (base of the polar wester- 
lies) is almost vertical, and it becomes entirely 
vertical above 18,000 feet at Manila. 

Comparison with conditions at 80° W and 
76° E: It is of interest to compare the upper 
wind field at 120° E as given by fig. 8 with 
conditions at 76° E (CHAUDHURY, 1950) and 
80° W (Hess, 1948). In so doing, we again 
make use of the earlier contention that our 
data are representative of mean winter condi- 
tions. It is to be noted that the cross sections 
of Chaudhury and Hess are both based on 
computations from pressure data with use of 
the geostrophic assumption, whereas fig. 8 
is constructed from actual wind observations. 

Outside the tropics fig. 8 is in good agree- 
ment with the section of Chaudhury. At 76° 
E one principal west wind maximum appears 
just south of the Himalayas and a second 
center is indicated north of the mountains. 
These jets correspond to those shown in fig. 8. 
Chaudhury has a third maximum in the far 
south, just where the subtropical line is indi- 
cated at 120° E. An explanation of this differ- 
ence is not offered. 

Comparison of the two east-coast sections 
shows that the heights of the principal jet- 
stream centers are roughly equal. The speed 
of the North American jet is somewhat small- 
er than that of the Asiatic jet, presumably 
because the latitude of the strongest wester- 
lies undergoes considerable changes over the 
United States whereas it is steady over China. 
Both sections show a secondary maximum in 
the north. But the separation of the two 
centers over North America is not as clearly 
marked as over eastern Asia, probably again 
on account of the variability of the jet stream 
position. This variability can also be seen 
from the fact that the isolines of equal speed 
have a lesser slope with height in the lower 
latitudes at 80° W compared to 120° E and 
that the intensity of the anticyclonic shear is 
much smaller (cf. also fig. 13). Finally, the 
subtropical ridge lies 4—5° latitude farther 
south over North and Central America and 
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its slope does not become vertical in the trop- 
ics as over the western Pacific. 

Figs. 5 and 8 and Chaudhury’s section, taken 
together, finally show that over Asia the air 
in the troposphere mainly flows around the 
mountains. In contrast, it largely moves over 
the North American barrier. This suggests that 
the mountains exert different dynamical effects 
on the upper-air currents of the two regions 
which must be allowed for in theoretical ana- 
lyses and in estimates of downstream influences 
such as positions of troughs and ridges. 


Fluctuations of position and intensity of jet 
stream and subtropical ridge 


In order to demonstrate further the high 
degree of steadiness of the upper-air flow 
pattern over the Far East in winter, figs. 9—10 
show the zonal wind distribution at 300-mb 
from October 1945 to January 1946 in synop- 
tic sequence. The west wind speeds given in 
these figures were obtained with use of the 
geostrophic formula from the daily 300-mb 
charts mentioned in the introduction. In the 
longitude intervals 110—120° E and 130—140° 
E the 300-mb height gradient was averaged 
on each day and the geostrophic wind pro- 
files were then combined in successive 5-day 
periods and plotted and analyzed in time- 
section form. Fig. 9 is representative of condi- 
tions in the sectors 110—120° E (center 115° E) 
and fig. 10 for 130—140° E (center 135° E). 

The most remarkable feature seen on both figures 
is the stationary position of the principal jet 
stream center beginning about October 15 at 
115° E and somewhat later at 135° E. Oscilla- 
tions of the center are within very narrow 
limits. The mean deviation from the mean 
position is less than 1° latitude at 115° E and 
less than 2° latitude at 135° E. The west wind 
speed in the central zone is more variable 
and also increases generally from October 
to the middle of the cold season. 

Farther north, the secondary belt of wester- 
lies evident in fig. 8 is present continuously 
over China (fig. 9) in November and De- 
cember 1945, when its mean latitude is 45° N. 
It is less distinct during January 1946, although 
still indicated in the data. This second jet 
cannot be found in fig. ro, and this is due to 
considerable variability of position and inten- 
sity as seen from inspection of the daily charts. 


Fig. 9. Time-section of the 300-mb geostrophic zonal 
wind along 115° E from October 1945 to January 1946. 


Fig. 10. Time-section of the 300-mb geostrophic zonal 
wind along 135° E from October 1945 to January 1946. 


The intensity of the northern maximum at all 
times is less than that of the southern one: 
highest speeds in the north are near 100 mph 
whereas they reach 175 mph in the south. 
This difference holds not only for the 5-day 
but also for the daily values. 

Coincident with the steadiness of the south- 
ern jet stream center we observe that the 
subtropical ridge undergoes only minor fluc- 
tuations in latitude over the tropical part of 
the western North Pacific Ocean. This may 
be seen also by inspection of time sections of 
daily rawin observations, especially at Manila 
and Guam. At both stations the resultant 
wind is nearly zero in the middle of winter 
(cf. also figs. s and 8). This small resultant is 
not produced by alternating high westerlies 
and easterlies which in the average cancel 
each other. Wind speeds in the upper tropo- 
sphere are light on all days. This contrasts with 
summer conditions when a train of large 
upper vortices keeps passing over the region 
(RIEHL, 1948 a). 

The inverse holds over the eastern Pacific 
Ocean, especially near the Hawaiian Islands. 
There the upper winds are most steady in 
summer and most variable in winter when 
large troughs and ridges move over the island 
chain. It is of interest that the rainy season 
occurs in winter in the eastern tropical Pacific 
and in summer in the western tropical Pacific, 
therefore in the season when the high-level flow 
is most unsteady. The dry seasons are associated 
with a relatively steady regime aloft. 

There is also a great difference between the 
upper-air circulation of the tropical West 
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Pacific and the Caribbean Sea. The subtrop- 
ical ridge drops below 10° N east of the Unit- 
ed States whereas it holds at 15° N in the 
western Pacific (cf. fig. 13). Moving troughs 
and ridges in the upper westerlies dominate 
the wintertime weather of the Caribbean. The 
westerlies have a north component and travel 
down the Antilles chain from the Gulf of 
Mexico. Wind speeds often are very high. 
Over Puerto Rico (18° N) WNW winds of 
80 mph are observed frequently and even the 
resultant winds aloft reach 45 mph (STONE, 
1942). 

At this time it is not possible to explain the 
striking difference in the low-latitude flow 
patterns cast of the two major continents. 
It is clear, however, that the existence of 
polar westerlies over the Caribbean prevents 
any possibility for tropical storms to develop 
in that region in winter (RIEHL, 1948b). How- 
ever, because of the persistence of the subtropical 
ridge near 15° N in the western Pacific, dyna- 
mical conditions aloft remain basically favorable 
for storm generation in that area, When the other 
synoptic conditions necessary for tropical cy- 
clogenesis are met, a typhoon can develop. 
Actually, several such cyclones form in most 
winters, as is well known. 


(28°) 


Latitude 


20 40 60 80 


O 


100 120 140 160 mph 


EIS a 
Fig. 11. Monthly mean profiles of 300-mb zonal wind 
along 115°E (a) and 135° E (b) averaged with respect to 
the center of the southern jet stream. 
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Monthly profiles of west wind speed and 
vorticity at 300-mb 


Figs. 11a—b summarize the data of figs. 
9—10 in form of monthly west wind profiles. 
Here, the averaging is done following the 
position of the southern jet stream in order to 
show the wind structure relative to the west 
wind maximum. It turns out that the mean 
latitade of the jet as 28° N at ris E and 
30° N at 135° E. At 115° E the highest speed 
rises from 75 mph in October to 135 mph in 
January. A secondary maximum is evident 
in the north in spite of the method of averaging. 
On the south side of the jet an east wind of 
to mph is found in all four months 15° lati- 
tude south of the center. At 135° E the speed 
of the jet is about 10—1$ mph higher than at 
15° E in all months except October. 

Figs. 12a—b illustrate the distribution of 
absolute vortivity (04) with latitude computed 
from figs. 11 a—b with the formula 

u 


on 

+ -tan®... I 
ro® r ( ) 
Here f is the Coriolis parameter, u the west 
wind speed, r the radius of the earth and ® 
the latitude. Use of this equation is justified, 
since the wind nearly always is due west. 
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Fig. 12. Monthly mean profiles of absolute vorticity, computed from fig. 11a—b, along 115° E (a) and 135° E (b). 


The general shapes of all vorticity curves 
are very similar. There is no abrupt change 
from month to month or from 115°E to 
to 135° E. A region of maximum vorticity 
appears 2—5° latitude north of the jet center. 
In January, this maximum attains a value as 
high as 1.3 x 10~* sect, almost the vorticity 
of the earth at the pole. South of the jet there 
is a broad belt in which the absolute vorticity 
is low and rather uniform. In December and 
January, its value is zero. A sharp gradient of 
vorticity that amounts almost to a discon- 
tinuity, extends across the region of west- 
wind maximum. These features are in accord 
with the description given by other writers 
(cf. PALMEN and NewTon, 1948). The small 
gradients of vorticity north of the jet stream 
are in good agreement with the vorticity pro- 
file postulated by RossBy (1947, 1949). 

Comparison of the four profiles along each 
meridian reveals an interesting change of the 
vorticity distribution with time. As we pass 
from summer to winter, there is a general 
increase of vorticity north of the jet and a 
decrease to its south. This feature probably 
holds on all meridians since the speed of the 
jet increases everywhere from summer to 
winter. Formally, the opposite change of 


vorticity in middle and low latitudes agrees 
with the requirements of the circulation 
theorem, as changes of wind speed in the 
equatorial zone are relatively small. Dynam- 
ically, however, it is a difficult problem to 
explain this change. A mechanism is required 
that is capable of transferring vorticity from 
low toward high latitudes against the gradient. 

Fig. 13 further illustrates the vorticity distri- 
bution with latitude and also compares the 
wind speeds as given by rawin data and 
geostrophic computation. In this diagram the 
latitude is plotted according to the cosine 
scale on the abscissa. According to SOLBERG 
(1936) “dynamic instability” exists in a zonally 
symmetric current, provided that 


We neglect here the effect due to curvature 
of the parallels because we are concerned pri- 
marily with low-latitude regions. Upon inte- 
gration with respect to latitude, 


u = —2wrcos®-+ const.... (3) 


where is the angular velocity of the earth. 
Equation (3) is represented by a straight line 
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Fig. 13. Comparison of mean 300-mb zonal wind of 
December 1945 and January 1946 along 120°E (thin 
solid line) and 65° W (heavy solid line). The abscissa is 
in the scale of cosine of the latitude. The straight dashed 
line is the line of zero absolute vorticity with u = o at 
15° N. The circles are the observed winds at 30,000 feet 
and the crosses indicate the 300-mb geostrophic winds. 


on a diagram with the coordinates of fig. 13. 
In our case, the constant of integration was 
determined by setting u= oat® = 15° N. 

The solid line of fig. 13 gives the wind 
distribution at 30,000 feet as taken from fig. 8. 
Crosses mark the corresponding values of 
geostrophic wind at 120° E as given by the 
300-mb charts. It is seen that actual and 
geostrophic wind values are almost identical 
south of latitude 30° N. The wind profile 
between latitudes 18° and 29° closely approxi- 
mates a straight line, indicating dynamic in- 
stability, in agreement with figs. ı2 a, b. 
It is of interest that this is the latitude belt charac- 
terized by extreme steadiness of the winds. 

Disagreement between actual and computed 
wind is again evident at Peking. Both possible 
solutions of the wind profile have been drawn 
as it is not possible to decide which is the 
correct one. 


For comparison, the corresponding veloc- 
ity distribution east of the United States 
(65° W) has also been entered. It is seen that 


TUSC HEN Gy YEE = 


the anticyclonic shear in this region is much 
weaker than in the western Pacific, that the jet 
stream center lies farther north and that the 
subtropical ridge lies farther south as discussed 
before. 


Summary of the results 


One of the most remarkable features of the 
winter circulation in the Far East is its abrupt 
onset manifested by the sudden appearance of 
strong winds. In 1945 it started abruptly in 
the middle of October over central and south- 
western China. The onset advanced regularly 
downstream at a speed of about 314° longitude 
per day. 

After establishment of the winter circula- 
tion there appear two belts of maximum west- 
erlies, one flowing around the southern and 
the other around the northern edge of the 
Tibetan Plateau. The southern jet stream is 
strikingly stationary in position while the 
northern one varies in latitude from day to 
day. The intensity of the principal jet increases 
from autumn to midwinter. Speeds also in- 
crease downstream beyond the border of the 
mountains and even beyond the edge of the 
continent. The maximum intensity is about 
160 mph in January at 135° E. 

The absolute vorticity south of this main 
belt of maximum westerlies is remarkably uni- 
form. In December and January this uniform 
absolute vorticity is practically zero. 

In spite of the existence of zero absolute 
vorticity the circulation above 15,000 feet in 
this part of the world is extremely stable. 
The wind direction deviates very little from 
the mean, which is due west. In contrast to 
this extremely stable upper-air current are the 
fluctuations of wind direction and speed asso- 
ciated with passage of small disturbances 
below 10,000 feet. Since these disturbances 
rapidly decrease upward in intensity the 
troughs must be warm and ridges cold rela- 
tive to the surroundings. 

Finally we shall mention a difference 
between the upper-air circulation of the trop- 
ical west Pacific and the Caribbean Sea. 
The subtropical ridge drops below 10° N 
east of United States whereas it holds at 15° N 
in the western Pacific. This difference in the 
structure of the basic current in these two 
regions helps to explain the observed forma- 
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tion of winter typhoons off the Chinese coast 
while such storms do not exist in the Carib- 
bean Sea in the colder season. 
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On the Influence of the Earth’s Orography 


on the General 


Character of the Westerlies 


By B. BOLIN, University of Stockholm 


(Manuscript received 15 May 1950) 


Abstract 


The upper westerlies in middle latitudes possess, in the mean, a wave-like character which has 
been explained as a result of thermal contrasts between land and sea. On the other hand, recent 
theoretical investigations by QUENEY, CHARNEY and ELIAssEn have shown that an obstacle of 
the dimensions of the Rocky Mountains generates a wave pattern downstream whose scale is 
comparable to the observed mean waves. In the present paper these theoretical studies are ex- 
tended, and an attempt is made to discuss, in a general way, the influence of the northern 
hemisphere mountains on the character of the westerlies. Finally the paper points out some of 
the climatic consequences of the proposed dynamic-topographic control of the prevailing 


flow patterns. 


I. Introduction 


The rapid increase of upper air observations 
over the northern hemisphere during the last 
few years has made it possible to construct 
circumpolar upper air charts which describe 
the conditions at higher levels in some detail. 
These maps show that the middle latitude 
westerlies are concentrated to a rather narrow 
band (‘‘jet-stream”’) on which large, wave-like 
perturbations are superimposed. The latitude 
of the jet shifts from one day to the next and 
the waves usually move from west to east. 
The formation and the location of the jet- 
stream, as well as the behaviour of the waves, 
appear to be determined by the distribution 
of continents and oceans and by the large 
mountain barriers which are located in middle 
latitudes (the Rocky Mountains and the large 
mountains in the interior of Asia in the northern 
hemisphere and the Andes in the southern hemi- 
sphere). This-is seen on the daily upper air maps 
and is made somewhat more obvious by the 
fact that well-defined perturbations are found 
also on the monthly mean charts for the mid- 
troposphere. The existence of such deviations 


D 
suggests that the influence of the surface of 


the earth is of fundamental importance for the 
character of the westerly flow. 

The location of the waves aloft indicates 
that they are connected, in an intimate manner, 
with the stationary cyclones and anticyclones 
which characterize the mean conditions at 
the surface of the earth in middle latitudes. 
These stationary high and low pressure areas 
have been explained as the result of the tem- 
perature difference between the air over the 
continents and the air over the oceans. The 
thermal influence of the earth has been con- 
sidered to be of conclusive importance. (See 
e.g.: V. BJERKNES and collab., 1933, pp. 686— 
689, SCHERHAG, 1948, pp. 29—30 RE. 
The reasoning is as follows: A solenoidal 
field is established along the coasts. In the 
equilibrium state the main part of this sole- 
noidal field is balanced by the wind change 
with height, but in order to compensate for 
the frictional loss of kinetic energy a direct 
circulation around the solenoids must occur. 


from a straight zonal flow on the mean maps 
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Fig. 1. The mean height of the soo mb surface as a func- 
tion of longitude. In summer the figure represents condi- 
tions in the latitude belt 45° N—so° N, in winter 35° N— 
40° N. The profiles have been computed from the hemi- 
“ spherical mean charts published by ScHERHAG (1948). 


The deviation from geostrophic flow because 
of friction must then be consistent with this 
solenoidal circulation. The picture fitting this 
requirement is well-known: Anticyclonic cir- 
culation at lower levels and cyclonic circula- 
tion aloft around regions which are colder 
than the environment, cyclonic flow in the 
surface layers and anticyclonic flow aloft 
around regions which are warmer than their 
surroundings. This pattern is superimposed 
upon the general westerlies around the hemi- 
sphere. Since the intensity of the zonal current 
increases with height the resulting pattern 
at upper levels should have a wave-like struct- 
ure with troughs over areas of low tempera- 
ture at the surface of the earth, and with 
ridges where the air is heated from below. 
No doubt this monsoonal circulation is an 
important process which effects the general 
character of the westerlies. There are, however, 
several features of the flow pattern aloft which 
seem to indicate that the monsoonal circula- 
tion is essentially a process in the lower tropo- 
sphere. We know that the continents act as 
cold sources and the oceans as warm sources 
in winter while conditions are reversed in 
summer. In spite of this fact certain basic 
features of the wave pattern at upper levels do 
not change essentially from winter to summer. 
This is brought out quite strongly in fig. 1, 
which shows the variation with longitude of 
the height of the soo-mb surface, in winter 
for the latitude belt between 35° N and 40° N, 
in summer for the belt 45° N—so° N. This 
diagram was constructed with the aid of the 
mean circumpolar maps which have been 
published by SCHERHAG. The seasonal displace- 
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ment of the westerlies has been taken into 
consideration through the choice of different 
latitudes for summer and for winter. Even 
though the intensity of the waves around 
the hemisphere is much less in summer than 
in winter the position of the troughs and 
ridges is almost the same. 

Considering now for a moment conditions 
in the southern hemisphere it is well known 
that a pronounced trough is located to the 
east of the Andes Mountains (see fig. 2) in 
Rea the same relative position as 
the North American trough east of the Rocky 
Mountains, in spite of the fact that the thermal 
conditions in the two cases are quite dis- 
similar (cf. Born, 1949). Both these troughs 
occır in winter as well as in summer. These 
facts point to a dynamic effect of the mountains 
rather than to a thermal influence of the surface 
of the earth on prevailing zonal air currents. 

Indications of the type.listed above suggest 
the importance of exploring, in further detail, 
the dynamic-orographic control of the pre- 
vailing upper level flow patterns. In the follow- 
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Fig. 2. Mean flow at 6 km in winter over South America. 
(after BOFFI, 1949). 
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ng, we shall attempt to extend the analysis 
of this control as given by CHARNEY AND ELras- 
SEN (1949), by considering also the north—south 
extent of the mountain ranges and the width 
of the zonal current. 


II. Basic assumptions 


The speed of the westerlies increases with 
height and has a maximum value at the tropo- 
pause in middle latitudes. The greater the 
kinetic energy of a current, the greater must 
be the force required to change its direction 
and character. Thus it seems reasonable to 
begin by considering the layer of maximum 
westerlies and then to assume that the flow at 
other levels (lower troposphere and lower 
stratosphere) is adjusted in accordance with the 
conditions prescribed by that level. At lower 
levels this adjustment must be in agreement 
“with regional differences in heating from 
below and with the restrictions imposed by 
ground friction. Thus we treat at first the flow 
pattern in the layer of maximum westerlies. 

The thermal influence of the ground on the 
atmosphere through radiation and convection 
decreases with increasing height, primarily 
because of two facts: 

(1) The distance from the heat and cold 
sources at the surface of the earth increases 
upward. 

(2) The increase of the westerlies with height 
reduces the time during which an upper 
level air mass remains above a region of pre- 
scribed thermal characteristics. 

This reduced influence of the ground on 
the upper atmosphere is shown by the fact 
that the annual temperature range of the 
upper atmosphere is less affected by the dis- 
tribution of land and sea than the corre- 
sponding range below (FLOHN, 1943). The 
dynamic influence on the large scale flow pat- 
terns should therefore be relatively more im- 
portant at higher levels than at the surface 
of the earth and may even be of greater im- 
portance than the thermal influences at the 
upper levels. In order to see to what extent it is 
possible to explain the wave character of the 
upper westerlies from a purely dynamic point 
of view we shall at first disregard completely 
heating and cooling from below and consider 
only the effect of the mountain ranges. The 
degree of similarity between the theoretical 
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results and the observed pattern will show to 

what extent these assumptions are justified. 
When treating the mean flow pattern by the 

procedure outlined above we have a priori 


made three assumptions, which must be kept. 


in mind. These are: 

(1) The variations from day to day of the 
mean current around the hemisphere are as- 
sumed to be small or at least do not represent 
any considerable deviation from the wave 
pattern computed from the mean zonal wind. 

(2) The effect of a mountain upon the wester- 
lies is assumed to be the same even if the 
current is not straight zonal, but non-station- 
ary troughs are superimposed upon it. 

(3) There exists no critical velocity, at which 
the zonal motion breaks down into vortices. 

We shall return to a discussion of these 
assumptions later on. 

Several theoretical investigations have been 
carried out in which the stationary flow pattern 
caused by an infinite mountain range has been 
computed. QuENEY (1948) has given a summary 
of the different types (spectrum) of waves which 
are generated in an unlimited and infinitely 
broad current striking a mountain range at 
right angles. In our particular study we are 
interested only in waves which are of the 
same order of magnitude as those observed 
in the mean pattern, that is the long planetary 
waves. It has been shown by CHARNEY (1948) 
that other types of waves (gravitational waves, 
inertia waves) may be eliminated from the very 
beginning of the analysis by assuming hydro- 
static equilibrium and geostrophic balance. 
Furthermore, CHARNEY AND ELIASSEN (1949) 
have shown that an equivalent-barotropic at- 
mosphere may be defined whose large-scale mo- 
tion reasonably well corresponds to the mo- 
tion of the real barocline atmosphere at about 
soo mb. We know from experience that the 
flow at soo mb may be considered as repre- 
sentative for the upper troposphere in general. 
By treating the equivalent-barotropic atmos- 
phere we thus get an idea of the flow in the 
layer of maximum westerlies. This is the 
reference level which we have chosen accord- 
ing to the discussion above. 

QUENEY took the vertical stability of the 
atmosphere into account in his discussion of 
the waves generated by a mountain range. 
CHARNEY and ELIASSEN, on the other hand, 
assumed a homogeneous atmosphere with a 
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constant height of the upper boundary and 
this simplified the computations considerably. 
In spite of the great difference between these 
two models of the atmosphere the principle 
features of the large-scale flow patterns are 
the same. This suggests that a fixed upper 
boundary of the homogeneous atmosphere 
in some measure may be interpreted as being 
equivalent to a certain degree of internal 
gravitational stability of the real atmosphere. 
In our analysis we shall follow CHARNEY and 
ELIASSEN and consider a homogeneous atmos- 
phere with a fixed upper boundary. 

As mentioned in the introduction we shall 
extend the analysis from one to two dimen- 
sions in the horizontal plane. Some calcula- 
tions have been made by STEwART (1948) for 
the purpose of determining the flow around a 
mountain with finite extension in the north— 
south direction. He studied the stationary 
horizontal waves generated by a cylindrical 
obstacle. The considerations below are an 
attempt to combine the results obtained by 
CHARNEY and ELIASSEN with those obtained 
by Stewart. Also in this analysis we shall 
assume that the basic current is infinitely 
broad. This is, however, a rather serious re- 
striction since we know that the westerlies 
by no means are broad compared to the north— 
south extension of the mountains. In sec. TV 
we shall try to get some idea of the mountain 
influence on quite narrow currents. The re- 
sults from II and IV will then be used in the 
last section for a qualitative discussion of the 
mountain influence upon the westerlies in 
middle latitudes in the northern hemisphere. 
III. Stationary flow pattern generated by 
a circular mountain in an infinitely 
broad westerly current 


In accordance with the discussion above we 
consider a homogeneous, barotropic and in- 
compressible atmosphere of depth D. A 


- straight infinitely broad westerly current of 


the intensity U (= const) is deformed by a 
mountain, the height of which is h = h (x, y). 
The x-axis is directed to the east and the y-axis 
to the north. For this atmospheric model the 
motion may be described with the aid of the 


vorticity equation: 


d (55) =o, f=2Qsing (1) 


dt\ D 
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¢ is the vertical component of the relative 
vorticity for an air column with the depth D, 
Q is the angular velocity of the earth and » 
is the latitude. The height of the upper bound- 
ary of the atmosphere is assumed to be con- 
stant and equal to Dy. Thus 


D--h = Dy = const! (2) 


The continuity equation then may be written: 


EidD . ou 4 ov 


Te y (3) 


where 4 is the component of velocity in the 
x-direction and v is the component of veloc- 
ity in the y-direction. The equations (1) and 
(3) form a set of differential equations, from 
which u and v are determined as functions 
of x, y and t. In order to solve this system the 
equations must be linearized. If u and v denote 
the perturbation velocities in the horizontal 
plane we thus assume 


Wi ean) 
h< D, (4) 
CE < fo 


These assumptions evidently imply that the 
height of the mountain, h, is small compared 
to Dy. Then the continuity equation is trans- 
formed to 


RE (s) 
DM dt Jie oye Dias À 


and if stationary conditions are considered the 
vorticity equation reduces to 


where 


1 In fact D, must be a function of y in order to get the 
proper pressure field for the basic current. An estimate 
shows, however, that no considerable error is introduced 
by treating D, as a constant. 
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We now introduce a new variable Y, defined by 
' = ne ce h | 
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The continuity equation (5) is identically satis- 
fied by this definition. If the expression for v 
is introduced in (6), this can be integrated 
with respect to x. For x = co we assume that 
G = 0, Y= 0 and ht 0. Then (6) may be 


written 
93V 
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If Y and h are assumed to be independent of y, 
this equation is reduced to the one treated 
by CHARNEY and ELIASSEN. It is possible to 
solve this equation quite general, taking the 
spherical shape of the earth into account. 
The computations become, however, rather 
tedious in this case and we shall restrict 
ourselves to a qualitative discussion, based 
on an exact computation for the very simple 
case of perturbations caused by a circular, 
solitary mountain in an unlimited zonal cur- 
rent. 

The solution of the non-homogeneous differ- 
ential equation (9) is obtained as the sum of 
the general solution to the corresponding ho- 
mogeneous equation and of a particular solu- 
tion to the non-homogeneous one. From (9) 
is seen that the homogeneous equation is 
obtained by putting h =o, ie. these solu- 
tions will describe the free perturbation which 
may exist independently of the orography of 
the surface of the earth. The particular solu- 
tion of (9), on the other hand, describes the 
stationary wave pattern caused by the moun- 
tain. It is possible to find a particular solution 
in a simple manner, if we assume that the per- 
turbation pattern is symmetrical with respect to 
the y-axis. This may be considered as correct 
to a first approximation in the vicinity of the 
mountain, (cf. fig. 4). Because of the assump- 
tion that h =h (r) where r? = x? + y? we 
transform (9) into polar coordinates. If a new 
independent variable o is introduced, we ob- 
tain 
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eV (11) 


A particular symmetrical solution of this 
equation may be written 

= P(e) + Q(o) sind (12) 
thus only considering the zonal harmonics of 
zero and first order, because of the structure of 
the function on the right hand side of equ. (10). 


P(e) and Q(o) are particular solutions of the 
two equations 
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(13) 
Thus P(o) and Q(o) are given by 


P(e) = — TA [a Meder (0) de — 
— Jo) (2) e* Yo(e) de] 
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where J, (0), Jı (0), Yo (e) and Y, (0) are ordi- 
nary Bessel functions. Now u and v can be 
computed from (8) and (12). As, however, 
diva v +0 the velocity field cannot be re- 
presented by a system of streamlines. A new 
function ® is then introduced, defined as the 
streamfunction for the volume transport of 
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Fig. 3. The theoretical flow pattern generated by a cir- 

cular mountain in a homogeneous basic current. Dashed 

lines are troughs and ridges and the dash-dotted line 
indicates maximum zonal velocity. 


the total layer. This is possible as the fluid 
is assumed to be incompressible. Consequently 


br - 


If the same approximations are made as in the 
linearization of th> basic equations, the func- 
tion ® is given by 


D u) 


(16) 


As seen from (15) the velocity vector is tangent 
to the streamlines for the function ®. Thus 
the scalar field ® may be used to illustrate the 
deformation of the current around the moun- 
tain. We will now first discuss the pattern 
thus generated by a simple circular mountain, 
the height of which is given by 


r 


2 (17) 


hehe 


En 
l is chosen so that h = — hy in a distance of 
e 


1,000 km from the centre of the mountain. 
The order of magnitude of the dimensions is 
the same as for the Rocky Mountains. 

A polar coordinate system has been used 
for the derivation of the perturbation pattern 
created by this circular mountain. This causes 
some difficulties for the applicability of the 
theory to real conditions because of the spher- 
ical shape of the earth. Furthermore, ß is nota 
constant but varies with latitude. We assume, 
however, that the radius of the earth is great 
compared with the wave-length of the waves 
and that the variation of 8 with latitude is 


Fig. 4. The observed flow pattern in winter at the 20,000 

foot level (after Normal Hemispherical Mean Maps). 

Streamlines are approximately drawn for every (Az 

sin 9) m instead for every Az m, in order to make 
the pattern comparable with fig. 3. 


small. We may then carry through the com- 
p'itations considering the surface of the earth 
as plain and ß as a constant during the integra- 
tion, but varying ß and taking the smaller 
circumference of the earth in high latitudes 
into consideration when constructing the final 
perturbation pattern. This method cannot be 
expected to give more than a first crude 
approximation of the exact solution. 

Another difficulty results from the fact 
that the basic current U is not homogeneous as 
assumed in the theory. We must then choose 
a reasonable mean value which may be con- 
sidered as representative for the actual atmos- 
phere. It is obvious that the appropriate value 
of U varies with season. According to CHARNEY 
and ErrAssen the proper value for U in our 
equivalent-barotropic atmosphere is about 
equal to the mean velocity at a level some- 
where between 600 and soo mb in the real 
atmosphere. We consider winter conditions 
and choose U = 20 m.p.s., which is approxi- 
mately the mean value at the soo mb level 
for a 10 degrees broad band centered around 
latitude 40° N. 

The flow pattern generated by this circular 
mountain is shown in fig. 3. Some interesting 
features of this figure shall be pointed out as 
well as some conclusions drawn from the 
basic expressions for the perturbations given 
by (14). This will be useful in the last section 
of the paper, when discussing the general 
importance of the orography of the earth. 

a) The pattern is rather sensitive to varia- 
tions of the horizontal scale of the mountain 
as seen from the expressions for P(g) and Q(e). 
This depends upon the fact that the integrand 
contains a factor 93. An estimate shows that if 
the horizontal dimensions of the mountains 


190 


decrease to one half, the values of P(o) and 
Q(o) are reduced to about one tenth. Thus 
only mountains of at least the dimensions of the 
Rocky Mountains are of importance for the gen- 
eration of planetary waves in the westerlies. 

b) The computed pattern may be considered 
to be a combination of two simple pattern, 
one in principle the same as the pattern QUENEY 
has derived as being generated by an infinite 
mountain range because of the flow over the 
mountain, the other similar to the pattern 
which Stewart has computed as an effect of 
the flow around a circular cylindrical obstacle. 
Both components are of almost equal impor- 
tance in the immediate vicinity of the moun- 
tain. The STEWART-pattern is, however, com- 
pletely negligible already one wave-length 
downstream. A rather good approximation of real 
conditions may therefore be expected from a treat- 
ment of the problem in one dimension in the hori- 
zontal plane. Modifications of considerable im- 
portance exist only in the vicinity of the mountain. 

c) The flow pattern varies with the value 
of the zonal velocity U in two respects. In the 
first place the wave-length is a function of U, 
as seen from equations (11) and (14): The 
wave-length is proportional to the squareroot 
of U. HAURwITZ (1940) and CHARNEY and 
ELIASSEN have shown that a correction is 
introduced if the disturbances are assumed to 
have a finite extension in north— south direc- 
tion, i.e. if the basic current is limited by fixed 
boundaries to the north and south. The wave- 
length is in fact very sensitive to the value of 
the width of the current. The real westerlies 
do have a finite meridional extension but it is 
difficult to choose a value which may be con- 
sidered as representative for this special pur- 
pose since no fixed boundaries exist. Thus no 
assumption of boundaries to the north and 
south have been introduced. In the second 
place small changes in the shape of the flow 
pattern in the vicinity of the mountain are 
found for varying U. For certain values of 
U the ridge at the mountain is somewhat 
stronger, for other values it is weaker than for 
the particular value of U used in the construc- 
tion of fig. 3. These variations are small and 
the principal features of the pattern are the 
same within the atmospheric range of U. 

d) The tilt of the trough and ridge lines 
increases with increasing distance from the 
mountain. 
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e) Two regions with intensified zonal flow 
are found theoretically, one just to the north 
of the mountain and the other on the south- 
eastern side of the first trough downstream. 


IV. Influence of a mountain range on a 
very narrow current 


It should be kept in mind that the theory 
advanced in the preceding section is based 
on the assumption that the width of the basic 
current is large in comparison with the north 
—south extent of the mountain. In the case 
of a very narrow current the mountain in- 
fluence may have an entirely different char- 
acter. When the current passes a mountain, 
horizontal divergence takes place resulting in a 
generation of anticyclonic vorticity within the 
current. If the mountain is large enough and 
the current not too strong, the anticyclonic 
curvature of the current will be so strong 
that the current is turned back to its source 
region. To demonstrate this particular oro- 
graphic effect we shall disregard the latitude 
variation of the Coriolis parameter and con- 
sider the extreme case of an infinitely narrow 
current striking a mountain range at right 
angles, though this last assumption is not 
essential. Furthermore we disregard the in- 
fluence of the environment completely. As 
before we assume that the vertical contraction 
of the air column is made up for by an equiva- 
lent transversal streching (= the velocity, c, 
of the current remains constant). As the rela- 
tive vorticity of the undisturbed current is 
zero and if fy denotes the constant value of 
the Coriolis parameter and D, the undis- 
turbed height of the air column, (1) may be 
written ; 


an 

a == ia (18) 
or 

Co es h (19) 
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Disregarding changes in the horizontal shear 
we get 


(20) 


where R is the radius of curvature of the cur- 
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Fig. 5. Deformation of a narrow current because of 

divergence when striking a mountain range. The upper 

part of the figure shows the horizontal flow pattern 

and the lower part a vertical cross-section across the 
barrier. 


rent. Introducing the analytical expression for 
R (20) becomes 


dy 
dx? mer 
HET = 
I 
dx 
which may be transformed to 
dy 
dx f 
‘(epee 
ig 
dx 
or 
d (sin x) = fo h dx (23) 


san 
a being the angle between the deflected jet 


and the x-axis. As the undisturbed current is 
assumed to be straight westerly we get 


ee 
Sina Br h dx (24) 


where x, is the x-coordinate for the point 
where the current strikes the mountain (cf. 
fig. 5). If the current is to cross the mountain. 
|x| must be less than 90° everywhere. Hence 


fo ax 
ope fiat (25) 
We may define a critical velocity 

re Pali ‘| h dx (26) 


where Ax is the width of the mountain. If 
CCE (27) 


the current is turned back to its source region. 
Puttnosj, = 10, sec RD — Tos ein. x 
1000 km and assuming a mean height of the 
mountain range of 2000 m we get ty, = 
20 m sec.t. 

If the variation of the Coriolis parameter 
with latitude is taken into account c,, will 
be modified. For a westerly current the crit- 
ical value will be smaller, for an easterly cur- 
rent greater than given in (26). This is easily 
seen from a qualitative discussion. It is be- 
lieved that the effect indicated above may be 
of importance in the atmosphere even if 
(26) merely gives the order of magnitude of 
the critical value. As shown by the numerical] 
example only mountains of the same order of 
magnitude as the Rocky Mountains are effect- 
ive. 


V. The importance of the orography of the 
earth for the general character of the 
westerlies 


The discussion in this section must necessa- 
rily be rather vague since a very simplified 
atmosphere has been used for the theoretical 
considerations on which it is based. Following 
the general procedure outlined in II we shall 
at first discuss conditions at upper levels (500 
mb), then consider the effect at the surface of 
the earth and in doing so also take account 
of the thermal influence from below. The 
degree of similarity between the theoretical 
results and the observed pattern will show to 
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what extent the hypothesis of a dynamic- 
orographic control of the upper level wester- 
lies may be considered as correct. 

As seen from III, a as well as from IV we 
need to consider only those mountains which 
are of the same order of magnitude as the 
Rocky Mountains. Hence it follows that the 
only mountains which may be of real importance 
in generating planetary waves in the westerlies 
are the Rocky Mountains and the mountains in 
the interior of Asia. 

Let us at first consider the general character 
of the waves around the northern hemisphere 
with regard to wave-length and the location 
of the troughs and ridges. It was pointed 
out in III that a rather good first approxima- 
tion may be obtained from a treatment of the 
problem in one dimension in the horizontal 


plane as made by Carney and EL1Assen. It is, 


essential, however, that we consider the belt of 
maximum westerlies and that the mountains 
are located in the centre of this current. It then 
seems reasonable to assume that the relatively 
weak currents in the environment ajust them- 
selves to the wave pattern in the belt of maxi- 
mum kinetic energy. Therefore CHARNEY’s 
and ELIASSEN’S computations are not quite 
representative since they considered the lati- 
tude belt between 40° N and 50° N, which is 
to the north of the highest parts of the Asiatic 
mountains. It is nevertheless of some interest 
to compare their results with actual condi- 
tions. Fig. 6 shows the soo-mb profile around 
the whole hemisphere computed by CHARNEY 
and ELIASSEN and the observed one obtained 
from the hemispherical upper air mean charts 
for winter, which have been published by 
SCHERHAG. The agreement between the com- 
puted and observed patterns is good, in fact 
“better than CHARNEY and ELIASSEN found them- 
selves using Normal Weather Maps published 
by U.S. WEATHER BUREAU (1944). 

A comparison between winter and summer 
conditions is of special interest when dis- 
cussing the relative importance of dynamic- 
topographic versus thermal factors for the 
control of stationary waves in the westerlies. 
The thermal influence from the surface of 
the earth is quite different in the two seasons 
and should thus be expected to bring about a 
rather far-reaching rearrangement of the wave 
pattern, as was already pointed out in the in- 
troduction. From fig. 1 is seen that such a 
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Fig. 6. Computed stationary profile around the hemis- 

phere of the soo mb surface (dashed curve) according to 

CHARNEY and ELIASSEN (1949) and observed mean con- 

ditions (solid curve) after SCHERHAG (1948). The curves 

represent winter conditions in the belt between 40° N 
and 50° N. 


change does not occur. On the other hand if 
one assumes that the large mountains are of 
importance in generating these waves, the 
observed similarity between summer and win- 
ter patterns is obtained as a necessary conse- 
quence. The stationary wave-length is propor- 


tional to YU/ß. The zonal flow decreases 
from winter to summer, but at the same time 
the belt of maximum westerlies is displaced 
farther to the north. These two facts counter- 
act each other and the stationary wave-length 
(measured in degrees longitude) remains almost 
constant. If we assume that the zonal wind 
decreases at the rate of so per cent and that 
the meridional displacement is 10° we obtain 


lets 

Feige (28) 
where L» and L, denote the stationary wave- 
length in winter and summer respectively. 
Since furthermore the waves are anchored 
by the mountains, there should be little or no 
difference in the location of the waves during 
the two seasons. This points to the correctness 
of interpreting the waves as generated mainly 
by orographic influences. As seen from (28) 
there is a slight decrease of the wave-length 
from winter to summer which in reality 
shows up as a tendency for a fourth trough to 
develop in summer off the American westcoast. 
This is in agreement with the fact that two 
polar fronts are often observed in the Pacific 
in summer. These two troughs can hardly be 
explained by assuming that the thermal in- 
fluence from below is of conclusive impor- 
tance. 
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Fig. 7. Theoretically computed profile of the zonal 
velocity in a cross-section 3,000 km downstream from the 
circular mountain. 


It is impossible to pursue the discussion 
further without considering the meridional 
variation of the height of the mountains and 
thus allow for flow around as well as over 
them. There is an essential difference between 
the mountains in America and in the interior 
of Asia in this respect. 

The Rocky Mountains extend all along the 
American westcoast and thus the current 
must always pass over the range in spite of the 
fact that the westerlies are displaced northward 
from winter to summer. Therefore the actual 
wave pattern should in principle be in accord- 
ance with the pattern obtained from the 
treatment of the problem in one dimension in 
the horizontal plane also in the vicinity of the 
mountains. It must be pointed out, however, 
that the Rocky Mountains are broadest and 
highest at about 40°N. According to the 
previous discussion this should result in a 
tendency for the current to flow around these 
parts of the mountain range. Some indications 
of this are also found by a comparison between 
fig. 3 and 4. The general structure of the mean 
flow pattern in the vicinity of the Rocky Moun- 
tains has several similarities with the pattern 


son 50°N 40°N 30°N 20°N 
Fig. 8. Zonal wind profile across the westerlies at the 
20,000 foot level. The solid curve is computed from Nor- 
mal Weather Maps along longitude 85° W in January 
and the dashed curve is obtained from the mean cross- 
section in winter along 80° W, which has been published 


by Hess (1948). 
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Fig. 9. Zonal wind profile in winter at 76° E across the 
westerlies at the 20,000 foot level (after CHAUDHURY, 1950). 


which we have computed for the flow around 
a circular mountain. Some indications of the 
tendency for the splitting of the basic current 
may also be found by comparing the com- 
puted and observed windprofiles across the 
current some thousand kilometers downstream 
from the mountains (fig. 7 and 8). However 
the dominating deformation caused by the 
Rocky Mountains results from the fact that 
the westerlies are forced to flow over these 
ranges. 

The mountains in the interior of Asia, on the 
other hand, do not have the same extent in the 
north-south direction. Hence the westerlies en- 
close the mountains in winter, while the mean 
jet is found well to the north of the Tibetean 
Plateau in summer. In winter the westerlies 
are split in two branches, one on each side of 
the mountains (fig. 9), which is in qualitative 
agreement with the computed effect (fig. 7). 
Further downstream a well developed wave 
pattern is found. In summer, on the other 
hand, the jet is affected only by the relatively 
low mountains further to the north. This may 
explain the fact that there exists only one 
mean jet at that time of the year and that the 
waves in the westerlies further downstream 
have small amplitudes (see fig. 1). The moun- 
tains obviously are of considerable importance, 
but it must be kept in mind that the Hima- 
layas in winter acts as a boundary between 
cold air to the north and warm air to the south. 
This certainly also is of great importance for 
the character of the westerlies in the vicinity 
of the mountains.! 


1 After this manuscript had been completed the writer’ 
attention has been called to a paper by YEH (1950), 
which is also published in this issue of Tellus. Yen’s data 
support most of the conclusions drawn in the last two 
paragraphs. Some similar results have also been obtained 
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The division of the westerlies in two bran- 
ches has one consequence which is essential 
when discussing intensity variations of the jet. 
SUTCLIFFE (1940) and Namias and CLAPP (1949) 
have pointed out the significance of the con- 
fluence mechanism for the intensification of a 
jet. This mechanism means simply the flowing 
together of two different air currents, one 
stream coming from the warmer southerly 
latitudes, the other from the colder north. As 
seen from fig. 3 a confluence pattern is gener- 
ated on the lee side of a mountain in a westerly 
current. If these two streams joining to the 
east of the mountain are of different thermal 
character, there should be an intensification 
of the current downstream. It is a fact that in 
the mean the most intense jets are found in the 
regions to the east of the large mountain 
areas (cf. Namıas and Crapp). The difference 
between summer and winter conditions may 
be explained by the seasonal differences in 
the confluence mechanism. These variations 
should thus depend a) upon the extent to 
which the mountains are capable of splitting 
the westerlies, b) upon the magnitude of the 
thermal contrast between the two currents. 

Until now we have only discussed the nor- 
mal upper air charts as representing a mean 
condition at upper levels. In reality consider- 
able deviations from this normal pattern may 
persist for periods as long as several weeks. 
It must be kept in mind, however, that the 
latitude and speed of the maximum wester- 
lies which have been used above are the nor- 
mal seasonal values; deviations from these nor- 
mal values must cause changes in the stationary 
patterns generated by the mountains. Further- 
more the effect discussed in IV may be of 
importance during low index periods. This 
may for example be one explanation of the 
stationary extended troughs, that sometimes 
are found over the Rocky Mountains. 

It must be emphasized that the solutions 
used in the preceding discussion are not 
unique; but they represent the particular flow 
patterns which are the immediate consequence 
of the mountain barriers. Besides these waves 
there may exist stationary waves which are 
not fixed by the orography, and which 


in two theses by Borrr and CHAUDHURY, which were 
carried through under the supervision of Dr H. Riehl 
at the University of Chicago. 
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disappear if the mean pattern 1s computed 
over a sufficiently long time. Such waves may 
influence the character of the mean pattern 
during periods of a few weeks duration. 

If one accepts the hypothesis that the moun- 
tains are of primary importance for the genera- 
tion of the mean stationary long waves in the 
westerlies, the ridge at the European westcoast 
and the trough over the White Sea must be 
considered as to a certain extent being formed 
in resonance with the waves at the Rocky Moun- 
tains (cf. fig. 1). On the whole, the positions of 
the ridge and trough are in agreement with 
the computations made in III if one allows 
for the correction that must be introduced 
because of the finite width of the current. 
In addition, it must be pointed out that the 
mean ridge to the west of Europe to some 
extent may be the result of the fact that this 
region appears to be a preferable place for 
blocking to start (ELLIOTT and SMITH, 1949; 
Rex, 1950). This effect has not been considered 
in our theoretical treatment. It is probable 
that blocking may be of some importance for 
the formation of the upper ridge but it should 
be noted that the ridge often is found even in 
the absence of blocking. It thus seems quite 
probable that the mean ridge to a considerable 
extent is caused by the topography of the earth 
further to the west. 

Thus we have found that there are several 
features of the upper flow pattern which indi- 
cate that the orography of the earth seems 
to be of considerable importance for the mean 
large scale wave pattern around the northern 
hemisphere. We shall now finally discuss the 
interrelation between this upper flow pattern 
and conditions at the surface of the earth. 
Following the general chain of reasoning 
outlined in the beginning of the paper we 
shall try to get an idea of the adjustment in the 
layer below soo mb considering the restric- 
tions which are impressed by the flow pattern 
aloft, the surface friction and the differential 
heating from below. Assuming then that the 
flow pattern above a certain height H (= 5 km) 
to a great extent is controlled dynamically it 
follows that the pressure distribution at this 
level is known apart from a constant. Adjust- 
ment below means a redistribution of mass, 
so that the pressure distribution at the ground 
fits the requirements of the monsoonal circu- 
lation in the lower part of the troposphere. 


ON THE INFLUENCE OF THE EARTH’S OROGRAPHY ON THE WESTERLIES 


This redistribution of mass in turn means a 
change of the mean temperature of the air 
column below H and consequently the mean 
temperature of this air column is at least in part 
controlled dynamically from above, in addition to 
the direct thermodynamic control from below. This 
conclusion, if supported by further investiga- 
tions, is of great importance to dynamic cli- 
matology (BERGERON, 1930). It implies that 
both pressure and temperature at the ground 
are in part controlled by the mean flow pattern 
aloft, which in turn are intimately related to 
the orography of the earth all around the 
hemisphere. At the present time it is difficult 
to make a qualitative estimate of this influence 
of mountain systems upon the climate of 
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far distant regions, but qualitatively we 
are justified in assuming that both the upper 
and surface high pressure ridge from the 
Azores to France in winter and the north- 
erly position of the mean storm tracks over 
the eastern Atlantic resulting in relatively 
mild winters in northwestern Europe are not 
exclusively to be looked upon as an effect of 
the warm ocean but also as a result of the in- 
fluence of the Rocky Mountains. 
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I. An Aerological Study of Blocking Action 
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Abstract 


A comparative study of five cases of blocking action in the upper westerlies 1s presented with 
the two-fold objective of obtaining a more complete description of this phenomenon and of 
attempting an explanation of its initiation, development and dissipation in terms of parameters 
derived from an hydraulic analogue. The effect of European blocking upon the zonal and merid- 
ional circulation indices and the röle of blocking in the evolution of hemispheric cir- 
culation patterns are discussed. In Part II, to be published later, the results of a preliminary 
semi-statistical analysis of blocking action, together with a description of the regional pre- 
cipitation and mean surface-temperature anomalies produced by blocking, are presented. The 
statistical results obtained are discussed in relation to recent climatic fluktuations. 


Introduction 


Within the sector lying between 120° W and 
60° E longitude and at the soo mb level, the 
normal circulation pattern may be described 
as a relatively broad westerly current (jet) 
circulating around the circumpolar vortex 
and bounded to the south by the sub-tropical 
high-pressure belt. Pronounced variations in 
the width and strength of these westerlies 
and in the extent and intensity of the cir- 
cumpolar vortex are observed from day to 
day and month to month. A particular longi- 
tudinal variation in this circulation pattern 
frequently appears on the synoptic chart ir- 
respective of season. This anomaly expresses 
itself as follows; the flow over the North 
American continent and western Atlantic, 
which is an essentially zonal motion at a 
relatively high energy level, subsequently 
“breaks down” over the European continent 
or eastern Atlantic into a more cellular mo- 
tion at a lower energy level. A series of wave 


disturbances are most generally observed super- 
imposed on the zonal motion, whereas the 
more meridional type is characterized by the 
appearance of a train of cyclonic and anti- 
cyclonic vortices; the energy released alon 
such a flow channel is absorbed in the large- 
scale horizontal vortices just mentioned. Im- 
portant characteristics of this pattern are that 
it is persistent and quasi-stationary and that 
the change in each case from an essentially 
zonal to an essentially cellular flow occurs 
abruptly at a rather well defined longitude, 
usually between 30°W and 10° E, and not as a 
progressive change from west to east. The 
trough-and-ridge diagram given by Hov- 
MOLLER (1949) for Nov. 1945 (a blocked period) 
shows these features quite clearly. Figure 1 
illustrates a well-developed, although not unu- 
sual, example at the soomb level of the circu- 
lation anomaly just described. It is interesting 
to note the similarity between this pattern 
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Fig. 1. An example of well developed blocking action over 
Eurasia. Absolute topography in dynamic decameters. 


and the “Viererdruckfeld” described by RAET- 
HJEN (1949) in connection with his treatment 
of cyclogenesis; the “Delta” used in his ter- 
minology corresponds to the point of jet 
division shown in Figure 1. 

A usual feature of the macro-turbulent 
region in well-developed cases is the forma- 
tion of a quasi-stationary warm ridge or anti- 
cyclone just downstream from the point of 
“breakdown”. This anticyclone may sub- 
sequently progress toward the east or move 
slowly retrograde, or it may remain essentially 
stationary in company with its associated high- 
pressure cell at the surface. Situations, showing 
the pronounced development of such a quasi- 
stationary anticyclone were first noted by 
GARRIOTT (1904), and are now commonly 
referred to as examples of “blocking action” 
(the surface anticyclone being called the 
“blocking high”). Considerable attention has 
been directed to them because of their pro- 
nounced influence upon other weather features 
even at a considerable distance up- and down- 
stream. 

In particular the work of Namras (1947) 
and ELLIOTT and SMITH (1949) should be noted. 
These writers have presented a comprehensive 
description of the phenomenon and its de- 
velopment in the lower troposphere and have 
discussed its effect upon synoptic developments 
at the surface. YEH (1949) has presented a 
theory of blocking which explains the more 
striking features of the process in terms of the 
retrograde motion of the longer wave com- 
ponents Ina) nondivergent baratropic at- 
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mosphere. In a more recent paper BERGGREN, 
Bouin and Rosssy (1949) have described the 
blocking process at the soomb level and its 
association with developments in the lower 
atmosphere. They point out the difficulty of 
explaining the blocking development in terms 
of the growth of wave-like perturbations and 
suggest that the standing-wave or hydraulic- 
jump, which may develop in open-channel 
flow under certain circumstances, resembles 
at least superficially the blocking phenomenon 
in the atmosphere. RossBy (1950) has further 
developed this analogy and, using a simplified 
model, has been able to derive an analytical 
expression for one necessary condition for 
locking. 

The objective of this paper (Part I) is to 
extend the aerological description of blocking 
action in the middle troposphere obtained by 
BERGGREN, BOLIN and RossBy (1949) and to 
attempt an explanation of its initiation, de- 
velopment and eventual dissipation in terms 
of parameters derived from the hydraulic- 
jump analogue. In a second paper (Part II), 
to be published later, the results of a pre- 
liminary semi-statistical analysis of blocking 
action (1933—1949), together with a descrip- 
tion of the regional precipitation and mean 
surface-temperature anomalies produced by 
blocking, are presented. These statistical re- 
sults are discussed in relation to recent 
climatic fluctuations. 

As a result of work to be reported in Part II, 
a definition of blocking action at the soomb 
level has been derived which will be used in 
the following discussion. A blocking case 
must exhibit the following characteristics: 


a) the basic westerly current must split into 
two branches, 

b) each branch current must transport an 

appreciable mass, 

) the double-jet system must extend over at 

least 45° of longitude, 

d) a sharp transition from zonal type flow 
upstream to meridional type downstream 
must be observed across the current split, 
and 

e) the pattern must persist with recognizable 
continuity for at least ten days. 


Q 


In a situation satisfying these requirements 
blocking is said to be initiated when (a) oc- 
curs and is said to have dissipated whenever 


CONTOUR CHANNEL AT THE 500 M8. LEVEL 
nn 46 JUNE 1949 


28 JUNE 1969 


Fig. 2. A series of plates showing the development of a blocking situation at the soomb level. The flow channel 

occupied by the hemispheric westerlies at 0300 GCT on 12, 16, 20, 24 and 28 June and on 2 July (1949) is 

indicated by the heavy black contour lines. The temperature field is shown in colour: dark blue< —25° C< light 

blue < —20° C< clear < —10° C< light red < —5° C < dark red. Contour heights in dynamic decameters. 
Temperatures in degrees Centigrade. 


(a), (b), (c) or (d) are no longer satisfied. Al- 
though such a definition is subjective in char- 
acter it has proven sufficiently definitive to 
permit the compilation of a catalog of block- 
ing cases. 

Five examples of blocking action, as defined 
above, have been analyzed in detail; four cases 
occurred over the eastern Atlantic and Europe 
(Jan. 1947, Feb. 1948, June—July 1949 and 
Jan. 1950) and one case occurred over the 
eastern Pacific and western N. America (Jan. 
1949). These examples, both seasonally and 


geographically scattered, were selected for 
detailed study from the 113 cases cataloged 
in Part II and are thought to be representative 
of typical well-developed blocking situations. 
In so far as the aerological features of blocking 
action are concerned, the results of the above 
mentioned studies are in remarkable agreement. 
Accordingly the presentation of one case 
description will suffice. 

The basic data used in these studies consists 
of the daily surface and soomb hemispheric 
analyses published in Daily Series, Synoptic 
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Weather Maps (U. S. Weather Bureau) for 
the January 1949 case and in Täglicher Wetter- 
bericht (Zentralamt Bad Kissingen) for the 
June—July 1949 and the January 1950 cases. 
Daily soomb charts prepared by Mr. B. Bouin, 
University of Stockholm and surface charts 
prepared at the Swedish Hydrological and 
Meteorological Institute were used for the 
February 1948 case. In order to obtain a 
continuous daily analysis over a sufficient 
longitudinal span (150° W—o°—90° E) for 
the January 1947 case, North American 
soomb and surface analyses published in 
micro-film form by the U. S. Weather Bureau 
were joined with similar European analyses 
published in Täglicher Wetterbericht (Zen- 
tralamt Bad Kissingen). In order to minimize 
subjective errors no individual analysis has 
been performed by the writer except that 
necessary to construct the composite 1947 
charts referred to above. 


The June—July 1949 Blocking Case 


The example selected for detailed descrip- 
tion developed over the eastern Atlantic and 
Europe during June 1949. As defined previous- 
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ly, the period in question began on 16 June 
with the active development of a high-pressure 
ridge located approximately at 10°W longi- 
tude and ended 22 days later on 7 July with 
the disappearance of this anticyclone over 
Scandinavia. The soomb chart for 24 June, 
illustrating this blocked case at near maximum 
development, is shown in figure 1. 

The set of colour plates shown in figure 2 
was prepared in order to more clearly demon- 
strate the sequence of events occurring at the 
soomb level in connection with the develop- 
ment of a “block”. In each instance the two 
height contours which most nearly bound 
the circumpolar westerly stream at soomb 
-vere selected and are represented by the heavy 
solid lines in the figure. These lines approxi- 
mately delineate the flow channel occupied 
by the upper westerlies, and are shown super- 
imposed on the observed temperature field at 
that level. It is clear that the isobaric (con- 
tour) pattern can give only an approximate 
picture of the actual circulation field; in view 
of the stability of the pattern in blocked cases 
however, we may expect the approximation 
to be a rather good one. The projection used 
in the colour plates is a transformed conic, 


zone for the June—July 1949 case. North—south distances measured in ° latitude from mid-point of the jet. 


Mid-points in upstream sector: 16 June — 50°N, 24 June — 43°N and 2 July — 52° N; in 


blocked zone: 


16 June — 52° N, 24 June — 57° N and 2 July — 58° N. 
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tangent at so’N — this projection is ortho- 
gonal but expands areas lying north of so’N 
latitude and shrinks those lying to the south. 

Prior to the development of the block (12 
June) we observe a regular pattern of shallow 
waves around the entire half-hemisphere, 
though they are somewhat confused over 
Europe due to the presence of a stagnant 
cold cut-off cyclone. On 16 June the fırst 
sign of a developing block appears in the 
form of a split in the westerly jet at about 
20°W longitude. Meanwhile the wave pattern 
persists with somewhat increased amplitude. 
By 20 June a well developed block is apparent 
which may be seen in each of the succeeding 
plates (24 June, 28 June and 2 July). In these 
last three illustrations the essentially zonal 
character of the westerlies over North Ame- 
rica as contrasted with the markedly cellular 
(macro-turbulent) character of the flow over 
Europe, as well as the almost stationary posi- 
tion of the blocking ridge, may easily be seen. 
The characteristic sharpness of the “flow dis- 
continuity’ across the block is also clearly 
shown. Active large-scale horizontal mixing is 
evident from an inspection of the temperature 
field in the blocked zone, and a consequent 
relaxation of the mean meridional tempera- 
ture gradient results. 

In order to show changes occurring during 
the development of the block in the tempera- 
ture and velocity fields at the soomb level, 
figures 3 and 4 have been prepared. The 
longitudinal sector between 100°W and 40°W 
is taken as representative of upstream condi- 
tions, and the sector between 10°W and so°E 
is taken as representative of the blocked zone. 
Figure 3 shows a series of mean zonal (space- 
averaged) geostrophic velocity profiles com- 
puted over the sectors just defined on 16 and 
24 June and on 2 July; figure 4 shows the 
corresponding temperature profile series. The 
following significant changes are observed 
within the upstream zone at the soomb level 
during the blocked period: 


1. The velocity profile sharpens and the velo- 
city maximum increases slightly during the 
first half of the period; later the jet broadens 
again. 

2. The latitude of the mean jet in this zone 
shifts slowly southward during the develop- 
ing stages (from so’N to 43° N). As the 
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block deteriorates the mean jet moves 
rather rapidly northward again (to 52° N). 
3. A small but relatively uniform warming-up 
of the air in this zone is observed over a 
wide latitudinal span (from 20° N to 85° N). 
This observation suggests a general sub- 
sidence effective over the whole upstream 
region and associated with the southward 
displacement of the jet during the de- 
veloping stages of the block. Later the tem- 
perature falls again in conjunction with the 
final northward movement of the jet. 


The following significant changes are ob- 
served within the blocked zone: 


1. The original single westerly jet, which was 
located at about 52° N latitude, is replaced 
during the development of the block by 
two westerly jet-branches, which are located 
approximately at 78° N and 36° N latitude 
with a zone of easterly mean flow between 
(at about 58° N). During the final stages of 
blocking the jet branches weaken and 
move further apart. It is perhaps signi- 
ficant to point out that these two jet 
branches are of approximately equal inten- 
sity throughout the active stages of the 
development — the velocity profile obser- 
ved is suggestive of the action of a momen- 
tum jump such as that described by Rossgy 
(1950). 

. The mean latitude of the blocked zone 
moves slowly northward throughout the 
period (from 52°N to 58°N). 
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Fig. 4. Space-averaged temperature profiles computed 
over the upstream sector and over the blocked zone for 
the June—July 1949 case. Distances measured as in fig. 3. 
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3. A pronounced warming (6° C at 64° N) in 
that part of the blocked zone lying to the 
north of 52°N latitude and a cooling 
(4°C at 45° N) in that part of the zone 
lying to the south of 52°N latitude is 
observed. The mean meridional tempera- 
ture gradient across the zone has there- 
fore been greatly weakened, or even re- 
versed, as a result of the large-scale hori- 
zontal mixing processes operating. 


In concluding this description of the deve- 
lopment of a typical block the changes pro- 
duced in the velocity field will be discussed. 
In figures 5 and 6 the development of zonal 
velocity as a function of time is depicted 
graphically. Values of the geostrophic zonal 
speed were computed for successive latitudinal 
bands, five degrees in width and extending 
from 30°N to 80°N latitude, at each tenth 
meridian. Mean values of these speeds, taken 
over the upstream zone (120° W—30°W) and 
averaged over overlapping five day periods, 
are shown as time functions at various lati- 
tudes in figure 5. A similar presentation com- 
puted over the blocked zone (30° W—30° E) 
is shown in figure 6. In both figures isolines of 
geostrophic zonal speed have been drawn 
for each 2 mps. The development of the double 
jet system and its eventual dissipation are 
clearly shown in figure 6. The changes in 
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mean jet width, strength and location which 
were discussed in connection with the velo- 
city profiles may again be seen in these figures. 

At this point it may be of some interest to 
discuss the effect of blocking action upon the 
zonal and meridional circulation indices. This 
is particularly true in view of the numerous 
conflicting statements which may be seen in 
current literature relative to index behavior 
during a blocking development. Both the 
zonal and the meridional geostrophic velocity 
components at the soomb level have been 
computed on a 10° grid for each day during 
the period 16 June—7 July. In table I are 
given the mean geostrophic zonal index va- 
lues obtained, in each ten-degree band of 
latitude from 30°N to 80° N, for the entire 
half-hemisphere (120° W—6o° E), for the up- 
stream zone (120°W—60° W) and for the 
blocked zone (30° W—30° E). Also given are 
normal index values averaged over the same 
zones and the departures of the observed in- 
dices from normal. For a discussion of the 
method of computing circulation indices, and 
their interpretation, the reader is referred to 
STONE and WILLET (1943). 

The normal values were computed from a 
new summer mean soomb chart (June, July 
and August) by FLOHN (1949). With regard to 
the half-hemisphere values it is significant to 
note from table I that only in the region from 
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Fig. 5. Time trend of mean zonal speed computed over the upstream sector for the June—July 1949 case. Westerly 
flow stronger than 10 mps indicated by coarse hatching. Easterly flow indicated by fine hatching. 
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Table I. Zonal Index — Mean Zonal Component at 500mb Level — 16 June—7 July 1949 
ihe) Seas ee ee ee ne ee eee 


120° W—60° E 


Latitude bes 
Zone Departure 
Index Normal from Index 
Normal 

70—80° N 4.8 1.5 33 3.0 
60—70° N 2.5 4.2 — 1.7 4-4 
50—60° N 6.3 8.2 — 1.9 10.6 
40—so- N 107, 1022 0.5 14.2 
30—40° N 8.0 7.6 0.4 4.0 
30—80° N | 07 | 6.5 | 0.2 | 


120° W—60° W 30° W—30° E 
Departure Departure 
Normal from Index Normal from 
Normal Normal 

0.3 2.7 6.0 2.0 4.0 

4.4 o 2.9 3.9 = 1) 

10.0 0.6 3.9 7.0 — 3.1 

13.2 1.0 7:9 7053 on 

7-3 3333 9.0 7-3 1.7 

7.0 | 0.2 | 5.9 | 6.1 | — 0.2 


(All values in meters per second — west wind positive.) 


50° N to 70° N latitude is the blocked period 
characterized by negative departures from the 
normal index values. Particularly at higher 
and to a lesser extent at lower latitudes, zonal 
velocities have greater than normal values. 
Thus, the zonal westerly index as usually 
defined (35° N to 55° N), although it includes 
the westerly maximum, maintains an essen- 
tially normal value during the blocking deve- 
lopment. We may explain the index anoma- 
lies obtained in the upstream zone in terms of a 
contracted and somewhat intensified jet during 
this period. The above-normal index obtained 
at high latitudes (70° N—80° N) in this zone 


MEAN ZONAL SPEED (30 W— 30E) 


no doubt results from the marked northward 
displacement of the jet stream over western 
N. America during the developing stages of 
the block. In the blocked zone we obtain a 
remarkably symmetrical anomaly profile with 
the maximum reduction in zonal How occurr- 
ing at about 50° N latitude. The net changes 
in index produced in each zone over the entire 
latitudinal belt (30° N—80° N) are insigni- 
ficant (+ 0.2 mps). We conclude therefore that 
mean zonal index anomalies, computed for 
blocked periods, are not only dependent upon 
the latitudinal boundaries chosen but also 
depend directly upon the longitudinal sector 
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Fig. 6. Time trend of mean zonal speed computed over the blocked zone for the June—July 1949 case. See figure 
5 for further explanation. 


BLOCKING ACTION IN THE MIDDLE TROPOSPHERE 


over which the computations are extended. 
This result is to be expected in view of 
the sharp quasi-stationary flow discontinuity 
shown in figure 2 at about 20° W longitude. 

In table II below are given the mean me- 
ridional index values obtained at every tenth 
parallel of latitude from 30° N to 70° N in the 
three sectors, 60° W—30° W, 30° W—o° and 
o°—30° E. Also given are the normal me- 
ridional index values and the departures of 
the observed indices from normal. 

Normal values were again computed from 
Flohn’s summer mean soomb chart. In the 
zone from 60° W— 30° W longitude we ob- 
tain small departures from normal except at 
70° N latitude. The strong positive anomaly 
observed at this latitude results from the strong 
southerly flow around the deep cyclonic vortex 
which stagnated over southern Greenland 
during the first part of the blocked period. 
In the zone from 30°W—0° longitude the jet 
split is clearly located between 40°N and 
50° N latitude with northward moving flow 
over higher latitudes and southward moving 
flow over lower latitudes. If we compute 
the mean magnitude of the meridional index 
anomalies over the entire latitudinal band we 


obtain a quantity (|A|) which may be used as a 
measure of “cellular growth” in the circula- 
tion pattern. These values, given in table II, 
show a marked increase in cellular growth in 
the two zones lying downstream from the 
block (30°W—o° and 0°—30° E). 


An Hydraulic Analogue 


Before proceeding with a discussion of the 
initiation, development and dissipation of 


203 


blocking, it may be desirable to review Ross- 
by’s reasoning (RossBy 1950) in a summary 
way and to state here again the particular 
quantitative results which will be used in 
this paper. Let us suppose that we have a 
simple 2-dimensional steady current of width 
2a and speed u in the positive x-direction, 


‘which is independent of y and vanishes in the 


surrounding environment, as depicted in fi- 
gure 7A. Further, if we take the x-direction 
as easterly and assume a constant Coriolis 
parameter, f, which is independent of latitude 
we note that the pressure field necessary to 
balance the motion will have a linear cross- 
current profile. In figure 7B this pressure pro- 
file is shown at the two normal cross-sections 
o and 1 as the solid lines MAN and LAO 
respectively. In this case it becomes immediat- 
ely evident that no pressure gradient exists 
in the direction of flow. If however we intro- 
duce a variable Coriolis parameter, f = fo + By, 
we note that the balancing cross-current pres- 
sure profile is no longer linear but rather is 
of the form MBN and LCO as shown by the 
dashed lines in figure 7 B. Thus we see that this 
second pressure field exceeds the first in mag- 
nitude at every point and further that this 
pressure excess, p’, reaches its maximum value 
at the center of the current. Rossby gives an 
expression for the value of p’ max. as follows: 


em (1) 


/ 
P max. 


where 0 is the density and V the volume 
transport, which can be expressed as follows: 


V=2au (2) 


Table II. Meridional Index — Mean Meridional Component at 500mb Level — 16 June—7 July 1949 
60° W— 30° W 30° W—o° 02 3008 
Latitude Departure Departure Departure 
Index Normal from Index Normal from Index Normal from 
Normal Normal Normal 
TON 9.0 1.8 72 7,2 4.1 ART == u 23 — 4.1 
60° N DAT 2,7. [6) 13.4 3.6 9.8 — 10.7 7 — 12.4 
50° N 1.2 Te o 5.9 DS 4.4 101 0.8 — 10.2 
40° N O 0.8 — 0.8 — 2.4 — 0.4 — 2.0 — 3.5 | — 1.2 — 2.3 
30° N —— 153 — 0,4 20.9 o 0.4 — 04 0 2.7 2.7 
ar EN 1 rn SEN EAS ee sd BE I OR ARE ER SEN PES 
HEN ee Ses en 


(All values in meters per second — south wind positive.) 
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The magnitude of the pressure-excess maxi- 
mum, p’max., increases with increasing 
volume transport or current width. Returning 
then to the simple case of figures 7 A and 7 B 
we see that a particle travelling along the 
streamline located at the center of the current 
will be acted upon by a negative pressure 
gradient and will thereby be decelerated in 
moving from cross-section 0 to 1. This follows 
from equation (1) and the fact that the volume 
transport must remain constant from cross- 
section to cross-section along the current. It 
may be shown in this case that all particles 
moving in the current, except at the bounda- 
ries, will be decelerated by an amount di- 
minishing toward zero as their distance from 
the current’s center approaches | a |. Rossby 
has shown further that the effect described 
above will exist in general, i.e. irrespective of 
the form of the current’s boundary or of 
its velocity profile, as long as symmetry about 
the center latitude prevails. He has therefore 
described an internal physical mechanism 
which is capable of reducing the speed of 
westerly current systems in the atmosphere 
and which suggests that under certain condi- 
tions a single current may tend to split into 
two branches. 

Returning again to figure 7A it may be 
shown that the expression for the momentum 
transfer, M.O.T., across a plane perpendicular 
to the current is given as follows: 


+a py? 
M.O.T. = @ fu (« + er) dy (3) 


which in this case may be easily integrated and 
with the aid of equation (2) expressed in the 
following form: 


MOT=PBV C on m) (4) 


24 u? 


It is evident therefore that for a given volume 
transport the momentum transfer will have a 
minimum value when the current speed is 
equal to: 


where this critical speed is denoted as too. 
Whenever the current speed exceeds the 
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Fig. 7. Atmospheric block model. 


critical speed, too, there are two permissible 
modes of flow, one at a higher energy level 
(speed) than the other. This may be compared 
with conditions existing up- and downstream 
from an hydraulic jump as described by 
STOKER (1948). A necessary condition for the 
development of a jump in the atmosphere 
analogous in character with that observed in 
hydraulic flow in open channels has thus been 
derived, and a plausible internal mechanism for 
producing the required reduction in speed 
across the jump has been described. It must be 
emphasized however that the velocity profile 
assumed in the integration of equation (3) is 
unrealistic, and therefore that the expression 
obtained for the critical speed (equation 5) 
can not be expected to accurately describe 
actual conditions in the atmosphere. RossBy 
(1950) has also given an expression for the 
rate of propagation of such an atmospheric 
jump which moves without change of shape as 
follows: 


C= tee | N— era (6) 


where the critical ratio, N, and the expansion 
ratio, n, are defined as follows: 


N=“ and n = À (7) 


Ugo Ag 


The subscripts used in equations (7) above 
have the significance indicated in figure 7. 
In order to obtain a partial test of the theory 
just summarized daily (0300 GCT) values of 
N and n (equations 7) at the soomb level have 
been computed for the five blocking cases 


| 
| 
| 
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before mentioned. In such an actual computa- 
tion some assumption as to the boundaries 
of the westerly how must be made; in view 
of the unfortunate sensitivity of the calcula- 
tions to stream width this selection of bounda- 
ries assumes a critical importance. In this case 
the contour lines lying furthest north and 
furthest south, yet still traversing the entire 
hemisphere, were defined as the flow boun- 
daries. One must nevertheless exercise a cer- 
tain amount of subjective judgement in placing 
these boundaries and as a result the values 
obtained are only crudely approximate. In 
view of the fact that the same technique has 
been used throughout it may be assumed 
that the trends obtained, if not the numerical 
values, reflect actual conditions. 

The values obtained for the expansion ratio, 
n, (ratio of down- to upstream widths) in 
three of the five cases analyzed are shown 
graphically as time functions in figure 8. 
The three cases given illustrate the extremes 
obtained. In order to permit a better compa- 
rison, a uniform time scale is used in this 
figure. Each case is “zeroed” on “B-day”, 
the day on which blocking was initiated as 
defined previously. The time trends for the 
1949 and 1950 cases show a slow rise until B- 
day; afterwards growth continues at a some- 
what accelerated rate until blocking is fully 
developed. The January 1947 case is in reality 
a rejuvenated cycle of a blocking situation 
which began in December 1946, as is shown 
in Part II. The time trend of in such a case 
would be expected to have the form shown. 
After maximum development, except for the 
1950 case which extends into February and is 
not shown, the values of n fall rapidly as the 
block deteriorates. Also shown, for compa- 
rison, in figure 8 is the time trend of n com- 
puted for the month of January 1949, an ex- 
cellent example of a persistent zonal-type flow 
pattern in the Atlantic-European area. We 
conclude from these results that an average 
value for the ratio of down- to upstream 
widths during blocking is about 2.5 whereas at 
other times it falls to about 1.5. A fixed defini- 
tion of up- and downstream sector limits is 
permissible for all cases of Atlantic-European 
blocking since in all such cases the jet split 
remains within a relatively narrow longitudi- 
nal span centered about 10° W longitude 
(Part II). 
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The values obtained for the critical ratio, N, 
in the three cases just discussed are plotted as 
functions of time in figure 9. The critical ratio 
at any particular time is defined as the ratio 
of the observed stream velocity to the critical 
velocity, too, at that time; it may be taken as a 
measure of the stream’s susceptibility to 
blocking development. In each of the three 
cases we observe at least a secondary maximum 
in N oa B-day. It is significant to note that 
large values of N are observed also prior to 
the initiation of blocking although there is 
some suggestion of an increasing trend (except 
for the 1947 case as before) up to B-day. 
The large values of N (> 15) were computed 
for very narrow stream widths, under condi- 
tions apt to render calculations of too inaccu- 
rate. They may therefore be largely discoun- 
ted. The time trend of N during January 1949 
is again shown in figure 9 for comparison. 


Computations of the movement of the 
block using equation (6) were made in each 
of the five cases analyzed. General agreement 
between the computed and observed direction 
of movement was obtained, but the computed 
speeds in almost every instance exceed the 
actual by a factor of two or three. In view 
of the assumption of constant shape, which was 
made in the derivation of equation (6) and 
since only approximate values for N and n 
were available, it is not surprising that the 
computed motion should be in error. 


a 
EXPANSION RATIO az 
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DATE 


Fig. 8. Daily values (0300 GCT) of the expansion ratio 
computed for the 1947, 1949 and 1950 blocking cases. 
Expansion ratios for January 1949 (zonal flow) are shown 
above for comparison. Upstream sector taken as 120° W 
— 30° W. Downstream sector taken as 30° W—30° E. 
Blocking development initiated on B-day. 
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The Initiation, Development and Dissipation 
of Blocking 


It is clear that, because of the normally low 
value of the critical velocity, oo, values of the 
critical ratio, N, normally exceed unity. The 
existence of supercritical velocities in the upper 
westerlies cannot therefore be used as a suf- 
ficient criterion for the initiation of blocking. 
It has been shown for the case of hydraulic 
jumps that such developments must be ini- 
tiated through the action of finite external 
impulses on the system. Further, it has been 
experimentally demonstrated by RUSSEL (1936) 
that the minimum energy required in any 
such impulse to precipitate a complete break- 
down of the flow in an open channel decreases 
rapidly as the ratio, N, increases. Accordingly, 
only extremely deep, slow-moving stream 
systems are stable to impulses of any intensity; 
whereas decreasing depth and increasing spee 
lower the system’s “‘threshold of sensitivity’. 
That is, successively less intense impulses 
are sufficient to initiate a breakdown of the 
simple flow pattern (an hydraulic jump results). 

In the atmosphere the westerly jet is sub- 
jected to an almost continuous barrage of 
external impulses of varying intensity. If 
one accepts the concept of lateral friction as a 
driving mechanism of importance in pro- 
ducing and maintaining the westerly jet, 
(RossBY 1947), it must be assumed that acce- 
lerations thereby produced in the jet stream 
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Fig. 9. Daily values (0300 GCT) of the critical ratio 
computed for the 1947, 1949 and 1950 blocking cases. 
Critical ratios for January 1949 (zonal flow) are shown 
above for comparison. Upstream sector taken as 120° W 
—30° W. Downstream sector taken as 30° W—30° E. 
Blocking development initiated on B-day. 


DANIEL F.4R EX 


~ 


are neither constant in direction nor in magni- 
tude, but rather vary markedly as mid-tro- 
pospheric cyclonic vortices move in space 
and vary in intensity. Regions of forced local 
acceleration are thus produced in the westerly 
stream which will in turn tend to produce 
difluent flow in adjacent areas immediately 
downstream. In such difluent zones one would 
expect the initiation of a momentum jump to 
be favored through the action of the pressure 
excess, p’, defined before (equation 1). In the 
colour plates of the June 1949 development 

igure 2) we note that the initial appearance 
and subsequent rapid development of the 
blocking ridge took place just downstream 
from an area exposed to the effects of an 
intense cyclonic vortex lying just to the north 
of the jet stream itself. Such a relationship 
has been observed in each case of blocking 
action inspected. We therefore postulate that 
the initiation of a block is impossible down- 
stream from any region in which the critical 
ratio, N, is less than unity and further, that 
such developments will be initiated whenever 
(and wherever) a sufficiently intense “cyclonic” 
impulse, for the then existing ‘critical ratio, 
is imparted to the stream. The appearance of 
blocking is therefore more probable when a 
strong and relatively narrow westerly current 
system aloft prevails over the upstream area; 
such a flow pattern is often associated with 
confluence as described by Namias and Crapp 
(1949). 

If the above-described theory is substan- 
tially correct we should expect to observe, at 
the “instant” of initiation in past cases of 
blocking, an inverse relationship between the 
magnitude of impulse intensity, I, and the 
critical ratio, N. In table III below, the values 
of I and N, observed in the 1947, 1949 and 
1950 cases of blocking, are given. 


Table III 
ees Critical Impulse Date of 
Ratio — N | Intensity — I Initiation 
1949 3.6 + 0.39 K 16 June 
1950 9.0 , +017 K 14 January 
1947 10.3 + o.11 K 17 January 


I= SK where S is defined by equation (8) and 
K is a constant of proportionality.) 


The date of initiation was determined in 
accordance with the definition of blocking 
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given before. The impulse intensity, for want 
of a better definition, has been assumed to be 
proportional to the average of the mean hori- 
zontal geostrophic wind shears, computed 
normal to the jet’s northern boundary on the 
soomb surface at three separate points around 
the cyclonic vortex lying nearest upstream 
from the developing blocking ridge. If the 
contour height in dynamic meters at a point, 
B, on the northern boundary of the jet is 
denoted as B, then the mean horizontal geo- 
strophic wind shear, S, between the points, 
A and C, is given as follows: 


ps. FH (8) 


where A and C denote the contour heights in 
dynamic meters at the points, A and C, si- 
tuated along the normal line at a distance, L, 
(in nautical miles) respectively north and south 
of point B. In the computation of table III 
the distance, L, has been taken as 300 nautical 
miles. Shear values computed from equation 
(8) are given in meters / second / 10 nautical 
miles; positive values indicating anticyclonic 
shear across the northern jet boundary at 
point B. We observe that the three values 
of N and I shown in table III are inversely 
related as we should expect. 

It may be of some interest to speculate upon 
the seasonal and geographic variations which 
block development should be expected to 
exhibit in accordance with this theory. With 
regard to seasonal influences we may state 
from observational evidence that during the 
summer the intensity of mid-tropospheric 
cyclonic systems, the average jet speed and 
the average width of westerly currents will 
be relatively diminished. We cannot therefore 
state directly any conclusion as to the seasonal 
frequency of block occurrence, since decreased 
impulse intensity is associated with an indeter- 
minate change in N. Evidence presented in 
Part II shows that the maximum frequency of 
occurrence is actually observed during the 
late winter or early spring. 

With regard to geographical influences we 
recall the observed fact that the strongest 
cyclonic vortices found aloft over the northern 
hemisphere (and within the sector 120° W— 
60° E) are most frequently situated over the 
eastern seaboard of North America and wes- 
tern Atlantic in association with strong polar 
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outbreaks down the Mississippi Valley and 
across the eastern United States. The collapse 
of such cold anticyclones in the lower tropo- 
sphere and their röle in the initiation of mid- 
tropospheric cyclonic systems has been dis- 
cussed by RossBy (1949). As a result we should 
expect to find a “preferred region”, of most 
frequent block initiation, over the eastern 
Atlantic. A similar preferred region, associated 
with the eastern Asiatic continent, should be 
expected to appear over the east-central Pacific 
and would produce blocking situations which 
would at times affect western North America. 
In Part II observational evidence substantiating 
the existence of these preferred regions is 
presented in agreement with the prior findings 
of ELLIOTT and SMITH (1949). 

The development of blocking. situations, 
both in meridional and longitudinal extent, 
can not be described in detail until the dynamics 
of the macro-turbulent zone itself are better 
understood. Meridional growth is most com- 
monly accomplished through the intensifica- 
tion and northward extension of the block- 
ing anticyclone aloft and is accompanied by 
a general increase in the amplitude of meri- 
dional flow components down stream. Meri- 
dional growth is typically followed by an 
extension of the blocked area downstream 
(longitudinally) through the addition to the 
train of more cut-off vortices. This train of 
cyclonic and anti-cyclonic cells may often 
extend more than halfway around the hemis- 
phere and in extreme cases may transform 
the entire hemispheric circulation pattern into 
a purely cellular type. Figure 10, showing the 
transformations produced in the hemispheric 
westerly pattern during the 1949 case pre- 
viously described, illustrates such an example 
of extended growth. It should also be men- 
tioned that the growth rate of blocks is often 
not uniform, but rather the development may 
pass through several cycles. The 1947 case, 
previously mentioned, illustrates the rejuvena- 
tion of a blocking situation which was initiated 
many weeks earlier. From observational 
evidence we may state that large values of 
the critical ratio, N, are associated with (favor) 
accelerated growth. The motion of blocks is, 
in general, irregular although a trend toward 
retrograde motion during the developing 
stages and toward progressive motion during 
dissipation may be detected. In general re- 
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trograde motion is associated with large values 
of N and vice versa. 

The eventual dissipation of blocking action 
has been preceded in each case analyzed by a 
decrease in the critical ratio to a very low 
value, less than or near to unity. We note in 
figure 9 that in the dissipating stages of both 
the 1947 and 1949 cases, values of N less than 
unity were observed. The critical ratio may 
be lowered by increasing the mean stream 
width and/or by decreasing the mean stream 
velocity in the upstream sector. Either or 
both effects are commonly produced by 
changes occurring in the lower troposphere 
over the upstream zone, i.e. a relaxed tem- 
perature gradient, decreased vertical stability, 
decreased confluence, etc. Occasionally the 
blocked region may extend so far down- 
stream as to affect its own upstream area and, 
in so doing, may destroy itself. An example of 
such a self-destructive block is shown in figure 
10 already referred to. 

In concluding this section, it may be desi- 
rable to discuss briefly two variations of the 
typical single warm-anticyclone blocking si- 
tuation. Blocking has previously been descri- 
bed and illustrated as originating with a warm 
anticyclone situated immediately downstream 
from the jet split. Instead in many cases the 
first vortical element in the developing block is 
a cold cut-off cyclone. Such a case in its more 
mature stages is illustrated in figure 11. In 
other cases, which were initiated with a typi- 
cal warm anticyclone, cold vortices are, during 
the life of the block, cut-off in advance of the 
blocking anticyclone and included in the 
macro-turbulent zone — see BERGGREN, BOLIN 
and Rosssy (1949), figure 26. In all cases how- 
ever, an anticyclone tends, in the mean, to 
dominate the area lying just downstream from 
the jet split due to the more stable character 
of a warm vortex. The anticyclone shown in 
figure 11 has partially replaced the already 
dissipating cyclone. 

Although it is more usual to find only one 
block in the hemispheric flow pattern, at 
times two blocking situations may exist 
simultaneously. Such a situation is illustrated 
in figure 12. Perhaps due to the type of self- 
interference already referred to such double 
blocks are normally short-lived. One or the 
other dissipates, leaving a single block which 
may persist for several weeks. 
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Fig. 10. An example of extended blocking development. 

Absolute topography of the hemispheric westerlies is 

shown at 10-day intervals. Contour heights in dynamic 
decameters. 
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Fig. 11. An example of a cold cyclone blocking case. 
Movement of troughs retarded as they approach the 
block; present position — double, preceding 24hrs. 
— dashed and preceding 48hrs. — solid. Pressure field 
over the hemispheric westerlies in millibars. 


Conclusion 


Certain apparent differences between the 
development of wave-like perturbations in 
the zonal current and the development of 
atmospheric blocks may be pointed out. 
Both processes produce strengthened north- 
south flow components and increased meri- 
dional exchange in the area affected. Either, 
under favorable conditions, may produce a 
complete breakdown of the hemispheric zonal 
westerlies. Wave development is typically 
associated with changes in the stream channel 
which are continuous both in space and time. 
The process, if continued, may lead to the 
formation of one or more cut-off cyclonic 
and/or anticyclonic vortices as has been 
described by PALMÉN (1949). The appearance 
of such vortices does not in general produce a 
flow pattern which satisfies the definition of 
blocking action as given herein; they com- 
monly appear quite independently of blocking 
phenomena. An example of the transforma- 
tion of the hemispheric circulation pattern 
into a distinctly cellular type, through the 
action of developing waves, is shown in 
figure 13. 

The development of a block in the upper 
westerlies is typically associated with the 
appearance of a persistent quasi-stationary 
discontinuity in the flow characteristics along 
the stream channel. Upstream from this 
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discontinuity, or “shock”, an essentially zonal 
flow at a high energy level (speed) is obser- 
ved; downstream the flow becomes markedly 
meridional at a lower energy level. A compa- 
rison of figure 10 with figure 13 will illustrate 
these apparent differences. Actually no sharp 
distinction may exist between these two types 
of development since in a 2-dimensional treat- 
ment both may be solutions, although of quite 


different appearance, to “ =o where the 


dt 

vertical component of absolute vorticity is 
represented as Z,. It is hoped that such a 
solution for a representative blocking deve- 
lopment may be undertaken with the elec. 
tronic computer; the results of such a computa- 
tion would permit a more rigorous comparison 
of wave and block development. 

In a recent paper RIEHI, YEH and La SEUR 
(1949) have computed time sections of the 
change in angular momentum at various levels 
in the upper air. They have shown that the 
transfer of angular momentum across the 
middle latitudes usually takes place from south 
(tropics) to north (arctic) but that during 
certain periods the direction of this transfer 
may be reversed. It is interesting to note that 
periods of blocking, as taken from the block- 
ing-case catalog, Part II, correlate well in time 
with periods during which a north to south 
transfer of angular momentum is observed. 
For example note the Au-momentum trends 
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Fig. 12. An example of a double block with the Atlantic 

“shock front” shown at ”A’” and the Pacific “shock 

front” shown at “’P’’. Pressure field over the hemispheric 
westerlies in millibars. 
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Fig. 13. An example of the transformation of the hemispheric circulation pattern into a cellular type through the 
action of wave development. Absolute topography of the hemispheric westerlies is shown at 4-day intervals. 
Contour heights in dynamic decameters. 


shown by Rrient, YEH and LA SEUR (1949, 
figs. 6 and 7) during February—March 1948 
and January—February 1949 (blocking cases 
discussed previously). The destruction of angu- 
lar momentum in the mid-troposphere is no 
doubt accomplished more effectively through 
the action of an atmospheric block than by 
any other available means. The development 
of such blocks acts as a kind of “safety-valve” 
on the atmospheric motion in that it limits 
the maximum value of angular momentum 
permissible in the upper westerlies. Blocking 
action therefore plays a fundamental réle in 


atmospheric dynamics in that it limits the accu- 
mulation of angular momentum in the upper 
westerlies and provides a mechanism for the 
southward transport of angular momentum 
across mid-latitudes. 

Various facts of observation relating to the . 
blocking process have been presented in this 
paper. Although no adequate theory of 
blocking is at present available an attempt 
has been made herein to fit these observations 
into a tentative theory and to discuss the 
blocking process as an element of the general 
circulation of the atmosphere. 
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Abstract 


A method for the determination of areas of rising and falling pressure in the free atmosphere 
is given, together with a brief discussion of how these pressure variations lead to the growth of 


isallobaric areas at sea level. 


I. Introduction 


The upper-air service in Norway was re- 
organized after the war by SVERRE PETTERSSEN. 
The current work is now performed along 
the lines described by S. PETTERSSEN and C. H. 
B. PRIESTLEY : “Preparation of Upper Air Anal- 

| a : : 
ysis and Forecast Charts.” (Reprinted in 
Norway in a somewhat modified form as: 
Vertjenestememorandum nr. 2.) We are 
making use of the rule of SUTCLIFFE, that the 
thickness lines approximately move with the 
speed of the gradient wind, and of the fol- 
lowing rules of PETTERSSEN: 

“(r) Where the thickness lines climb up the 
pressure surface the pressure must fall, and the 
fall increases with the isobaric thermal wind 
and the slope of the pressure surface. 

(2) Where the thickness lines are parallel to 
the contour lines (i.e. no climb) the pressure 
remains stationary. | 

(3) Where the thickness lines descend the 
pressure surface, the pressure must rise, and 
the rise increases with the isobaric thermal 
wind and the slope of the pressure surface. 

In applying these rules, due allowance must 
be made for the non-conservation of the thick- 
ness in regions where there is indication of 
subsidence or ascending motion, and also in 


regions where the air is being appreciably 
heated, or cooled, by the underlying surface. 
In particular, in outbreaks of polar air the non- 
adiabatic heating will be appreciable in the 
lower troposphere, and the change of the 
thickness from 1,000 to 700 mb. is often not 
higher than so per cent of that which would 
be expected of sheer advection.” 

In the following paragraphs we shall study 
such weather situations where the upper-air 
charts show a distinct grid between contour 
lines and thickness lines, as in the rule of 
PETTERSSEN, but we shall consider the problem 
from a different point of view. 


II. Variations in the height of the upper- 
air surfaces of pressure 


We shall try to judge these variations in 
height from the movement of the thickness 
lines and from the barometric tendencies at 
sea level. As the tendencies are valid for three 
hours, we must determine the movement of 
the thickness lines for the same interval. At 
the weather service in Oslo only the 700 and 
soo mb. upper-air charts are drawn, and only 
twice a day, at 4h and 16 G.M. T. 
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. Intersecting contour lines and thickness. lines. 
Further explanation in the text. 


We shall first study the variation in the 
height of the 700 mb. surface by means of 
some examples. In order to simplify our 
discussion, we assume that the contour lines 
and the thickness lines are straight, fig. 1. As 
a first illustration we assume that the baro- 
metric tendency has been zero at the ground 
during the last three hours. We may then, 
with sufficient accuracy, say that the 1,000 mb. 
surface has been stationary during this interval 
of time. With the help of the “‘gradient ruler” 
we can determine how far the thickness line 
through point a, fig. 1, has moved in three 
hours. Let us say it has reached the position 
of the foregoing thickness line through point 
b; then the 700 mb. surface has sunk 40 m at 
point b during three hours (assuming that the 
thickness lines are drawn at intervals of 40 m 
and that cold air is moving forward). 

If a point a in the 1,000 mb. surface has 
moved to b, the corresponding point in the 
700 mb. surface has moved to b’. A column 
of air through point a, originally vertical, can- 
not remain vertical throughout the movement. 
But, during a short interval of time we may, 
as a first approximation, think of the columns 
as remaining vertical during the movement, 
at least as far as the single layers between the 
chosen pressure surfaces are concerned. 

As a second illustration we assume that the 
surface barometric tendency has been + 1.0 
mb. over the area in question. See fig. 1. 
Then the 1,000 mb. contour line abc has 
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moved 1/5 of the distance between c and d, 
i.e. c has moved to c’ (approximately). At 
the point b the 1,000 mb. surface consequently 
has risen 8 meters. Now let the thickness line 
through point a reach b in 3 hours. The fall of 
the 700 mb. surface at b will then be 4o—8 = 
= 32 m. If the barometric tendency had been 
+ 5 mb. instead of 1, in the same manner we 
would have found, that there had been no 
change in the height of the 700 mb. surface. 
Further if the tendency had been greater than 
s mb., we would have found a rise instead 
of a fall of the 700 mb. surface, assuming the 
speed of the thickness lines has been the same 
in all cases. 

In a similar manner we can determine the 
change of height when warm air is streaming 
forwards. 

If we possessed absolutely reliable charts, we 
could determine the rise or fall of the 700 mb. 
surface over any chosen point at sea level by 
the method described above. In the same 
manner we could later determine the rise and 
fall of the soo mb. surface relative to the 
700 mb. surface and hence to the sea level (or 
1,000 mb. surface). With current charts we 
can only hope to determine the main features 
of the isallobaric areas in the free atmosphere, 
and only during great changes in the weather 
situation. 

When the rising and falling areas have been 
determined — as far as possible — we can 
take into account accelerations in the upper air, 
which we have put aside up to now. We can 
thus in a qualitative way decide if a decrease 
or increase or no change in the sea level 
isallobaric areas is to be expected. In order to 

et a better understanding of such accelera- 
tions we shall study them with the help of 
some simple computations. 


III. Movement of a particle of air in a 
variable field of pressure. 


In the following discussion we can only 
make use of point dynamic relationships and 
can take into account only horizontal move- 
ments. In order to make the calculations as 
simple as possible we shall consider only 
straight, parallel and equidistant isobars and 
isallobars. In general the isobars and isallobars 
will intersect, and we shall take the x-axis 


in the direction of the isallobars, see fig. 2. 
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G is the gradient force with components Gy 
and Gy. The change in the isallobaric gradient 
in { seconds is G, ¢. In fig. 2 this vector is 
chosen vertical to the x-axis. The specific 
volume & and the Coriolis parameter 
A=2 sin are regarded as constants. The 
horizontal components of the wind velocity 
are u and v. Regarding a unit mass, we have 
the following equations of motion: 


d?x dy a 

a EE CG (1) 
d2 dx 
ee 6 


Potime = 4, and 1 = vy by fap we 
have the solutions: 


| (1 + =) cos À f + 
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Fig. 2. Coordinate system with the x-axis along the 
isallobars. : 


In these equations the first square bracket 
shows the influence of the original field of 
pressure, and the second square bracket gives 
the influence of the isallobars. The last term in 
the first square bracket represents the geo- 
strophic wind, the two first terms giving the 
oscillations caused by deviations from the 
geostrophic wind at t = o. Generally we do not 
know these deviations, we must use our few 
direct, wind observations to correct the pressure 
gradient, when drawing the contour lines. 
Since the geostrophic wind itself cannot con- 
tribute to the development of rising or falling 
areas, we must limit our research to the in- 
fluence of the isallobars, i.e. we may only 
make use of the second square bracket in each 
of the equations (3) to (6). In the two last 
equations we can determine C, and C, by 
putting x and y=o, and so we get the 
following expressions: 


aG, 


= Ti sin À t — ra t (7) 
v= a cos At — a (8) 
n= Fh cos At a + LE (9) 
de aM sin At — a t (ro) 
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We cannot expect exact results with so 
many simplifications, but practice seems to 
indicate a marked influence of the isallobars, 
when large changes in the field of pressure 
take place. 


Table ı 
Art cm/sec km 
in t — | — 
degrees u v N | y 
se 14.5 min] 0.025 0.35 

HOP 29 » 0.06 I.I | 0.0003|0.006 
Se 44 > 0.3 2.5 | 0.0017|0.022 
207 58 » 0.5 4.4 | 0.0044]0.051 
25° |r hour 13 » 1.0 6.8 [0.011 |0.099 
ro LE N. » EF 0:70 | 0/0239104172. 
Ae a By eT » 2.7 13.1 | 0.048 [0.271 
AO Es » 4.0 17.0 | 0.071 [0.401 
er Pe PO CEE » ely 21.3511 0.203 1055068 
Sor 2, DES » 77 25.9 | 04172 10.774 
5S 12 ie » | 70,2 31.0 | 0.249 |1.022 
OOo ESS pe | rs 30.3. 0,3 Sa maus 
05 1300 Ne 16:6 41.9 | 0.480 [1.656 
FO» |3. ide » | 20.5 47.8 | 0.64I |2.047 
RSS OS 38 » 124.9 53.8 | 0.839 [2.491 
SOUS 2? 5B » 129.9 60.0 | 1.078 |2.987 
Ss tl Oa gore) | 35-4 66.3 | 1.362 |3.538 
OO rss pee: » 141.4 72.6 | 1.691 |4.144 
OS OS NET » 148 79 2.100 |4.806 
100° 14 » ST » 55 85 25:90 WSS 21 
HOSS T'S) ee 5 » | 63 91 3.018 |6.293 
HO IS ae 20 POlE7E 97 3.639 17.116 
RES WS. 40350 180 103 4.338 [7.992 
CR 2 102 IE |1o9 |4.966 |8.918 


In order to compute u, v, x and y we will 
put « = 784 cm, i. e.airato°C and 1,000 mb., 
Ve peitoswane G7 mb./1,000 kn in3 
hours, i. e. the steepness of the pressure surface 
has changed 8 m/r,000 km in 3 hours. The 
result is found in table 1 and figs. 3 and 4. 
Since u, v, x and y are proportional to G,, 
we have only to multiply by a factor to find 
the corresponding result for other values of 
G, and further since « G, is approximately 
constant in the atmosphere, we may also use 
sea-level calculations for any surface in the 
upper air. 

We shall use these computations to form 
an estimate of the influence of the isallobaric 
areas, assuming that the isallobars can be drawn 
with great accuracy for a short interval of 
time. 

In a previous paper it is shown that we 
can use the length of the isobars to determine 
the mean pressure gradient and the mean 


215 


3! 53 min 


2600 
3! TETE 


2200 
1600 
1400—  } 


1000 


600 


100 500 1000 m 


Fig. 3. Path of a particle moving under the influence of 
straight, parallel isallobars. 
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Fig. 4. Same path as in fig. 3, but for a longer interval of 
time. 
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geostrophic mass transport in a horizontal unit 
layer. If we divide the intermediate space 
between two neighbouring isobars into squares 
— “isobaric squares” — the kinetic energy is 
the same for each square, because the wind 
speed is inversely proportional to the distance 
between adjacent isobars. Since this last con- 
dition holds for isallobars too — equations (7) 
and (8) — we deduce that v - As (and u - As) is 
constant for each isallobaric square, where 
As is the length of a side. Counting the number 
of squares we may thus calculate the horizontal 
flow across the isallobars and therefrom deter- 
mine the horizontal divergence. (We may 
also calculate the wind circulation around an 
isallobar, and the mean mass transport produced 
by the isallobars; the mean isallobaric gradient 
may be determined by measuring the length 
of the isallobars. Here only the divergence 
shall concern us.) In order to get an acceptable 
approximation, we must use a time interval so 
short that an arbitrary particle of air can be 
supposed to stay within an isallobaric square 
during this interval. If the geostrophic wind is 
strong at f— 0, we will have to choose a 
time interval considerably shorter than 3 hours. 
At present we have only vague ideas of the 
isallobaric field and, therefore, cannot yet make 
use of a method based on the isallobaric 
squares in attempting a computation of the 
mass transport across the isallobars. We must 
content ourselves with a more qualitative 
determination. More accurate charts in the 
future will surely lead to a more accurate 
determination of this mass transport. 

Now let the schematical fig. 5 represent a 


Fig. 5. Rising and falling area in an upper-air pressure 
surface, schematical. 
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rising and falling area in the — say — soo mb. 
surface. From the rising area air will flow out 
in every direction and into the falling area air 
will flow in from every direction. We should 
consequently expect to find a falling area at 
sea level under the area of a rising upper-air 
surface, and a rising area at sea level under the 
area of a falling upper-air surface. The in- or 
outflow depends not only upon the isallobaric 
gradient, but also very much upon the time 
interval chosen for the isallobars, as is shown 
by table 1. For instance: Small tendencies at 
sea level coincident with strong isallobaric 
areas in the upper air indicate that the flow 
across the isallobars in the upper air has gone 
on only for a short time, and that an increase 
in this flow is therefore to be expected, which 
will in turn lead to a strengthening of the 
sea-level isallobaric areas. 


IV. Further remarks on thickness lines 


and isallobars 


The rising or falling of the upper-air surfaces 
depends upon the sea level tendencies and upon 
the movement of the thickness lines. To 
estimate the movement of the thickness lines 
is, however, no easy task, as we cannot make 
use of the rule of SuTCLIFFE in a servile way. 
Fig. 6 shows the 700 mb. chart of 4/11 1948, 
o4h G. M. T. The two thickness lines AB and 
CD are lying approximately motionless at the 
intersections. The speed of the gradient wind 
north-east of Iceland is indicated by the arrow, 
but it is not likely, that a corresponding dis- 
placement of the thickness line AB will take 
place in this region. The cold stream of air 
is relatively narrow, and we cannot expect 
that the cold air will move within the contour 
lines as in solid tubes. It will try to sink and 
spread out along the ground, and we therefore 
should expect a retardation in the speed of 
the thickness lines. But how great this retarda- 
tion will be one can only guess. In the stream 
under consideration (near the arrow) I should 
probably have said at least 30 %. Retardation 
resulting from warming at the sea-surface is 
also to be taken into account. PETTERSSEN and 
PRIESTLEY (l.c.) have pointed out that we may 
expect a retardation up to so %, if cold air 
moves across a warm sea-surface. In cold 
arctic or polar air moving from Labrador 
towards the southeast the thickness lines at 
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Fig. 6. 700 mb chart of Nov. 4, 1948, 04h G.M.T. 


times appear to be retarded even more than 
so %. In the Norwegian Sea I should estimate 
that 25 °4 would be a more reasonable value, 
but accurate estimates cannot be made at 
present owing to the scarcity of observations 
both in space and in time. The wind speed and 
the difference between air and sea temperature 
must play a röle too. Cooling from the ground 
is of lesser importance than warming. A 
cooling from sea level will not effect the height 
of the 700 mb. surface in any marked way 
within the short interval of 3 hours. 

The heat of condensation may also effect 
the height of the upper-air surfaces in a very 
pronounced way, when we consider great 
areas of precipitation such as in warm-frontal 
rain areas. If the rainfall amounts to I mm an 
hour and the condensation takes place in a 
vertical column between 2,000 and 6,000 m, 
this column will expand about 30 m in 3 hours. 
We should expect therefore a rapid rising of 
the upper-air surfaces over and near any large 
area of heavy precipitation and that outflow 


from this area would consequently produce 
a strong falling area at sea level. In fact the 
most marked falling areas at sea level seem to 
be closely associated with areas of heavy rain. 
If we therefore find it likely that a rain area will 
increase, we may also conclude that an increase 
of the falling area at sea level will probably take 
place. Such a development must at present be 
forecast by indirect methods. The upper air 
humidity observations are too scarce to permit 
drawing of any conclusions. 

If the thickness lines and the contour lines 
are running parallel, no great changes can take 
place in the height of the upper-air surfaces. 
The changes are then determined by the sea 
level tendencies, but in such cases the sea level 
tendencies are apt to be small too. Distinct 
grids between thickness and contour lines may 
be found in the vicinity of strong cols, 
where cold and warm air masses are drawing 
together. Fig. 7 shows a schematical example 
of such a col. The thick, even lines represent 
1,000 mb. contours, the broken lines are lines 
of equal thickness between this surface and 
the 700 mb. surface. These lines have been 
drawn arbitrarily in the figure and afterwards 
the 700 mb. lines have been constructed in the 
usual way. Cold air is streaming southwards 
across the Norwegian Sea, warm air is stream- 
ing northwards in the Mid-Atlantic. This 
movement will lead to an accumulation of 
thickness lines near the col and gradually create 
a front there. 

A vertical plane through the straight line 
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Fig. 7. Schematical illustration of two pressure surfaces 
near a col. 
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‘Fig. 8. Vertical section through the line AB in fig. 7. 


AB, fig. 7, will intersect the pressure surfaces. 
In fig. 8 ACB represents its intersection with 
the 700 mb. surface, and DEF with a pressure 
surface near the ground. Let us assume — 
in order to simplify the case — that the sea- 
level tendencies are small enough to be neg- 
lected. The 7co mb. surface will then rise, 
whenever the thickness lines are moving 
northwards in the warm air, whereas it will 
fall over the cold air moving southwards. The 
line ACB, fig. 8, will gradually change into 
the line A’CB’. E represents the col near the 
ground, C the col at the 700 mb. level. The 
isallobaric areas, which are created in the 700 
mb. surface, will therefore cause a mass 
transport across the isallobars. The movement 
will be at first slow but will steadily increase in 
speed. See table r. 

In the same manner we may study the soo 
mb. surface and determine how this surface 
will rise or sink in relation to the 700 mb. 
surface. In most cases the two surfaces will 
move qualitatively in the same way, if great 
air masses are involved in the process. 

As soon as outflow from the rising areas 
in the warm air is perceptible, a falling area 
will develop below, at sea level, and as the 
outflow in the upper-air is accelerated, the 
falling area at sea level will increase. Similarly 
we should expect an increase in flow into the 
upper area over the cold air atid at sea level 
an increasing rising area. 

The utility of upper-air charts is therefore 
evident in cases with zero tendencies at sea 
level, and with distinct grids between thickness 
lines and contour lies In these cases the 
development of isallobaric areas at sea level 
may be foreseen or at least suspected when 


Fig. 9a. Schematical illustration of “open waves” on the 
contour lines near a col at the beginning of a cyclone 
cycle. 
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Fig. 9c. 


Illustrating the form of the thickness lines of a 
great central cyclone. 


the sea-level chart alone would give no sure 
indications. 


In some situations the thickness lines may 
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initially have the shape of broad, open “waves”, 
see fig. 7. As the warm air moves northwards 
and the cold air southwards, the amplitude of 
these “‘waves’’ will increase. Fig. 9a shows 
a thickness line at the beginning of the process 
and fig. 9 b in a later phase. During phase b 
we must imagine that there is a well developed 
falling area at sea level in the warm-air tongue. 
This falling area will accelerate the air in the 
ground layers, and a flow into the falling area 
will take place there. This inflow will partially 
counteract the building of the falling area, but 
the magnitude of this reaction is not easy to 
estimate. In fig. 9 b it is suggested that cold 
air is being drawn into the falling area from 
the northwest. If this cold air is intersected 
by many thickness lines, a sinking of the 
upper-air surfaces is likely to take place during 
the movement. This may cause an increased 
outflow from the warm-air tongue towards 
the west. We may then have an increasing 
sea-level rising area in the cold air west of the 
tongue, a continuation of the falling area in 
the tongue, but a decrease in the sea-level 
rising area east of the tongue, for the reason 
that more air is now moving westwards from 
the upper-air area in the warm tongue. We 
cannot expect that the process will always 
follow the same scheme, each case must be 
studied individually. The warm tongue will 
at any rate have a tendency gradually to be 
longer and narrower and to disappear towards 
the end of the process (together with the 
dying cyclone). 

The falling area at sea level in the warm 
tongue will try to build a cyclone there or 
to move an old cyclone in that direction. As 
we at the same time will have a rising area 
at sea level in the cold air east of the col, we 
may point to a practical rule given in a 
previous paper: A deepening cyclone will 
generally have a pronounced rising area east 
of the cyclone. This rule was developed at a 
time when no upper-air charts were available. 
By means of such charts we now have a 
better understanding of what takes place. 

A cyclone may ultimately develop into a 
great central cyclone. In this case the warm, 
upper-air tongue may take the shape suggested 
in fig. 9c; the warm air has a tendency to 
circle around a central cyclone. A small dying 
cyclone may, on the other hand, be carried 
away by the main upper-air stream and the 
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thickness lines will then take the shape of a 
relatively narrow warm tongue above the 
occlusion. This warm tongue seems to dis- 
appear simultaneously with the disapperance of 
the rain area. 

The building, development and death of a 
cyclone may require many days, and we can 
generally make a forecast for only a part of 
the process, at any given time. We may make 
use of the rising and falling of the upper-air 
pressure surfaces to get a qualitative view of 
the process and to forecast an increase, a 
decrease or no-change in the sea-level isallo- 
baric areas. Since the method is based on the 
movement of the thickness lines, it is generally 
very difficult to make forecasts in this way for 
more than 12 hours at a time. This is easily 
understood from a look at fig. 7. North of a 
line Tromso—Jan Mayen we have generally 
no observations and should be forced to draw 
conclusions from the extrapolated thickness 
lines in this region. We have similar difficulties 
in the Atlantic Ocean. In some cases it is, 
however, possible to suggest the pattern of 
development over periods even longer than 
one day, as we have experience with the sea- 
level charts to assist us. 


Example 


Conclusions, drawn from upper-air charts, 
will generally confirm what may be already 
expected from the sea-level chart together 
with an application of indirect aerology. We 
may sometimes find cases however, where the 
upper-air charts furnish us with new informa- 
tion not given by the sea level chart. An 
example from my cwn experience will be 
given below. 

On the afternoon of the 7th Oct. 1948 the 
13-hour sea-level chart had to be used instead 
of the 16-hour chart, because new ship observa- 
tions were lacking from the Atlantic, where 
the development was of most importance. 
The sea-level chart, fig. toa, led to the 
following conclusions: 

The small low near Iceland had not deepened 
much during its existence in the last 15 hours, 
and it was supposed to be of minor importance. 
The great low in the Mid-Atlantic had existed 
for several days. It had been moving slowly 
from the west, and its falling areas had 


apparently weakened. We then probably 
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Fig. 10c. $oo mb chart of Oct. 7, 1948, 16h G.M.T. 


would expect a slow-moving central cyclone, 
and the falling areas would therefore have a 
tendency to circle around the center and to 
gradually decrease in size and intensity. The 
rising area cast of Iceland was suspiciously 
strong, but it was supposed to diminish as the 
central cyclone entered its dying phase. 

The upper-air charts, fig. tob and toc, 
pointed to quite a different conclusions as to 
the cyclone in the Mid-Atlantic. Northwest, 
west and southwest of Iceland we had a good 
grid between contour lines and thickness lines, 
and because the falling areas at sea level 
were mostly small, we should expect a pro- 


Fig. tod. Sea-level chart of Oct. 8, 1948, 072 G.M.T. 


nounced rising of the upper-air surface in this 
region, which would soon lead to an increase 
in the falling areas at sea level. The rising area 
from Jan Mayen to Scotland would probably 
not diminish, but move eastwards, because 
there was a distinct grid between contour lines 
and thickness lines over the Norwegian Sea 
and only little change in the sea-level pressure 
over northern Norway. — Fig. 10d shows 
the weather situation at 07% G. M. T. 8/10 1948. 
We see a strong falling. area southwest of 
Iceland and a strong rising area near the Lo- 
foten Isles. 

It must be emphasized that clear-cut cases 
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— such as that shown in figs. 10a to fig. 10d of practical use. It is necessary however to 
— may be comparatively seldom seen, but in learn their limitations and proper application 
my Own experience as a weather forecaster I through practice. 

have found the above-mentioned ideas to be 
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Abstract 


The following article deals with an investigation of the local variations in the horizontal 
temperature distribution in the Uppsala district during different times of the year, of the 
day and under changing weather conditions. The article treats mainly the methods of 
investigation, especially the choice of observation stations, but also comments briefly 
upon some of the results. An attempt is made to place the established temperature 
contrasts in relation with simple climatological elements such as cloudiness, wind veloc- 
ity, temperature and absolute humidity. The article is a preliminary report of an in- 
vestigation, the complete results of which are shortly to be published in Geographica 


(Papers from the Geographical Institute, University of Uppsala, No. 22). 


Introduction 


It very often, happens that a thorough treat- 
ment of a particular problem requires a basic 
knowledge of the purely local variations of 
the climatological elements. This applies, for 
example, to problems of mechanical weather- 
ing, questions concerning the location of 
plant boundaries and plant colonies, and the 
problem of the physiological importance of 
climatological factors. 

The observation series of meteorological 
stations give at best only a good picture of 
the climate just at the place and in the location 
where the meteorological instruments are 
placed. Unfortunately they give us little 
indication of the purely local deviations in the 
climatological elements. And yet these devia- 
tions are often of such an order of magni- 
tude that one must move many scores of 
miles in order to obtain corresponding changes 
in the ‘“‘macroclimatological”. elements. Our 
only possibility of finding a solution to this 
problem is in obtaining through special in- 
vestigations an increased knowledge of the 
distinctive features of local climate. 


A determination of the influence of urban 
areas on the climate has for a long time and 
for many reasons been an important task of 
inquiry. A compact city area with its mass of 
houses, industrial structures, streets and parks 
constitutes a pronounced interruption in the 
natural formation of the landscape. The changes 
in the surface configuration alone are suffi- 
cient to produce a modification of the climato- 
logical elements. More significant, however, 
are the consequenses of the commercial and 
industrial activity which characterize the city. 
Combustion processes and other transforma- 
tions of energy make the city a heat generator 
of great magnitude. These processes are 
simultaneously accompanied by a considerable 
production of dust particles, which surround 
the city and influence, at first hand, the radia- 
tion balance of its climate. Not only is the 
radiation balance a typical in and near a city, 
but the other climatological elements as well 
e.g. temperature, humidity, wind and precipi- 
tation. 


The clearly noticeable influence of an urban 
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area on climate makes desirable city planning, 
which seeks to eliminate the undesirable 
effects of a community upon healthful cli- 
mate. A good knowledge of urban climate is 
also necessary for the study of other climatolo- 
gical problems. Our oldest and most valuable 
meteorological observation-series were in ge- 
neral obtained from cities. It follows that our 
understanding of secular climatic changes can 
be influenced to a certain degree by the growth 
of the urban community. It is impossible to 
discuss the magnitude of this effect without 
closer investigation. 


Methods. General 


In order to determine the distinctive fea- 
tures of the climate of a city, comparisons are 
usually made between parallel meteorological 
observation series taken inside and outside of 
the city area in question. The disadvantage 
of this method, is the difficulty of obtaining 
fully comparable values. Very seldom are 
stations located in such a way that the influence 
of the community expresses itself clearly in 
their record. In addition, there is of course 
uncertainty in the given numerical values. 


It was therefore a great step forward, when. 


SCHMIDT (1927) in Vienna and PEPPLER in 
Karlsruhe (1929) developed a method of using 
an automobile equipped with meteorologi- 
cal instruments. This method has since come 
to be used in a number of investigations in 
different countries. 

An extensive literature concerning city cli- 
matological problems has appeared during the 
last decade. ALBERT KRATZER in 1937 published 
a survey, Das Stadtklima, which gives a sum- 
mary of all literature on city climate published 
prior to that time. This work contains an 
exhaustive description of all research results 
achieved in addition to a bibliography of 250 
titles. Among the more modern works one 
should mention an investigation by MIDDLE- 
TON and Mitar (1936). In this investigation 
the temperature distribution during a winter 
night was examined along a street in Toronto 
with the help of a resistance thermometer 
attached to an automobile. The temperature 
profiles registered show, among other things, 
the importante of topography in the hori- 
zontal temperature distribution. 

During recent years two important works 


223 


have been published concerning temperature 
and precipitation conditions in the Bergen 
district in Norway; one by SPINNANGR (1942) 
and the other by GoDsk£ (1943). Finally one 
may cite an important English investigation in 
city climatology by BALCHIN AND PYE (1947), 
concerning temperature conditions in the 
Bath district. 

For the investigation of temperature condi- 
tions in the Uppsala area, the following three 
alternative methods were available: 

I: A number of self-registering instru- 
ments, thermographs could be placed at points 
inside and outside the city. 

Il: Fixed stations equipped with instruments 
designed for regular reading could be estab- 
lished at corresponding points. 

II: Observations could be made with help 
of a thermometer mounted on an automo- 
bile according to the method of Schmidt and 
Peppler. 


Also, one can conceive of a combination 
of the three methods in some suitable way. 

The advantage of the first two methods in 
comparison with the third is that one can 
obtain a larger number of simultaneous tem- 
perature determinations from separate points. 
The first method has however the disadvan- 
tage that the values obtained are relatively 
unreliable, since errors in thermograph regi- 
strations are often as great as + 0,5° C. With 
the second method it is difficult to select espe- 
cially interesting weather situations for closer 
investigation, since all readings must necessa- 
rily occur at previously determined times. 
Moreover, both of these methods possess the 
inconvenience of requiring large and expen- 
sive equipment and a number of dependable 
observers. It should be pointed out, further, 
that unless preliminary investigation is carried 
out, difficulties will arise in the selection of 
instrument station locations. In order to ob- 
tain representative data one needs a priori a 
good knowledge of the local temperature 
variations. 

Because of these circumstances, the third 
method was chosen as the only one which 
would permit a single person to obtain relati- 
vely correct results without too great an ex- 
penditure of funds. The limitations of this 
method have been reduced as far as possible, 
as will be explained later. 
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Instrumental equipment 


An electrical resistance thermometer, con- 
structed by Professor Hilding Köhler, was 
used for the observations. Fabrication of the 
instrument had been accomplished at the Cam- 
bridge Factories, in England. This thermo- 
meter satisfied all reasonable demands which 
were placed upon it. 

Fig. 1 shows a sketch of the instrument ’s 
construction, and in fig. 2 the instrument 
appears mounted on an automobile. A copper 
tube was wound in a spiral so that the coils 
formed a mantle surface to a truncated cone. 
Its wider opening was turned forward, thereby 
improving ventilation. Inside of the copper 
tube is a double, isolated platinum wire, 
which constitutes the electrical resistance wire. 
The resistance thermometer was provided 
with a double radiation shield. The dimensions 
of the outer shield are: length 18 cm and 
diameter 10 cm. The thermometer was placed 
on the automobile in such a way that its center 
was situated 150 cm above the ground and 
15 cm from the edge of the automobile roof. 

A guide wire ran from the resistance thermo- 
meter to aWheatstone bridge and galvanome- 
ter located inside of the automobile. The 
Wheatstone bridge was manufactured by the 
Rubicon Company, Philadelphia, USA. (Por- 
table Wheatstone Bridge Test Set No. 1352.) 
It is equipped with built-in batteries and the 
galvanometer is not sensitive to jans and jolts. 
The galvanometer is a powerful pointer-gal- 
vanometer with a sensitivity of one micro- 
ampere per scale millimeter and a period of 
two seconds. In the bridge are four resistors 
of altogether 9 x (1000 + 100 + 10 + 1) 
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Fig. 1. Sketch of the construction of the resistance 


thermometer. A coiled copper tube. B electric wires to 
the Wheatstone Bridge. C inner radiation shield. 
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The instrument mounted on the automobile 
roof. 


ohms and a multiplier instrument with adjust- 
ments for 0.001, 0.01. 0.1, I, 10, 100 and 1000. 
The resistors are made of manganese and there- 
fore have an extremely low temperature co- 
efficient. The battery is a 4.5 volt dry battery 
composed of three standard cells. In the pre- 
sent investigation, one should be able to deter- 
mine the resistance used to 0.1 % or some- 
thing less, which corresponds to a tempera- 
ture uncertainty of less than 0.2° C. 

Ventilation of the thermometer was achiev- 
ed through the movement of the automobile. 
The speed of the car was kept in general at 
about 4o km/hr, but in the city it was, of 
course, often necessary to drive slower. In the 
country the speed was sometimes higher, nev- 
er greater, however, than 60 km/hr. If one 
considers a speed interval of 20—60 km/hr 
(about 6—17 m/sec.) it is obvious that the 
speeds used were fully sufficient to secure 
adequate ventilation. At the same time no 
account need be taken of dynamically condi- 
tioned temperature changes. Readings were 
always taken with the automobile in motion 
or immediately after stopping. 
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The thermometer apparatus was carefully ex- 
amined for sources of error such as the con- 
duction resistance of the guide wire, radiation 
influences, effects of heat from the automobile, 
unsuitably adapted reaction speed, insufficient 
ventilation, etc. No detailed account of these 
investigations will be presented, but it can 
be said summarily that the average error in 
the readings was estimated at about 0.1° C. 


Representativity of the observation points 


A prime problem in an investigation of 
this type is the choice of the observation 
points. Not all points in any given area will 
give reliable or representative values. A cer- 
tain variation in the reliability or representa- 
tiveness of the points chosen is unavoidable. 
One should therefore attempt a definition of 
what is meant in this case by representativity, 
so that a final selection of observation points 
can be made with the help of this definition. 

If one choses a limited terrain area not all 
points at any certain level above the earth’s 
surface have the same temperature at the 
same moment. One may assume that the value 
which at a certain moment best represents the 
temperature of the whole area at this EVE] is 
the mean value of the temperatures at a ırge 
number of points evenly distributed over the 
area. This mean value is indicated Oy and 
the temperature at a certain point P Op, 
then we have the relation: 


Op = Ox ae Ap (1) 


Ap is a quantity which indicates the deviation 
of the temperature at point P from the mean 
value for the whole area. Its value, of course, 
changes with time and varies from point to 
point. For a large number of observations, 1, 


we have: 
2Op _ 2 Ou | ZA» (2) 


n n n 


Ap 


constitutes the average value of the de- 


viation and is subsequently indicated by kp. 
kp is therefore expressed by the relation: 


n n 


kp 
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kp can of course take both positive and nega- 
tive values. A positive value kp means that 
the point P has on the average a higher tempera- 
ture than the mean value for the whole area. 
Even if a point of observation shows a large 
numerical value of kp, it can at least theore- 
tically give reliable values of the temperature 
of the whole area if one knows the magni- 
tude of kp and if Op and Oy vary ina syste- 
matic manner. 

The variations in 4p are therefore unsuit- 
able for determining the representativity of a 
point. Instead kp + dp may be substituted 
for Ap thus obtaining: 


Op = Ox ETES Op; (4) 


where kp is a constant typical for the point P 
and 6p is a variable. It follows from the defini- 
tion of kp that Z'ôp for a large number of 
observations is equal to o. 

If one has in some way determined kp, then 
a reading of the temperature at P gives a good 
understanding of the temperature of the area 
as a whole independent of the magnitude of 
kp, provided that 6p is small in magnitude. 
Indeed, one can say that the smaller the mag- 
nitude of dp, the more representative is the 
point P. Hence variations in 6p can be said to 
give a measure of the representativity of the 
point P. The statistical expression for standard 
deviation has been chosen as the quantitative 
expression for representativity. 


n 


Oo = 


(5) 


where n is the number of observations. It 
therefore holds that: the less o, the better the 
representativity. 

The above reasoning and definitions are of a 
statistical nature and may be applied without 
regard for the form or size of the area for 
which the statistical calculations are made. 
In practice the determination of the size and 
the limits of the area to be investigated is a 
problem of decisive importance. In certain 
instances a large meshed-net of reference 
areas is necessary, in others one may perhaps 
choose very small areas. 

For the subject investigation the problem 
was primarily that of finding observation 
points of two kinds: those expected to give 
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Godk. för publ. av Rikets allmänna kartverk. 


Fig. 3. General map of the area of the investigation. 

Level curves have an equidistance of 5 meters. Built- 

up area is lined. Wooded area is indicated by the 

conventional signs for forest (rings = hardwood; 
stars = pine forest). 


representative values for the city and those 
expected to give representative values for the 
surrounding countryside. In addition, the 
observations must give abundant and uniform 
numerical data suitable for statistical analysis. 

During December 1947 and January 1948 
preliminary attempts were made to find good 
points of observation. The city was criss- 
crossed in several directions and temperature 
profiles were constructed in each case for the 
area in question. In this way an approximate 
understanding of the general temperature 
distribution was obtained. On the basis of 
this experience an extensive route was selected 
which included all of the main highways from 
Uppsala and extended out 4 to. 5 kilometers 
from the city. The total length of this route 
was nearly five Swedish miles (about 30 Eng- 
lish miles) and the number of observation 
points along the whole stretch was seventy- 
nine. 

From about the middle of January 1948 
until the middle of April observations were 
made along the entire route some 20 times. 


Godk. för publ. av Rikets allmänna kartverk. 


Fig. 4. Isotherm map based on mean temperatures 
for 20 observations made during the period January— 
April 1948. 


s 


These observations were carried out at different 
times of the day and under different weather 
conditions in order to avoid a subjective pre- 
ference for any particular weather type. It 
was intended that the observations taken 
should compose as good a picture of the 
“average” weather conditions as possible. 

In arranging the material from each obser- 
vation run it was necessary to take into con- 
sideration temperature variations during the 
time that observation was in progress. One 
must, in fact, reduce all observations to the 
same point in time in order to obtain the 
temperature distribution at that particular mo- 
ment. Since every run took one and a half 
hours on the average, it is obvious that signi- 
ficant temperature changes occurred. Four 
readings were made at Stora Torget, the City 
Common, during each run (cf. the survey 
map, fig. 3). In addition, at the beginning 
and at the end of each run, temperature 
values were taken from the thermogram re- 
corded at the Meteorological Institute, 1 km 
west of the center of the city. The temperature 
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Godk. för publ. av Rikets allm. kartverk. 


Fig. s. Map showing the representativity of the ob- 
servation points. The continuous line indicates the 
route finally chosen. Representativity is indicated 
by the value of o according to the following scale 
(Greg. Si: 

0.61°—o.70°: ® 
0.71°—0.80°: @ 


0.81°—0.90°: & 
0.91 —-1.00°: & 


range was afterwards computed as the mean 
value of the temperature range at the City 
Common and at the Meteorological Insti- 
tute. A linear variation has been assumed 
wherever the thermogram showed no spe- 
cial irregularity. With a knowledge of the 
general temperature range and the time of 
each observations all observations may be 
reduced to the same point in time. It was there- 
fore possible to construct isotherm maps for 
each observational run. Examples of these 
maps will be demonstrated later. 

The next phase in the investigation of repre- 
sentativity concerns the determination of mean 
temperature at each point for all 20 series of 
observations. The result appears in graphic 
form in fig. 4. Isotherms are drawn in this 
figure with an equidistance of one half of 
one degree. The figure clearly shows a posi- 
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Godk. för publ. av Rikets allm. kartverk. 


Fig. 6. Map showing represcntativity of the observa- 
tion points when the entire area of the investigation 
is treated as a single reference area. 


tive temperature anomaly in the settled urban 
area, but also indicates to a certain extent an 
influence of topography on the temperature 
conditions. The largest positive temperature 
anomaly occurs, as we should expect in the 
central parts of the city. One degree isotherms 
generally follow the outer contours of the 
city but extend both north and south far out- 
side of the inhabited area; this probably re- 
sults from topographical conditions. The majo- 
rity of the observation runs were taken during” 
the winter or carly spring when temperature 
inversions are likely to occur. Therefore it is 
to be expected that points in the north and 
south, situated at higher levels, should have 
higher temperatures than the surrounding 


terrain. 


The isotherm map in fig. 4 were used 
in subdividing the general area for the 
studies of representativity. The more or 
less concentric lines around the center of 
the city roughly indicate the degree of in- 
fluence of the settled city-area upon the tem- 
peratures. The areas thus selected consist of 
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Godk. för publ. av Rikets allmänna kartverk. 


Fig. 7. Isotherm map for February 25, 1948, around 

0600 hours. Cloudiness: 0. Wind: SW 1.5 m/sec. 

Absolute humidity: 1.9 mm. Relative humidity: 91 %. 
Ground completely covered with snow. 


a central city area with three more or less 
concentric surrounding rings. It is obvious 
that the variations in size and form of the sev- 
eral areas somewhat affect the values for 6. 
In spite of this, such a division of the general 
area is preferable to the more artificial one, 
e.g. in form of a regular square net. 

Fig. 5 shows the representativity of the 
different station points. It is noticable that 
“all of the observation points near the central 
parts of the city show very small standard 
deviations. This results from the fact that tem- 
perature variations at all points inside the area 
occur in phase, a fact favoured to some ex- 
tent by the small dimensions of the area, but 
also clearly indicating that almost all central 
points are representative for their area. 

If one proceeds from the center of the city 
towards the surrounding countryside, points 
in the boundary zone between the city and the 
countryside show especially large deviations. 
The explanation is perhaps that the horizontal 
temperature gradient in this boundary zone is 
often especially large. A large temperature 


Godk. för publ. av Rikets allmänna kartverk. 


Fig. 8. Isotherm map för August 30, 1948, around 
2000 hours. Cloudiness: Acu, ci 0. Wind: ENE 
1.0 m/sec. Absolute humidity: 8.9 mm. Relative 


humidity: 88 %. 


gradient means that these observation points 
are sometimes found inside of the effective 
area of the city and sometimes outside of it, 
depending upon the prevailing wind direction 
and velocity. This promotes great variations 
in the value of 6 and of o. 

Outside of the effective area of the city 
proper one sees immediately the influence of 
the topography. Points having a typical loca- 
tion (high or low elevation) show large rela- 
tive standard deviations. 

A third factor can also be pointed out. In 
the wooded section south of the city high 
values are probably partly conditioned by 
the vegetation which gives this district a 
more maritime temperature climate than the 
adjacent terrain. 

The best points of observation in the country 
are found where neither topography nor vege- 
tation gives the area a special character. It is 
therefore completely consistent that the obser- 
vation points located on the plateau south-east 
of the city show the best representativity. 
The terrain is especially level, no mounds or 
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depressions occur, and the entire area forms es- 
sentially one homogeneous, cultivated surface. 

A final observation route was chosen based 
upon these results and is shown in fig. s. 
This route embraced a stretch of about 1 
Swedish mile (6 English miles) and included 
15 observation points. Some of the points 
chosen were representative of the urban area 
and some of the country. The effect of topo- 
graphy and vegetation was extremely slight 
in these points. 

Fig. 6 shows for comparison the result ob- 
tained when the entire area under investiga- 
tion is treated without subdivision. The most 
characteristic feature of the preceding figure 
also recurs here; the pronounced influence of 
topography and vegetation. A significant differ- 
ence remains however; the influence of the 
city is expressed in an entirely different man- 
ner. Points situated in the center of the city, 
certainly representative for the central part, 
are in no way representative for the larger 
area embracing the city and for the countryside. 
Accordingly observation points inside the 
city now show the greatest deviations. The 
unsuitableness of placing a meteorological 
station in a city follows naturally from this. 


Isotherm maps 


From a study of isotherm maps one can 
obtain a good understanding of the general 
temperature distribution and of the factors 
which determine it. Some of the most im- 
portant observations will be stated. 

The temperature in the settled area is in 
general higher than in its surroundings. Signi- 
ficant differences never arise during the day, 
but at night differences can often become very 
large. Positive city-countryside temperature 
differences are favoured by clear and calm 
weather. In such weather situations an espe- 
cially strong temperature gradient can often 
be established in the boundary zone between 
the populated and unpopulated areas. The in- 
habited area is enveloped by a limited, warm- 
air mass. The center of this warm-air mass is 
generally displaced somewhat in the direc- 
tion of the wind. The concentration of iso- 
therms, which often appears in those parts of 
the city outskirts which face the wind, can 
then be understood as a wind effect. 

The connection between topography and 
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temperature distribution can be shown very 
easily. The vertical temperature gradient is in 
this respect, the most important factor in 
determining of the size of topographically- 
conditioned differences. When the tempera- 
ture decreases with height, those points situa- 
ted at higher levels show somewhat lower 
temperatures than those located in the surround- 
ing area. The differences, however, are not 
large, usually only a few tenths of one degree, 
provided that the differences in level do not 
exceed ten meters or so. If the vertical tempe- 
rature distribution is reversed, however, i.e. 
if we have an inversion, the topographically- 
conditioned temperature differences may reach 
values of many degrees. 

It was possible to study the influence of vege- 
tation on the temperature distribution in the 
forest south of Uppsala. A persistant effect 
on the daily temperature trend was ascertained: 
the forest climate is more maritimely stressed. 
This observation is in general agreement with 
known conditions. Vegetation reduces outgo- 
ing nocturnal radiation and therefore prevents 
temperature inversions near covered surface 
from becoming as pronounced as those above 
open fields; as a result temperatures in the 
forest are higher during the night than 
those over the surrounding terrain. The situa- 
tion reverses itself, however, during the day. 
Incoming radiation. is obstructed by the vege- 
tation and the temperature can not rise to 
values as high as those observed over the open 
areas around. 


Temperature profiles 


In certain interesting or especially typical 
weather situations measurements were car- 
ried out at short time intervals along the short- 
er observation route. These temperature pro- 
files made it possible to study how the city- 
countryside temperature contrast is developed 
and how it changes with time. Two examples 
are shown in figs. 9 and 10; the text with 
these illustrations describes briefly the effect of 
the meteorological factors on the temperature 
contrasts. 


The relation between temperature contrast 
and simple climatological elements 


The isotherm charts and temperature pro- 
files clearly indicate that an intimate connec- 
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Fig. 9. Temperature profiles through Uppsala Octo- 


ber 28—29, 1948. In the lower part of the figure a 
raised profile shows the topographical conditions along 
the route (Cf. the continuous line in fig. 5). Built-up 
area is schematically shown by the sign for building. 


Profile Time Cloudiness Wind Humidity 
I I10O [6) NNW 1.3 3.4 (80) 
II IAIS o N 1.8 3.8 (73) 
Ill 1605 o N 1.6 3.9 (78) 
IV 1820 o NNW 1.5 3.5 (88) 
V 0040 I NNW 2.1 2.9 (89) 
VI 0405 10 N 222 Rat (OP) 


tion exists between the temperature differen- 
ces and the simple climatological elements; 
especially cloudiness and wind velocity. Parti- 
cularly in the analysis of the temperature pro- 
files it became apparent that small variations 
in cloudiness and wind velocity were also 
clearly reflected in the local temperature con- 
trasts. One may therefore suppose that there 
is a possibility of deriving a quantitative ex- 
pression from which the temperature diffe- 
rence between two arbitrary points could be 
computed. 

For a period of about 1 year, from May 
1948 until May 1949, it was possible to deter- 
mine more than 200 values of the tempera- 
ture contrast between city and country. The 
measurements were, as far as possible, dis- 
tributed equally over the seasons and time fo 
day. During the entire period every effort 
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Fig. 10. Temperature profiles through Uppsala Janu- 
ary 18—19, 1949. Ground completely covered with 
snow. 
Profile Time Cloudiness Wind Humidity 
i 1435 Otek, W 183 3.2 (92) 
II 1645 Stcu 8 W 2.6 3.I (94) 
Ill 2005, St TO NW 4.0 3.2, (02) 
IV O100 o SW 16 2.1 (89) 
V 0820 Acu 3 ENE 1.9 1.8 (91) 


was made to include all possible weather 


types. 
he mean value of three readings at the 
City Common with the reading at the nearb 
observation point situated south-east of the 
Common has been used as the value for the 
temperature in the centre of the city. The 
temperature of the open country has been 
defined as the mean value of the readings taken 
at the two most distant points of observation 
along the route. The readings have been re- 
duced to the same time in all cases. The diffe- 
rence between the mean temperatures com- 
puted in this way was taken as the observed 
temperature contrast between city and 
country. 

The simultaneous values for the simple 
climatological elements; cloudiness, wind ve- 
locity, temperature, as well as absolute and 
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relative humidity; were obtained partly by 
the investigator’s own observations and partly 
from recorded data taken at the Meteorolo- 
gical Institution. 

If in a diagram one represents the observed 
temperature contrast as a function of one of the 
simple climatological elements, e. g. wind 
velocity, one obtains a result having very 
great diffusion at the points. One can certainly 
discern a general tendency toward lower tem- 
perature contrasts with higher wind velocities, 
but such a technique cannot provide the foun- 
dation for a quantitative calculation. It is 
therefore impossible to isolate immediately a 
single factor and treat it separately. Instead one 
must attempt to treat in some way the whole 
complex of causes simultaneously. 

An attempt was made to solve this problem 
by calculating the partial correlation co- 
efficients for the relation between the 
temperature contrast and a particular meteoro- 
logical element. Thereafter the regression 
equations were derived in which the tem- 
perature contrast is represented as a function 
of the elements included in the calculation. 

Fundamentally, the material should be divid- 
ed into several groups with respect to season, 
time of the day, etc., and then each group 
should be subsequently treated individually. 
Such a subdivision is limited, however, by the 
fact that soon the numerical data in each sub- 
group becomes so small that it does not per- 
mit statistical analysis. In view of this difficulty 
it was decided to divide the material into two 
large groups; observations made during the 
day and observations made during the night. 
In this way the material is divided into those 
two groups which with respect to radiative 
phenomena can be expected to be more or 
less uniform. Since the number of observa- 
tions remaining in each group amounts to 
96 and 111 respectively, the material permits 
certain quantitative conclusions. In this paper 
only a short summary of the statistical result 
(without more detailed discussion) is given. 

The temperature contrast D between the 
central parts of the city and the open country 
is assumed to be expressed by the linear rela- 
tion: 


D=A+bN+bU+b:-0+b,e; (6) 


where N is cloudiness in the ten-degree scale, 
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U the wind velocity, © the temperature and e 
the absolute humidity. 

Insertion of the numerical values obtained 
from the statistical calculations gives the equa- 
tions the following form: 


Dpay = 1.4° — 0.01 N — 0.09 U — 


— 0.01 9 —0.04e (°C); (7) 
Dwight = 2.8° — 0.10 N — 0.38 U— 
— 0.02 8 + 0.03 e (°C); (8) 


The coefficients for the different meteoro- 
logical elements may not of course be consid- 
ered as certain in the final numerical values. 
The order of magnitude should however be 
essentially correct. The coefficients of correla- 
tion between Dop and Dea were calculated 
at + 0.49 + 0.05 and + 0.66 + 0.03 for day 
and night respectively. 

It is necessary to clearly realize that the corre- 
lation relations obtained are of a formal, ma- 
thematical nature, and that the statistical 
computations do not in themselves answer 
questions concerning the causal complex. 
They can however be of guidence in further 
study and of some use in practical investiga- 
tions. 

Summarily it can be said that equations 7 
and 8 approximately represent the tempera- 
ture difference between city and country in 
Uppsala expressed in simple climatological 
elements. The lack of good agreement between 
the calculated and observed differences, espe- 
cially at high and low values, seems to in- 
dicate that the assumption of a linear relation 
between the wind velocity and the tempera- 
ture difference is incorrect. No appreciable 
variation in the repreSentativity of the equa- 
tions was observed during different seasons of 
the year. Equation 7 appears to give values 
for the temperature difference which are some- 
what too low near the times of sunrise and 
sunset. Very small or negative values for the 
difference were observed only in relatively 
high solstice. Very large differences were as- 
sociated with situations having well developed 
temperature inversions. 


The effect of the growth of the inhabited 


community 


The foregoing has sought to indicate that 
local city climate often exhibits great tem- 
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Fig. 11. The daily temperature amplitude in Stock- 


holm (broken line) and Uppsala (continuous line) for 
January and July. Mean numbers for decennial periods 
from 1871 to 1940. 


perature contrasts between city and the sur- 
rounding countryside. In a rapidly growing 
city the influence of the inhabited area upon 
a place of observation in.the city can be 
expected to increase continually. An interest- 
ing example of this effect can be seen in a 
comparison between the temperature values 
from Stockholm and from Uppsala. 
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In fig. 11 decennial mean values for the daily 
temperature amplitude in Stockholm and in 
Uppsala are given for the months of January 
and July over the period 1870—1940. The 
amplitudes were calculated as the difference 
between the daily maximum and minimum 
temperatures. Uppsala’s temperature ampli- 
tude shows a relatively horizontal course. 
One can perhaps detect a very slight tendency 
towards decreasing amplitude. The condition 
in Stockholm is otherwise; for both months 
one notes a significant and continual decrease 
of the temperature amplitudes. This diminu- 
tion for the month of July goes to more than 
2°C in comparison with the Uppsala values. 
It can be stated parenthetically that the changes 
in amplitudes were produced both by an 
increase in the minimum temperature and by a 
decrease in the maximum temperature. The 
cause for these changes is difficult to explain; 
possibly it is related to changes in the vegeta- 
tion around the place of observation, but it is 
also probable that the development in Stock- 
holm of the populated area itself plays an 
important röle. 

It then follows that even observations from 
very good stations must, under such condi- 
tions, be most carefully interpreted. Especially 
when one wishes to discuss secular climate 
changes it becomes particularly important 
that the influence of a developing city or 
changes in the local surroundings of the obser- 
vation point be not interpreted as real climatic 
changes. 
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Geostrophic Departures in the Jet Stream 


By VICTOR P. STARR, Massachusetts Institute of Technology! 


(Manuscript received 2 April 1950) 


Abstract 


The effect of gradients in the cross-current eddy motions on the geostrophic balance 
is illustrated by means of an integrated model. It is suggested that gradients of the 
Reynolds stress associated with such eddy motions in the atmosphere and in the oceans 
may lead to sensible average geostrophic departures which appear as a normal inertial 


effect in eddying currents. 


In the past many investigators have devoted atten- 
tion to the study of the degree to which actual 
wind currents may be measured by the geostrophic 
approximation. Some of these studies have been 
directed toward an examination of effects which 
might be expected on rational grounds, while 
others are entirely empirical in nature. In view of the 
fact that research workers interested in the large- 
scale currents in the atmosphere are usually forced 
to resort to wind estimates obtained from pressure 
data, the nature of the approximations involved in 
this procedure is of much practical concern. The 
purpose of this discourse is to bring attention to 
one systematic effect which might be expected on a 
rational basis in those regions of the atmosphere 
where the circulations attain their greatest intensities, 
although lesser manifestations of the same kind 
might be found to exist at other levels in the at- 
mosphere as well as in the oceans. 

The effect to be discussed cannot be regarded as 
any new concept in the field of hydrodynamics, 
having long ago been recognized in the treatment 
of such topics as the turbulent flow of fluids in 
channels (see for example GoLDSTEIN 1938). How- 
ever, it appears that the meteorological applications 
of the concepts involved should be emphasized. 

Although the subject is not dependent upon the 


1 This research was performed under contract with 
the U. S. Air Force Cambridge Research Laboratories. 


existence of any particular integrated model of flow, 
we shall find it convenient to introduce the basic 
arguments through the examination of one such 
specific model. In terms of cartesian coordinates in 
which x positive is taken eastward and y positive 
northward, let us study the two-dimensional flow 
of a fluid which satisfies the following equations of 
motion and continuity: 


Ou, du ou I dp 

ru v— = fv — — — 
dt ox dy 0 0x 
dv dv dv I Op 

| ) == 
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Here u is the eastward and v the northward particle 
velocity, p is pressure, f time, f the Coriolis para- 
meter and © density. We assume that @ is constant 
and uniform and that 


Key (2) 
where fo and B are constants. 


The three relations (1) form a closed system of 
partial differential equations for the determination 
of the three dependent variables u, v, p in terms of 
the independent variables x, y, f. It has been shown 
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by NEAMTAN (1946) that this system of simultaneous 
equations possesses an integral in which 


— Al sin k (x—ct) : cos Î y, 
(3) 


+ sin | y, (4) 


u = U + V cos (xy+e) 


v = Ak cos k (x—«t) 
provided that 
eye EAN (5 


and U, V, A, k, I, e are arbitrary constants which we 


take as being real. 
This (particular) solution is completed by writing 


that 


Alc sin k (x —ct) : cos ly 


I ; : 
— — Aa? sin? k (x — ct): sin? ly + 
2 


+ À sin k(x— ct): sin ly [fo + By + 
+ Vasin (xy +e)] 


V . 
— fo ar (fo +ßy) sin (ay + le 


A v(4-2) THON 


a 


I 
SE Pr 
2 


. (xy +6) 
+ F (6), (6) 


obtained through the integration of the first two 
equations of (1) with the aid of (3), (4) and (s). 
The last term in (6) is an arbitrary function of time 
alone and may ordinarily be taken to be an appro- 
priate constant without affecting the other charac- 
teristics of the solution. In view of the form of the 
continuity equation contained in (1), it is possible 
to use a stream function y for the velocity compo- 
nents, of the form 


7 


y=— Uy—— sin (xy + e)+Asink (x—ct) : sin ly. 
& ‘ 


(7) 


It is apparent that the motion consists of a basic 
zonal flow on which periodic disturbances are 
superposed. 

Letting u, represent the geostrophic wind toward 
the east, we write that 


uqg=uUu—u, | (8) 
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so that ug is the nongeostrophic part of u. Also, 
introducing the space average of a quantity ( ) 
over a wave length L in the: eastward direction by 


the relation 


——— I 


i ) ax, (9) 


we have that, 


iy == DP The 


(10) 


Either by using the geostrophic formula with (6) 
and comparing with (3), or more directly by sub- 
stitution in the second equation of (1) from (3) 
and (4) and averaging, we find that 


= A2k?] 
Ua = — 
A 7 


sin ly + cos ly. (11) 


It is thus seen that # contains a systematic nongeo- 
strophic component which varies in a periodic . 
fashion with the independent variable y alone. 

Since v vanishes for y = O and for y = x/l, 
it will tend to clearness if we visualize the flow as 
contained in a channel between two parallel walls at 


these ordinates. The departure ug then vanishes 
according to (11) at the two walls and also at the 
= 7/2 |. In the southern half of the channel, 
between y = O and y = 2/21, ug is negative 
while in the northern half it is positive. On the 
average over the whole channel the conditions are 
such that the flow is geostrophic. 

By way of physical explanation of this phenome- 
non we might note that the minimum value of the 


center, y 


Reynolds stress — ovv is found at y = 2/2 1. This is 
interpreted as a net transfer of positive y-momen- 
tum from the southern half to the northern half. 
This drain must be made good in the lower half by 
an excess of the mean northward pressure gradient 
force over and above the requirements needed to 
balance meridional Coriolis forces on the fluid in 
this region. The reverse process takes place in the 
upper half. 

Viewed somewhat differently the process may be 
explained physically by noting that individual par- 
ticles enter the southern half with negative y- 
momentum and leave with positive y-momentum. 
They must therefore be accelerating northward on 
the average while in the southern half and south- 
ward while in the northern half. From this one 


may infer that ug cannot be zero everywhere, 
unless indeed the motion is purely zonal. 
It is to be noted that the results described do not 
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depend upon the profile of u but only upon that of 
ovv. We may, for example, place a maximum or 


»je of u in the southern half at y =27/4 1 by taking 
an appropriate value U together with V > o and 
proper values of &, e. 

Much of what has been said follows directly if 
we write the second equation of (1) in the form of a 
continuity equation for y-momentum with the aid 
of the third and average so that 


dv 


AT Jovy _ Eu 44 
dt dx dy 3 dy 
The results stated above follow immediately, since 
the first two terms vanish. Moreover we see that 
similar results are to be expected in any two dimen- 
sional model with à quasi-steady regime of turbu- 
lence or eddy motion. For all such cases the mean 
geostrophic deviation is determined by the cross- 


stream gradient of the Reynolds stress — ovv. 
The direct application of the model described to 

the atmosphere would lead to the presence of a 

westward directed mean geostrophic departure 


equatorward from the latitude at which ovv is a 
maximum in each hemisphere, and an eastward 
mean departure poleward. However, the effect is 
in actuality superposed on other sources of geo- 
strophic departures such as the curvature of the 
latitude circles which produces a westward depar- 
ture even in the case of purely zonal flow. Another 
factor is that for the actual atmosphere even in the 
simple cartesian equation (12) an additional verti- 
cal term must be added so that, taking z positive 
upward and w as the upward particle velocity, and 
neglecting horizontal density variations in the:x- 
direction, 


= 1 (dovv , dovw 
HT — — x 
Sd Pe), 


It is difficult to estimate the effect of this vertical 
term, although if (13) be integrated also with 
respect to height through the whole atmosphere, 
the total contribution of this added term must re- 
sult from the presence of a net meridional surface 
stress at a given latitude. There is not much a priori 
reason to expect that such net meridional surface 
stresses around latitude circles are of importance. 


(13) 


(7 
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The latitude at which ovv reaches a maximum 


is probably displaced farther poleward in each he- 


misphere as compared with the latitude at which u 
reaches its maximum, although observational stu- 
dies should be made concerning this point. Further- 
more, there is no reason to suppose that this quan- 
tity tends to zero as the pole is approached. For 


purposes of orientation let it be supposed that vv 
increases by 107 cm?/sec? for an increase of 2 000 
km in y. Neglecting horizontal density variations, 
the first member on the right hand side of (13) 
then gives a contribution of— 6.7 meters per se- 
cond to ug, if f = 0.75 x 10”? (corresponding to 
about 30° N latitude). This may give some notion 
of correct orders of magnitude in the vicinity of 
the jet stream. At lower levels the contribution 
may be one fifth or one tenth as large. 

In an interesting paper published recently, LOEwE 
and RADOK (1950) have presented a meridional cross 
section of the atmosphere in the southern hemisphere 
through the jet stream. These writers make a com- 
parison of geostrophic versus actual wind measure- 
ments in the jet. The actual winds show velocities 
much lower than those computed from the iso- 
baric contours. Although many questions might 
arise concerning the representativeness of the data, 
it may also be that a part of the discrepancy is due 
to the manifestation of the phenomenon here 
described. It is to be noted that local climatolog- 
ical time averages such as those used by LOEWE 
and RapoK would tend to approximate space 
averages except for the possible presence of stand- 
ing semipermanent disturbances. 

Finally it should be mentioned that the considera- 
tions used here may be of some consequence in the 
charting of mean ocean currents by means of dy- 
namic velocity computations, in those cases where 
there exists a gradient in the cross-current eddy 
motion. This may, for example, be significant in 
the case of such currents as the Gulf Stream. 
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Letter to the Editor 


Late Cenozoic Climatic Changes as Recorded by the Equatorial 


Current System. 


Dear Sır: 


Since the preparation of the above indicated paper 
(Tellus, 2, 2) additional petrographical and palaeonto- 
logical information has become available which gives 
a still more detailed picture of the marked stra- 
tification mentioned therein. With regard to one 
detail the statements of the previous paper must 
however be changed. On page 85, column 2, 
line 9 from the bottom it is stated as an added 
remark, that at 5° S the biogenous sedimentation 
increases. Control calculations and increased sample 
frequency reveal however that on the contrary 
the biogenous sedimentation slightly decreases at 


this latitude, the average number of foraminifera 
per gram of the sediment being approximately 
1200 as compared with approximately 1500 at the 
Equator. The suggestion of a secondary conver- 
gence at 5° S is therefore unnecessary. This modi- 
fication does not influence other statements and 
conclusions made in the article. 


Sincerely yours, 


Gustaf Arrhenius 


Kagghamra, Grödinge, 
September sth, 1950 


